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Presentazione 


Il 13 ottobre 1989 ha avuto luogo a Torino, nella sala di Consiglio 
di Facoltà di Ingegneria del Politecnico, la Giornata di studio a ricordo 
di Gianni Jarre, improvvisamente ed immaturamente scomparso l’8 ot¬ 
tobre 1988. 

Erano presenti numerosi colleghi ed amici, che hanno voluto manife¬ 
stare la stima, l’affetto ed il desiderio di mantenere vivo il ricordo di 
Lui, all’insegna di quelle Scienze Meccaniche, in cui Egli aveva prodi¬ 
gato le doti migliori del Suo intelletto. 

A quell’appuntamento ero presente anche io, che fin dal lontano 1948 
avevo con Lui iniziato a percorrere l’arco della carriera universitaria, 
prima come assistenti di Aerodinamica, sotto la guida del nostro comu¬ 
ne, sommo Maestro Carlo Ferrari, poi come docenti in rami diversi del¬ 
la Meccanica. 

Mi è stato affidato il compito di presentare questo volume che racco¬ 
glie i contributi della Giornata di studio: che profonda commozione mi 
rivolgo ai gentili lettori, nel segno della grande amicizia che mi legava 
a Gianni Jarre, e della grande stima che riponevo nella Sua acuta ed ar¬ 
guta intelligenza. La Sua figura di scienziato,, di docente, di uomo, è 
illustrata nel modo più autorevole nella commemorazione di Carlo 
Ferrari. 

Testimonianze importanti e significative vengono date da Rolando Ri- 
gamonti, Pietro Buzano, Carla Oggiano e Vittorio Marchis. 

Segue poi la parte più strettamente attinente al tema specifico della 
Giornata, con le memorie di Massimo Germano, Ettore Antona, Gino 
Moretti, Maurizio Pandolfi, Luigi G. Napolitano e Antonio Viviani, Pla¬ 
cido Cicala, Ario Romiti, oltre che del sottoscritto. Tutti devono essere 
ringraziati caldamente per essersi prodigati con impegno e tempestività 
in questo tributo corale di stima pe l’Amico scomparso, ma soprattutto 
per averlo fatto con una forte carica spirituale ed affettiva, a testimo¬ 
nianza del sincero rimpianto per il nostro indimenticabile Gianni. 

Torino, Aprile 1991 


Silvio Nocilla 




Commemorazione di Gianni Jarre 

(1924-1988) 


Carlo FERRARI 


Commemorare un caro Collega è sempre triste compito, ma se que¬ 
sto Collega è stato Allievo, Collaboratore intelligente ed attivo dalPini- 
zio della Sua carriera accademica sino alla fine, amico carissimo sempre 
vicino nei momenti tristi e lieti della vita, alla tristezza si aggiunge un 
sentimento più profondo, più intimo, è come se qualcosa dentro di noi 
si sia spezzato. Comprederete perciò il mio stato d’animo. 

Io ho conosciuto Gianni Jarre nel 1948; mi fu presentato da Gino 
Moretti, che era mio Assistente per la Cattedra di Meccanica razionale 
e aveva deciso di lasciare Torino per correre la sua avventura in Argen¬ 
tina. Moretti e Jarre si erano conosciuti a Vevey (Svizzera), dove erano 
espatriati per sfuggire ai nazi-fascisti. 

Jarre era a quell’epoca allievo in questo Politecnico, e avendo potuto 
seguire a Vevey Corsi tenuti da illustri professori italiani pure rifugiatisi 
in Svizzera, fra i quali Gustavo Colonnetti, al Suo ritorno in Italia, alla 
fine della guerra, era riuscito a conseguire la laurea in Ingegneria Mec¬ 
canica già nel 1948. Il risultato di quel colloquio fu che, quando Moret¬ 
ti ci lasciò, Gianni gli subentrò come Assistente, e da quell’epoca non 
mi abbandonò più, seguendomi quando io, dopo il ritiro del Prof. Pa¬ 
netti, passai all’insegnamento della Meccanica applicata. La Sua bril¬ 
lante intelligenza e la Sua profonda preparazione lo spinsero subito a 
svolgere una intensa attività scientifica, e con Lui e con altri giovani stu¬ 
diosi di altissimo valore (giovani allora, valentissimi sempre), Silvio No- 
cilla, Aldo Muggia, Ario Romiti, si costituì un gruppo animato dal più 
grande fervore e desiderio di affrontare i problemi alla ribalta dopo il 
digiuno scientifico prodotto dai lunghi, dolorosi anni della guerra. Tale 
attività lo portò nel 1955 al conseguimento della libera docenza in Ae¬ 
rodinamica, e nel 1958 alla cattedra prima di Aerodinamica e poi di 
Gasdinamica presso il Politecnico di Torino. Qui divenne direttore dal 
1973 al 1978 dell’Istituto interdisciplinare di Meccanica applicata, Ae¬ 
rodinamica e Gasdinamica, e dal 1974 direttore del Centro Studi di Di¬ 
namica dei fluidi del C.N.R. Fu inoltre Presidente del Consiglio scientifico 
del «Centro Studi per le ricerche sulla Propulsione e sull’Energia» del 
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Politecnico di Milano; Membro del Consiglio scientifico dell’INSEAN 
(Vasca Navale) di Roma, e del CETENA (Centro Studi Tecnica Navale) 
di Genova. Per la sua rinomanza di Tecnico e Scienzato insigne fu ricer¬ 
cato per consulenza da industrie privati ed Enti pubblici nazionali ed 
internazionali: così, fu Impiantista industriale termomeccanico dal 1948 
al 1960, Consulente della C.E.E. nel 1968-70 per un’indagine sull’indu¬ 
stria aerospaziale europea. Nel 1962 fu accolto per un periodo di due 
mesi come «visitatore» dall’Accademia sovietica delle Scienze, dove fe¬ 
ce una proficua indagine su svariati problemi di Aero e Gasdinamica. 
L’intenso impegno nell’attività scientifica non lo distolse da un intenso 
impegno anche nell’attività didattica: tenne, oltre all’insegnamento del 
Corso di cui era titolare, quello di «Meccanica applicata alle Macchine 
e Macchine» per allievi-ingegneri civili dal 1960 al 1978; profondamen¬ 
te sensibile ai problemi della Scuola, e non soltanto di quella universita¬ 
ria, svolse per diversi anni corsi di aggiornamento tecnico-scientifico per 
i Docenti di formazione professionale, tra i quali corsi citerò uno dedi¬ 
cato alla Matematica di base, ed un altro, in collaborazione con Gian¬ 
franco Chiocchia, alla Meccanica applicata. 

Per l’opera da Lui svolta nei campi tecnico, scientifico e didattico non 
gli mancarono i riconoscimenti: fu nominato Socio corrispondente del¬ 
l’Accademia delle Scienze di Torino nel 1960 e Socio Nazionale della 
stessa Accademia nel 1983; fu insignito della «Medaglia d’oro di bene¬ 
merito della Scienza, della Cultura e dell’Arte» da parte del Ministero 
del P.I. nel 1984; era Membro del Consiglio direttivo del Centro di Astro- 
dinamica «Giuseppe Colombo». Altri riconoscimenti avrebbe certamente 
meritati. 

Per quanto si riferisce alla Sua attività scientifica, questa riguarda es¬ 
senzialmente tre campi: Macchine a fluido; Gastermochimica; Flussi di 
strato limite; e le sue caratteristiche principali sono: ingegnosità e, al¬ 
meno per i primi due campi sopra menzionati, pionerismo. 

Le ricerche sulle Macchine a fluido furono le prime da Lui effettuate 
nel periodo 1948-1952: alla fine della seconda guerra mondiale erano 
apparsi i primi aeroplani militari a reazione, come si diceva impropria¬ 
mente allora, a getto secondo la denominazione attuale. 

La propulsione a getto era problema da lungo tempo studiato, ed era 
stata discussa pure nel celebre «Convegno Volta» del 1935, presso l’Ac¬ 
cademia d’Italia, ma la sua realizzazione richiedeva la possibilità di co¬ 
struire compressori rotativi di elevato rendimento: da qui la necessità 
di uno studio delle leggi che determinano il moto di un gas in una turbo- 
macchina, più approfondito di quello corrispondente agli schemi tra¬ 
dizionali di calcolo delle macchine stesse, allora in uso. Detti metodi 
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tradizionali furono molto duri a morire, particolarmente in Italia: ora 
Jarre nel 1948 (e quindi a soli 24 anni) nella nota «Un nuovo tipo di 
stadio di compressione, per compressori assiali» studia il tipo di com¬ 
pressore detto «a vortice libero», per il quale la circuitazione è costante 
lungo le pale, e per detto tipo indica un procedimento di calcolo per il 
dimensionamento del compressore stesso, dimostrando come la adozio¬ 
ne dello stadio proposto consenta una notevole riduzione del numero 
totale di elementi, rispetto ai normali compressori assiali multipli. Nel 
1952 poi scrive una monografia «Compressori assiali a vortice libero», 
che è in realtà un trattato sull’argomento, che per la completezza e la 
profondità, colla quale i problemi aerodinamici e termodinamici ad es¬ 
so relativi sono considerati, è tuttora pienamente valido. In effetto, in 
detta monografia, tra altri risultati interessanti: 1) è dedotta e discussa 
la dipendenza del rendimento del compressore dall’efficienza dei profili 
alari utilizzati come sezioni delle pale; dai gradi di azione e di reazione 
della sezione periferica, definiti come rapporti del salto di energia rea¬ 
lizzato nello statore e nel rotore, al salto totale di ogni singolo stadio; 
dal rapporto di funzionamento, definito come rapporto tra la velocità 
assiale e la velocità rotatoria periferica del rotore, ed infine dal rappor¬ 
to fra raggio interno e raggio esterno di ogni stadio; 2) è analizzata l’in¬ 
fluenza della compressibilità, e quindi del numero di Mach, e della 
viscosità, e quindi del numero di Reynolds, sulle prestazioni del com¬ 
pressore. Sempre al campo delle macchine a fluido appartiene un’altra 
Nota del titolo: «Moto di un fluido compressibile in una girante radiale 
di turbomacchine» (1950), particolarmente notevole per l’eleganza del¬ 
la trattazione matematica: il problema era stato studiato da Busemann 
nel 1928, per il caso del fluido incompressibile non viscoso, per il quale 
la soluzione può essere ottenuta in modo relativamente semplice col me¬ 
todo della trasformazione conforme, inapplicabile però nel caso ora in¬ 
dicato. 

Jarre linearizza le equazioni del moto colla semplice ipotesi che l’au¬ 
mento di entalpia sia essenzialmente prodotto dal lavoro della forza cen¬ 
trifuga, e riesce ad esprimere, almeno in prima approssimazione, il 
potenziale del flusso compressibile attraverso la girante per mezzo del 
potenziale del flusso incompressibile. 

Pioneristiche ed ingegnose sono pure le ricerche svolte nel periodo 
1955-1965 sulla gasdinamica delle miscele di gas inerti o reagenti; e co¬ 
me quelle precedenti sulle macchine a fluido furono determinate dall’av¬ 
vento della propulsione a getto, così queste ricevettero stimolo ed impulso 
dal grande interesse suscitato dai nuovi fenomeni, che il raggiungimen¬ 
to delle alte velocità consentito da detto sistema di propulsione propose 
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alla stadio dei fluidodinamici. Questi fenomeni sono quelli connessi colle 
elevate temperature, che si generano in una corrente ad alto numero di 
Mach quando essa è bruscamente frenata o da un’onda d’urto o dall’at¬ 
trito a contatto di una parete solida nello strato limite adiacente ad es¬ 
sa. Ed anche in questo campo Jarre sviluppò in una magistrale 
monografia «Elementi di gasdinamica delle miscele» (1955) ì fondamenti 
di questa nuova branca della Dinamica dei fluidi, dando ad essi notevo¬ 
li contributi originali, formulando differenti espressioni della equazio¬ 
ne dell’energia, delle quali diede interessanti interpretazioni, ed indicando 
come le ipotesi della termodinamica dei processi irreversibili, ed in par¬ 
ticolare il teorema di Onsager, siano implicitamente contenute nella teoria 
cinetica dei gas. In una Nota successiva «Un modello semplificato di mi¬ 
scela multipla di gas» (1956), pubblicata dall’Accademia dei Lincei, nelle 
ipotesi che i coefficienti di termo-diffusione siano nulli, i coefficienti di 
mutua diffusione uguali, i numeri di Prandtl e di Schmidt siano uguali 
a uno per tutti i costituenti la miscela, dimostrò che il trasporto di mas¬ 
sa del costituente i-esimo, a parte la sua produzione chimica di massa, 
è prodotto dal gradiente di frazione di massa soltanto di detto costituente, 
mentre, nel caso di flusso di strato limite, il trasporto conduttivo, diffu¬ 
sivo e viscoso di energia è ridotto ad un fenomeno di opposizione al so¬ 
lo gradiente di energia totale (ossia, termica, chimica e cinetica), ed infine 
che, sempre per fenomeni di strato limite esiste l’integrale primo dell’e¬ 
nergia totale: questo, che generalizza per miscele multiple di gas reagen¬ 
ti il classico integrale primo di Busemann, può ben a ragione chiamarsi 
«integrale primo di Jarre». Sempre nell’ipotesi già detta relativa ai nu¬ 
meri di Prandtl e di Schmidt, considerando una corrente uniforme di 
miscela di gas e di vapore surriscaldato lambente una superficie piana 
parallela alla corrente stessa e sulla quale il vapore può cambiare stato, 
nel caso in cui siano costanti le velocità tangenziali, le entalpie e le con¬ 
centrazioni di vapore sia nella corrente libera che sulla superficie, Jarre 
dimostrò resistenza di altri integrali primi esprimenti l’analogia dei campi 
(di velocità, di entalpia, di concentrazione), per mezzo della quale le leggi 
di variazione della velocità, della entalpia e della concentrazione sono 
definite da una unica funzione, e l’analogia degli scambi (di quantità 
di moto, di calore e di massa), che mette in relazione la tensione di attri¬ 
to alla parete, la quantità di calore trasmessa per unità di superficie nel¬ 
l’unità di tempo, e la portata di vapore emessa dalla parete per unità 
di superficie: l’analogia dei campi sussiste fra le grandezze misurate per 
unità di massa di miscela, mentre l’analogia degli scambi sussiste fra le 
grandezze misurate per unità di massa di gas. Dette analogie, che costi¬ 
tuiscono la generalizzazione delle classiche analogie di Reynolds e di 
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Colburn al caso di miscele binarie, in quanto esse valgono per velocità 
comunque elevate, per concentrazione di vapore comunque grande, per 
valori diversi dei calori specifici del gas e del vapore, e per tutti i regimi 
di moto, dall’ordinatissimo moto laminare di Couette ad disordinatissi¬ 
mo moto completamente turbolento, possono a buon diritto essere chia¬ 
mate «analogie di Jarre». Una applicazione molto interessante di esse 
è stata da Lui fatta nello studio del raffreddamento evaporativo per con¬ 
trastare il riscaldamento aerodinamico e superare la cosiddetta «barrie¬ 
ra termica», studio i cui risultati sono da Jarre sintetizzati con questa 
immagine: un mezzo fornito di un sistema di raffreddamento evapora¬ 
tivo ed un mezzo fornito di un sistema di iniezione di gas normalmente 
alla superficie lambita si comportano rispettivamente come bulbo umi¬ 
do e come bulbo secco di uno psicrometro atmosferico in una atmosfe¬ 
ra riscaldata per attrito. 

Altre fondamentali applicazioni della gasdinamica di miscele di gas, 
da Lui fatte, riguardano la determinazione della struttura di un’onda 
d’urto retta intensa in un gas puro biatomico, pubblicata nel numero 
del «Journal of Applied Mathematics and Physics» dedicato a Jacob 
Ackeret in occasione del suo 60mo anno, e del flusso nello strato limite 
di una parete adiabatica in corrente ipersonica, sempre di un gas puro 
biatomico. In entrambi i casi si presenta il problema della formulazione 
della velocità di reazione corrispondente ai processi di dissociazione e 
di ricombinazione conseguenti alle alte temperature provocate dalla brusca 
decelerazione della corrente ipersonica; difficilissimo problema di cine¬ 
tica chimica per il quale non esiste tuttora una soluzione generale soddi¬ 
sfacente. Tuttavia per la particolare semplice reazione sopra indicata, 
Jarre applicando la teoria della cinetica chimica assoluta, che postula 
1’esistenza di uno stato intermedio fra la molecola e i due atomi, chia¬ 
mato complesso attivato, le cui proprietà determinano la velocità di rea¬ 
zione, perviene ad una formulazione di detta velocità ammettendo che 
il complesso attivato abbia le caratteristiche che minimizzano la sorgen¬ 
te di entropia della reazione. In conseguenza del valore finito dalle va- 
locità della reazione chimica, sia a valle dell’onda d’urto sia a valle 
dell’arresto sul bordo di attacco della lamina piana si ha uno stato di 
non equilibrio, e quindi un processo di rilassamento durante il quale il 
fluido si raffredda progressivamente col progredire della reazione verso 
l’equilibrio, e i risultati delle ricerche cui prima ho accennato si possono 
così sintetizzare: a valle dell’urto il fenomeno di rilassamento è molto 
veloce, e può pertanto essere considerato come raggiuntimento quasi 
istantaneo dell’equilibrio; a contatto della lamina adiabatica il rilassa¬ 
mento è molto più lento e può pertanto in generale essere considerato 
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in prima approssimazione come persistenza indefinita dello stato di non 
equilibrio, ossia di flusso congelato, secondo la terminologia usata in 
gasdinamica. 

Ed ancora deve essere ricordata la fondamentale Memoria «The in¬ 
ternai structure of thè shock wave — Multicomponent mixtures», nella 
quale dimostra che lo spessore dello strato d’urto nel flusso di una mi¬ 
scela a N componenti inerti, qualunque sia il numero di detti compo¬ 
nenti, non è mai più grande di quello corrispondente a una miscela binaria 
dei soli gas più leggero e più pesante, e questo in contraddizione a quan¬ 
to si deriverebbe dall’osservazione che ogni nuovo fenomeno irreversi¬ 
bile, come quello della diffusione di un nuovo componente della miscela, 
dando origini ad un nuovo termine della sorgente di entropia, dovrebbe 
pure aumentare lo spessore dello «strato d’urto». 

Un terzo gruppo di ricerche fu dedicato da Jarre, come già prima detto, 
a partire dal 1966 allo studio aerodinamico e termodinamico del flusso 
laminare nello strato limite bidimensionale in condizioni di grande ge¬ 
neralità per quanto riguarda sia la corrente esterna allo strato, sia le con¬ 
dizioni alla parete delimitante lo strato stesso. Sono problemi questi che 
già erano stati trattati da illustri Cultori della Meccanica dei fluidi, qua¬ 
li ad es. Blasius, von Karman, Betz, Gòrtler, Lighthill, Schlichting; per 
altro i procedimenti da loro proposti o sono molto complessi, e questa 
complessità rende difficile la valutazione dell’approssimazione nella quale 
i vari problemi sono risolti, o non sono validi in presenza di singolarità 
del flusso, quali ad es. si hanno se si produce distacco della corrente dalla 
parete lambita. L’uso del calcolatore, colla possibilità che questo dà di 
ottenere soluzioni numeriche anche di equazioni di grande complessità, 
può far sembrare le analisi teoriche in questo campo superate, ma que¬ 
sto non è. Il calcolatore dà solo risultati numerici in ogni singolo pro¬ 
blema, e pertanto non è facile da essi dedurre proprietà di carattere 
generale; indica come un fenomeno avviene, ma non perché avviene, e 
la presenza di singolarità insieme colla non linearità delle equazioni ren¬ 
de il loro uso molto delicato. Mi ricordo di avere inteso da Giuseppe 
Colombo, che di Jarre fu grande amico e stimatore, che Egli era solito 
suggerire ai suoi Allievi e Collaboratori che prima di fare ricorso al cal¬ 
colatore, prendessero regolo e matita e vedessero se con questi fosse pos¬ 
sibile dedurre qualcosa che aiutasse a comprendere il problema fisico. 
E un concetto analogo mi pare esprima Gino Moretti quando in «To- 
wards a closer cooperation between Theoretical and Numerical Analy- 
sis in Gasdynamics» suggerisce di sottoporre ad una analisi teorica i 
fenomeni descritti dalle equazioni prima di procedere alla loro discre- 
tizzazione, e quindi all’uso del calcolatore. Ora le ricerche di Jarre, cui ho 
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accennato, si prestano mirabilmente a detta investigazione perché in es¬ 
se, con procedimenti tanto semplici quanto ingegnosi, si perviene ad il¬ 
lustrare chiaramente i fenomeni che avvengono nei flussi di strato limite. 
Io mi limiterò a segnalare tra esse solo quelle relative alla determinazio¬ 
ne delle tensioni tangenziali e della trasmissione termica alla parete per¬ 
ché particolarmente significative ai fini sopra indicati. Esse conducono, 
per mezzo di operazioni non lineari molto semplici, ad una determina¬ 
zione della tensione di attrito e della trasmissione termica alla parete in 
ogni punto di questa, che mostra in maniera molto espressiva la loro 
dipendenza dalle leggi di variazione della velocità e della entalpia ester¬ 
ne allo strato, e della entalpia alla parete in ogni punto di questa, e pre¬ 
cisamente, secondo le parole usate da Jarre, dal presente e dal futuro 
rispetto al punto in esame, dell’evolversi della velocità e della entalpia. 
In particolare, per la trasmissione termica, nel caso di entalpia di parete 
variabile, il metodo di Jarre mette bene in evidenza il sistematico dupli¬ 
ce meccanismo di ispessimento e di anticipazione, ispessimento dello strato 
limite lungo la parete, anticipazione nel senso che la trasmissione termi¬ 
ca in un punto generico avviene come se il valore della entalpia in detto 
punto fosse uguale all’entalpia in un punto spostato più a valle, e que¬ 
sta anticipazione risulta proporzionale ad una lunghezza, che memoriz¬ 
za quanto avviene a monte del punto stesso. L’applicazione del medoto 
di Jarre sia alle correnti «simili» rappresentate dalle equazioni di Falk- 
ner e Skan, e di Fage-Falkan, sia al flusso in corrispondenza del punto 
di arresto della corrente e sia alla corrente attorno ad un profilo circola¬ 
re, mostra un’ottima concordanza coi risultati ottenuti con soluzione nu¬ 
merica da Flartree e da Curie e Rosenhead ed inoltre, con detto metodo 
si trova subito la singolarità di Goldstein nel punto di separazione della 
corrente, mentre si toglie il paradosso di flusso termico nullo in questo 
punto, quale risultava invece dalle trattazioni di Lighthill e di altri. Ed 
ancora: a) si determina subito il punto di inversione termica, caratteriz¬ 
zato dall’annullarsi della trasmissione termica in detto punto e si dimo¬ 
stra che in esso non si ha alcuna singolarità; b) si ottengono le espressioni 
asintotiche dei coefficienti di trasmissione termica per numero di Prandtl 
tendente a zero (caso dei plasma, dei metalli liquidi), o tendente all’infi¬ 
nito (caso degli olii). 

L’attività scientifica di Jarre non si è limitata ai campi di cui ho fatto 
cenno; i Suoi interessi scientifici e la Sua cultura erano molto vasti. Aveva 
una mente termodinamica, ed aveva così bene assimilati i concetti, spesso 
astrusi, della termodinamica (si pensi al paradosso di Gibs relativo alla 
entropia di miscele di gas), che sapeva farli comprendere anche a chi 
aveva con essi poca familiarità. È stato detto da taluni che Egli non era 
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amico della Matematica: nulla di più falso. Aveva un grande rispetto 
della Matematica, ed un’ottima cultura di essa, come appare chiaramente 
dai Suoi lavori nei quali aveva saputo applicare con assoluto rigore e 
con eleganza metodi e concetti matematici molto elevati. Quello che non 
amava era il narcisismo della Matematica, o per dirla ancora colle paro¬ 
le di Giuseppe Colombo, non amava quella Matematica «che si limita 
a un raffinamento di vecchi strumenti di calcolo o lo sviluppo di forma¬ 
lismi, sovrastrutture e generalizzazioni non necessarie». 

Accanto alla figura di Gianni Jarre scienziato, si deve considerare quel¬ 
la di Gianni Jarre docente. Era docente efficace e brillante; ambidestro, 
sapeva disegnare sulla lavagna contemporaneamente con entrambe le ma¬ 
ni, così che gli allievi del terzo anno potevano seguire ammirati e stupiti 
l’ingranamento di due ruote dentate, rappresentate sulla lavagna con¬ 
temporaneamente, e seguire lo spostamento del punto di contatto dei 
due profili coniugati dai denti in presa durante l’ingranamento. Sapeva 
usare un linguaggio colorito, atto a descrivere con maggiore evidenza 
il fenomeno in istudio: così ad es., a commento della proprietà da Lui 
trovata e di cui ho fatto prima cenno, che nel flusso laminare su parete 
con entalpia di questa variabile, si ha un sistematico doppio fenomeno 
di ispessimento dello strato limite, e di anticipazione proporzionale ad 
una lunghezza che memorizza il passato, Egli dice: «lo strato limite di¬ 
venta previdente coll’età». Spesso sapeva usare anche l’arguzia, conna¬ 
turale in Lui, dotato di forte spirito di humour; così ad es. nelle Sue 
lezioni di «matematica di base» per i docenti di formazione professio¬ 
nale, parlando del concetto di «infinito» così si esprime: «in tutti i casi 
prima considerati, il vero, pieno, autentico infinito è stato approssima¬ 
to ma mai raggiunto; è rimasto un lontano traguardo intoccabile, sol¬ 
tanto immaginato come passaggio al limite, quasi per un naturale bisogno 
di tranquillante completezza. Eppure c’è un caso in cui il vero infinito, 
proprio tutto, infinitamente completo è a portata di mano, concreto e 
tangibile: un docente universitario riceve uno stipendio fisso dallo Stato 
per insegnare e ricercare; pratica invece una impegnativa professione, 
che gli impone il full-time extra muros. La sua paga oraria, con zero 
ore prestate intra-muros, fa pieno centro nel pieno infinito. Egli cioè 
riesce a totalizzare paga oraria infinita, seppur statale». 

Gianni Jarre scienziato, docente, uomo: era Uomo di grande genero¬ 
sità d’animo, e questo lo posso ben testimoniare io, che ne ho avuto tante 
prove in tanti anni di comune militanza in questo Politecnico; Uomo 
di pensiero, riteneva che ogni attività umana doveva tendere al progres¬ 
so spirituale e materiale dell’uomo. Per questo apprezzava grandemen¬ 
te il motto Veritas ed Utilitas della nostra Accademia, della quale era 
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fiero di farne parte, e per la quale non aveva risparmiato tempo e fatica 
per dare il Suo contributo alla celebrazione del bicentenario di essa; per 
questo, nella Sua magistrale Relazione «Il contributo deirAccademia allo 
sviluppo delle Scienze Tecniche» nel Convegno tenuto a ricordo dei pri¬ 
mi due secoli deH’Accademia stessa, e così pure nella altrettanto magi¬ 
strale Relazione scritta in collaborazione con Vittorio Marchis 
«Accademici o Tecnologi?» pubblicata nel volume «Tra Società e Scien¬ 
za», ed ancora nella bella Conferenza, da Lui scritta e proprio qui, in 
questa sala, letta purtroppo a Sua morte già avvenuta, aveva messo be¬ 
ne in evidenza quanto l’Accademia aveva fatto per onorare il suo motto. 

Uomo molto sensibile ai problemi umani e sociali, oltre che a quel¬ 
li scientifici, riteneva che nel passaggio dalla Veritas alla Utilitas il con¬ 
tributo deìì’Homo Faber è da considerarsi di pari dignità di quello 
del VHomo sapiens. Per questo nel suo opuscolo delle lezioni al Corso 
di aggiornamento dei docenti professionali, col pretesto di ricavarne il 
soggetto per un problema, aveva riportato per educare gli animi, oltre 
che istruire la mente, il testo della lapide posta nel cortile dell’Universi¬ 
tà di Padova: 

Dal 1938 al 1942 — muratori e manovali padovani — in 150.000 gior¬ 
nate lavorative — costruirono e rinnovarono — dalle fondamenta — 
queste auguste case di Minerva — gli studenti nell’Ateneo — ricordino 
ed onorino — con il lavoro dell’intelletto — tanta nobile sudata fatica. 

Poiché credo nell’insegnamento di S. Agostino, che un morto non 
è morto se la sua presenza è in noi, lasciate che mi rivolga per un mo¬ 
mento a Gianni per dirGli: Caro Gianni, nel dedicarmi una Tua opera 
nel Natale del 1982 Tu scrivesti: «All’amico Carlo confidando nella sua 
benevolenza». Le stesse parole ora io rivolgo a Te: hai sentito le mie 
parole, confido nella Tua benevolenza, perché se le mie parole sono sta¬ 
te inadeguate, inadeguato non è il sentimento che le ha ispirate. 

Con questo stesso sentimento rinnovo il più vivo cordoglio per la scom¬ 
parsa da questa vita terrena dell’indimenticabile amico Gianni, e rinno¬ 
vo le più commosse condoglianze alla gentile Signora Marina e ai suoi 
Figli. 




Testimonianza di Rolando Rigamonti 


Nel pomeriggio di uno dei primi giorni del luglio 1970 il Corpo Acca¬ 
demico del Politecnico era riunito per l’elezione del Rettore, a seguito 
dell’improvvisa scomparsa del prof. Capetti. Due erano i candidati: il 
prof. Giovanni Jarre, portato da un gruppo di colleghi, e chi vi parla, 
che controvoglia (ormai posso dirlo)! cedendo alle pressioni di un altro 
gruppo aveva accettato di presentarsi. 

Jarre si era preparato coscientemente a quel momento ed aveva scrit¬ 
to il suo programma in una decina di pagine che lesse al Corpo Acca¬ 
demico. 

Quando, terminato lo scrutinio, risultò che, per bontà della maggio¬ 
ranza dei colleghi, io ero stato prescelto, Jarre, con molti spirito sporti¬ 
vo, pubblicamente mi augurò buon lavoro ed aggiunse: «Rigamonti, ti 
faccio omaggio del mio programma» e mi porse il suo manoscritto, che 
mi fu prezioso. 

Per tale suo gesto, e per mostrare, a coloro che non mi avevano vota¬ 
to, che apprezzavo il loro candidato, qualche giorno dopo chiesi a Jarre 
se volesse fare il prorettore. Egli accettò senza riserve. Successe però che, 
anziché essere io ad aver bisogno del suo aiuto, fu lui che si rivolse a 
me più volte per parlarmi del Collegio Universitario. 

Infatti, il prof. Renato Einaudi era ancora presidente di tale istitu¬ 
zione; aveva però, per motivi di salute, lasciato Torino da sei anni, tra¬ 
sferendosi a Pisa, alla cattedra di meccanica razionale di quella Università. 
Non potendo perciò seguire giornalmente l’andamento del Collegio, ne 
aveva affidata la cura al prof. Jarre. 

Ma purtroppo, nel 1969-70 era scoppiata la contestazione studente¬ 
sca, che dagli Atenei aveva raggiunto il Collegio, e Jarre dovette sop¬ 
portarne l’urto. 

Fu un urto violento, con occupazioni del Collegio che obbligarono 
anche a far intervenire la polizia. Jarre se ne preoccupò moltissimo, cer¬ 
cando tutti i modi per accontentare, entro i limiti del possibile, le richie¬ 
ste degli studenti: accolse certe richieste di miglioramento delle mense; 
cercò di ospitare un maggior numero di studenti nelle sezioni del Colle¬ 
gio, acquistando sacchi a pelo e letti e scovando posti ove sistemarli, 
sia pure senza tutte le comodità di una camera singola; giunse a mettere 
a disposizione altri posti letto presso l’albergo San Silvestro, sito in piazza 
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Statuto ed ora scomparso, tramite una convenzione con l’amministra¬ 
zione del medesimo. 

Jarre veniva spesso a parlarmi di questi suoi problemi, che lo assilla¬ 
vano ed ai quali dedicava con passione tutta la sua attività, e mi chiede¬ 
va se faceva bene a compiere certi passi, a prendere certe decisioni; c’era 
in lui una continua preoccupazione per far rispettare l’ordine delle isti¬ 
tuzioni e nel contempo cercar di soddisfare certe richieste degli studen¬ 
ti, che riconosceva giuste. Andammo anche, insieme, qualche volta, a 
parlare al prof. Einaudi all’albergo Suisse (ora anch’esso scomparso, pres¬ 
so la Stazione), dove egli alloggiava quando veniva a Torino. 

Ma gli studenti non capirono (o non vollero capire!) lo sforzo che 
Jarre faceva per venire loro incontro, e questa misconoscenza destò nel 
suo animo una profonda amarezza, che si ripercosse sulla sua salute. 

Mi accorsi, che il suo spirito stava vacillando e perdendo il controllo 
di se stesso, una sera, dopo cena; eravamo da poco usciti dal Suisse, 
dopo un incontro con il prof. Einaudi e, ritornando verso casa, stava¬ 
mo ancora discutendo dei problemi del Collegio e di come risolverli, quan¬ 
do improvvisamente, preso da un impeto di riconoscenza, mi disse: 
«Rigamonti, ti ringrazio del tuo aiuto; lascia che ti abbracci» e compì 
tale gesto. 

Questo episodio, e le circostanze in cui avvenne, mi colpirono pro¬ 
fondamente e mi fecero in un lampo vedere quanto dolore avesse pro¬ 
vocato in lui, e lo avesse sconvolto, l’ingratitudine di coloro che cercava 
di aiutare, e con quanto abbandono egli si attaccasse a chi gli dava una 
mano. 

Poco dopo, il suo stato depressivo era tale da costringerlo a parecchi 
mesi di riposo. 

Ritornò più tardi al suo lavoro al Politecnico; al Collegio la contesta¬ 
zione era finita e nel 1972 era subentrato un nuovo presidente, l’ing. Gi¬ 
no Castelnuovo. Sopraggiunse però alla fine del 1973, una nuova 
disposizione legislativa, in virtù della quale il prorettore entrava a far 
parte di diritto del Consiglio di Amministrazione dell’Ateneo. 

Jarre partecipò, così, a qualche seduta del Consiglio; ma un giorno 
venne a trovarmi ed a dirmi che non si sentiva di affrontare le questioni 
amministrative sulle quali il Consiglio doveva discutere e deliberare, e 
mi pregava pertanto di accettare le sue dimissioni dalla carica di proret¬ 
tore che gli avevo affidato. 

Sembrava, dalle sue parole, che avesse perfino timore di affrontare 
tali questioni. Forse era il ricordo, ancora cocente, di quanto aveva sof¬ 
ferto al Collegio; o forse era il desiderio, od il proposito, di tornare ad 
interessarsi pienamente della ricerca e degli studenti, di quegli studenti 



TESTIMONIANZA 


per il bene dei quali aveva agito al Collegio e, prima ancora, in occasio¬ 
ne della riforma del 1960 degli studi di ingegneria, quando aveva dedi¬ 
cato molto del proprio tempo ad una revisione meticolosa dei curricula 
degli allievi, affinchè molti potessero laurearsi già nell’anno accademi¬ 
co 1960-61 con i nuovi titoli accademici, come vi dirà più diffusamente 
tra breve il prof. Buzano. 

Non riuscii a dissuaderlo dalla sua idea di dimettersi da prorettore: 
così lo persi come collaboratore; non però come amico, ed è per tale 
motivo che ho accettato di ricordarlo con questo mio breve intervento. 




Testimonianza di Pietro Buzano 


Gli interventi precedenti hanno già messo ampiamente in luce le Sue 
doti di scienziato, di docente, di uomo: perciò mi limiterò al ricordo di 
tre particolari momenti di miei rapporti con Giovanni Jarre. 

1948. Giunto da poco al Politecnico (proveniente nel ’45 dalla Facol¬ 
tà di Scienze) mi trovo costretto con tutti i Colleghi nell’unica sede del 
Valentino dove quello che dovrei chiamare «il mio Istituto» consta di 
una sola stanza con qualche scaffale i cui vetri sono sostituiti da fogli 
di compensato. Un tramezzo sottile separa quella stanza da un’aula at¬ 
tigua dove tiene le sue lezioni Giovanni Jarre: la sua voce sonante supe¬ 
ra quella barriera e giunge fino a me, che mi trovo ad ascoltarlo — prima 
involontariamente — ma poi sempre più attratto dal suo modo brillan¬ 
te, ma profondo, di affrontare difficili argomenti di meccanica. Così 
in certo senso posso considerarmi anch’io «allievo» di Giovanni Jarre. 

1958. Nell’inaugurare la nuova sede di corso Duca degli Abruzzi, il 
rettore Antonio Capetti annuncia la nomina di Jarre a professore straor¬ 
dinario di Aerodinamica: come tale Jarre è anche segretario del Consi¬ 
glio di Facoltà per il triennio 58-61 che segna una svolta importante 
nell’ordinamento degli studi di ingegneria. Infatti sta per arrivare in porto 
il progetto di riforma elaborato da una commissione presieduta dal prof. 
A. Capocaccia. 

Il 31 gennaio ’60 esce il D.P.R. col nuovo ordinamento che reca im¬ 
portanti novità: corso di studi quinquennale in tutte le sedi, istituzione 
di nuove lauree (in luogo delle tre sole preesistenti: civile, industriale, 
mineraria), possibilità per le Facoltà di attivare indirizzi, alleggerimen¬ 
to degli sbarramenti. Jarre si mette all’opera per adeguare al D.P.R. il 
nostro Statuto: già il 22 marzo ’60 una prima bozza è all’esame della 
Facoltà, che affida a Jarre la stesura definitiva, approvata il 15 maggio 
’60. Il novembre successivo vengono definite le norme transitorie. A Jarre 
il compito di ricostruire le carriere scolastiche degli studenti optanti per 
il nuovo ordinamento: gli siamo a fianco io e qualche altro Collega, ma 
Jarre assume in prima persona l’onere maggiore e lo assolve con tanto 
impegno che, nella sede torinese, l’opzione risulta totalitaria e con pie¬ 
na soddisfazione degli studenti. 
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1968. Sono stato chiamato da poco alla Presidenza della Facoltà di 
Ingegneria quando nel marzo ’68 esplode anche da noi la contestazione 
studentesca che già covava in altri Atenei. Il 4 marzo ’68 gli studenti 
avanzano richieste che alcuni Colleghi giudicano pesanti: formazione di 
commissione consultive, istituzione dei semestri, modifiche di Statuto 
che diano luogo a maggiori possibilità di scelte individuali. Il Ministero 
si mette in disparte con una circolare che dice e non dice, concede e non 
concede. Io propongo alla Facoltà di accogliere nei limiti del possibile 
le richieste degli studenti, essenzialmente per dar prova che non inten¬ 
diamo trincerarci dietro formalismi procedurali. Così il 28 marzo il Con¬ 
siglio decide di procedere sollecitamente a modifiche di Statuto: ma è 
una decisione sofferta per vari docenti che ne sono poco convinti. In 
tale delicata situazione Jarre fa parte di un gruppo di Colleghi che si 
schierano al mio fianco sorreggendomi con la loro approvazione e coi 
loro consigli. Anche per questo appoggio io gli confermo oggi la mia 
riconoscenza. 



Testimonianza di Sandra Oggiano 


Desidero ricordare il prof. Jarre, che per 14 anni è stato Direttore 
del Centro Studi sulla Dinamica dei Fluidi del CNR e la cui improvvisa 
scomparsa è stata per tutti noi, ricercatori e tecnici del Centro, una per¬ 
dita tanto più dolorosa in quanto inaspettata. 

Cercherò di improntare il mio intervento alla brevità, e questo non 
già perché non abbia di Jarre molti ricordi, ma per rendere omaggio al 
carattere essenziale della Sua personalità, e cioè l’assenza di retorica. 

Chi ha conosciuto il prof. Jarre ricorda certamente come egli fosse 
la negazione della prolissità: un discorso di Jarre era generalmente bre¬ 
ve perché egli rispondeva alle domande rivoltegli prima ancora che que¬ 
ste venissero completamente formulate. 

Di una grande lucidità di pensiero egli prevedeva talora anche la pos¬ 
sibile, ulteriore risposta del suo interlocutore, cosicché quest’ultimo per¬ 
cepiva immediatamente il fascino di un modo così intelligente di 
comunicare. 

Questa caratteristica, quando il discorso si trasferiva sul piano scien¬ 
tifico, poteva rappresentare anche una notevole difficoltà per il suo in¬ 
terlocutore; si trattava di allenare il cervello ad un grado di elasticità 
sufficiente per poter cogliere in questo suo rapido dialogare l’enorme 
densità di informazioni che vi erano contenute. 

Entrando nell’ufficio del prof. Jarre sono sempre stata colpita dalla 
grande varietà di testi che occupavano, disordinatamente, la sua scriva¬ 
nia e che spaziavano nei vari campi dell’ingegneria manifestando come 
un bisogno di rinnovarsi intellettualmente studiando a fondo certi temi 
sui più svariati argomenti, disinteressandosi completamente (o quasi) del¬ 
l’aspetto utilitaristico costituito da prospettive di pubblicazioni o di pre¬ 
sentazioni a congressi ed esprimendo piuttosto il desiderio di allargare 
le proprie conoscenze e soddisfare le proprie curiosità, quasi una impli¬ 
cita esortazione a perseguire lo studio come fine a se stesso. 

E proprio questa molteplicità di interessi permetteva al prof. Jarre 
di inquadrare i più vari problemi con immediatezza, esprimendo al ri¬ 
guardo giudizi rapidi ed essenziali, ma sempre precisi (ed a volte pun¬ 
genti) sia sulle persone che su una certa maniera piuttosto superficiale 
di condurre la ricerca. 

Soprattutto negli ultimi anni mi è parso che il prof. Jarre fosse costretto 
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dai molti impegni a cui non si era potuto sottrarre, a lottare per conqui¬ 
starsi lo spazio necessario per approfondire i propri campi di interesse, 
lontano comunque e sempre da una cultura prefabbricata o improntata 
alla moda del momento. 

Ora che non sentiamo più risuonare il suo passo affrettato per i cor¬ 
ridoi del Politecnico, quella sua eterna fretta pare colorarsi nel nostro 
ricordo di un triste significato, quasi presagisse per le molte cose che 
egli aveva ancora da dire e da fare il pochissimo tempo a sua disposizione. 

A noi tutti del Centro Studi rimane il rimpianto per quello che il no¬ 
stro Direttore avrebbe ancora potuto insegnarci. 



Testimonianza di Vittorio Marchis 


Incontro Gianni Jarre — e con questo intendo dire ho incominciato 
a conoscerlo — intorno ai primi anni 80 quando due fatti concorrono, 
in maniera alquanto fortuita: il dibattito sulla formazione del Diparti¬ 
mento che avrebbe dovuto nascere dal vecchio Istituto di Meccanica Ap¬ 
plicata, Aerodinamica e Gasdinamica, e la preparazione della mostra 
su Alberto Castigliano. 

Io, allora, sono semplicemente un apprendista nel difficile mestiere 
dell’universitario, e Gianni Jarre (che un po’ assomiglia nell’aspetto a 
Prospero Richelmy, professore di idraulica, uno degli artefici di questo 
Politecnico più di un secolo fa) subito ha molte cose da raccontarmi ed 
insegnarmi. 

Non vorrei essere travisato, ma ci troviamo quasi l’uno come l’ap¬ 
prendista, l’altro lo stregone. L’ingegneria sì: i lunghi discorsi, di vivaci 
dialoghi, sulla termodinamica dei gas, sul «tirare a secco» l’aria degli 
impianti frigoriferi, sempre un caloroso, epidermico aggancio con la real¬ 
tà, con la fisica dei tecnici e degli artigiani, con la «geometria dei pia¬ 
strellisti», con la simmetria dei cristalli e delle varie maniere di sezionare 
un cubo; la complicità con cui un sabato riusciamo a convincere un an¬ 
ziano soffiatore di vetro dell’Istituto di Chimica a costruirci — vera e 
di vetro — una bottiglia di Klein, sino ad allora vista solamente sui libri. 

Ho parlato dell’ingegneria, ma non solamente... L’arte, che forse in¬ 
vidiava in altri che gli erano vicini, la letteratura, da cui era più schivo 
e di cui ebbi solo fugaci lampi, la storia in cui ritrovava l’entusiasmo 
del detective. 

Il 6 ottobre dell’anno scorso, di sera, ancora discutevamo delle vi¬ 
cende, peraltro oscure e ambigue, di un fuoriuscito piemontese che Bab- 
bage si portò a Torino, a quando fu invitato a Corte, nella prima metà 
dell’800. La storia di Antonio Prandi è ancora tutta da scrivere... 

Di queste brevi note, l’Accademia delle Scienze è una delle quinte, 
forse anche gran parte del fondale. Nell’Accademia, che ho imparato 
a conoscere come timido frequentatore introdotto proprio dal profes¬ 
sor Jarre, ho scoperto, anche tra gli scaffali polverosi e sui pavimenti 
che i recenti restauri avevano denudato lasciando libere le volte, le mol¬ 
te facce del sapere e della scienza e i tortuosi sentieri del suo rinnovamento. 

Con la Mostra del Bicentenario dell’Accademia lavoriamo di nuovo 



VITTORIO MARCHIS 


assieme a lungo, ma sempre separatamente, per poi discuter animosa¬ 
mente alla sera delle congetture e dei parziali risultati raggiunti. Molti 
dei problemi, e non solo storici, che quasi per sfida Gianni Jarre mi pro¬ 
pose, devo ancora affrontarli. 

Gianni più volte mi esternò la sua soddisfazione nel rileggere i necro¬ 
logi — i più belli — degli accademici illustri, che ricevevano l’onore del 
ricordo. Questa giornata è stata pensata anche in questo senso. 

Il ricordo e la fantasia spesso camminano insieme. Per chi è vissuto 
accanto a Gianni Jarre, me lo auguro di cuore, avranno un significato 
tutto particolare nomi come Giannone, Zumstein, Quintino Sella, Ger¬ 
mano Sommeiller, Diodata Saluzzo, il padre Beccaria, Plana, Carlo Igna¬ 
zio Giulio, ma anche Vito Volterra, Sclopis, Bepi Colombo e molti altri 
che ora non ricordo. 


Oscillazioni, instabilità e caos 


Silvio NOCILLA* 


Riassunto. Vengono richiamati brevemente i concetti basilari relativi alle 
oscillazioni non lineari, alla stabilità e all’instabilità. Vengono poi presen¬ 
tati i moti caotici, segnalando difficoltà concettuali e problematiche ine¬ 
renti alla loro comprensione, e mostrando con esempi la loro struttura e 
la loro formazione. 

Summary. The basic concepts regarding thè nonlinear oscìllations, thè sta- 
bility and unstability are briefly recalled. The chaotic motions are then in- 
troduced, indicating conceptual difficulties and problems for understanding 
them, and showing with examples their structure and formation. 


1. Vibrazioni e stabilità secondo Liapounov 

Si legge nei testi di storia della scienza che Galileo, osservando nel 
duomo di Pisa le oscillazioni dei lampadari appesi al soffitto, abbia sco¬ 
perto la legge dell 'isocronismo secondo cui il periodo delle oscillazioni 
non dipende dalla loro ampiezza. Nella prima giornata del «Dialogo delle 
nuove scenzie» Galileo fa pronunciare a Salviati queste testuali parole: 
«Però notate: slargato il pendolo del piombo, v.g., cinquanta gradi dal 
perpendicolo e di lì lasciato in libertà, scorre, e passando oltre al per¬ 
pendicolo quasi altri cinquanta, descrive l’arco di quasi cento gradi, e 
ritornando per se stesso indietro, descrive un altro poco minore arco, 
e continuando le sue vibrazioni, dopo gran numero di quelle si riduce 
finalmente alla quiete. Ciascheduna di tali vibrazioni si fa sotto tempi 
eguali, tanto quella di cinquanta, di venti, di dieci e di quattro; ...» 

Oggi scriviamo in forma quantitativa tale legge nel modo seguente: 


(1) T=T 0 = 2'K^J7g 


essendo / la lunghezza del pendolo e g l’accelerazione di gravità. 


* Dip. di Ing. Aeronautica e Spaziale, Politecnico di Torino, C.so Duca degli Abruzzi, 24, 10129 
Torino, Italia. 
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Oggi sappiamo anche che strettamente parlando tale legge è vera so¬ 
lo al limite quando l’ampiezza delle oscillazioni tende a zero poiché 
la legge rigorosa (in assenza di forze passive) è rappresentata nel grafico 
di Fig. 1, che presenta tangente orizzontale nelForigine, ed un asintoto 
verticale per = ir. Questa osservazione fondamentale di Galileo può 
essere considerata l’atto di nascita della meccanica delle vibrazioni. Sap¬ 
piamo oggi che i sistemi dinamici lineari posseggono tanti periodi pro¬ 
pri di vibrazione T r (r= 1,2, ..., n ) quanto è il numero n di gradi di 
libertà del sistema, ed il calcolo di tali periodi (o dei loro inversi, le fre¬ 
quenze proprie) è il primo compito dell’analisi modale, capitolo della 
meccanica di grandissima attualità e importanza. 




Fig. 1 


Nei sistemi non lineari l’analisi modale diventa molto più complicata per¬ 
ché, come già oservato in Fig. 1 per n= 1, non abbiamo più n frequenze ben 
distinte, bensì n frequenze dipendenti in modo essenziale dalle ampiezze di vi¬ 
brazione delle diverse coordinate lagrangiane q r (r= 1,2, ..., n ) che caratte¬ 
rizzano istante per istante la configurazione del sistema. Inoltre queste 
vibrazioni non sono più sinusoidali come nel caso lineare, ma, nella ipotesi 
che siano periodiche, posseggono una loro specifica forma d’onda, che 
dal punto di vista analitico può essere espressa mediante serie di Fou- 
rier, ossia somma di infinite armoniche con frequenze crescenti in pro¬ 
gressione aritmetica. Sono molto espressive le seguenti rappresentazioni 
di un moto periodico x(t), v. Fig. 2: 
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a) sul piano delle fasi (x, x) 

come linee caratteristiche chiuse in se stesse; 


b) come spettro di Fourier 

discreto con frequenze in progressione aritmetica; 

c) come mappa di Poincaré 

ossia dall’insieme dei punti sul piano delle fasi corrispondenti ad istanti 
t h = t Q + h T(h= 1,2, ...) intervallati di un periodo T. Nel moto pe¬ 
riodico la mappa di Poincaré è costituita da un solo punto. 



Se il moto è quasi-periodico lo spettro di Fourier è ancora discreto, 
ma le frequenze non sono più in progressione aritmetica. 

Se il moto vibratorio non è periodico, le caratteristiche sul piano del¬ 
le fasi non sono più linee chiuse, lo spettro di Fourier è continuo, e la 
mappa di Poincaré è costituita da un insieme più o meno disordinato 
di punti. Soprattutto nei moti caotici la mappa di Poincaré costituisce 
uno strumento assai valido di osservazione e di studio. 

L’ambiente naturale per una visione completa dei moti di un sistema 
con n gradi di libertà è lo spazio delle fasi con 2 n dimensioni e coordi¬ 
nate (< 7 j, ..., q„; q x , ..., q„). Proprio in tale spazio delle fasi si introduce 
il concetto fondamentale di stabilità secondo Liapunov, ossia rispetto 
ai dati iniziali, mediante il procedimento seguente, vedi ad es. [2], Sia: 

(2) x=F(x, t), xtR m 

l’equazione del moto del sistema, e: 

x=x 0 


( 3 ) 
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il valore assunto da x all’istante iniziale, e sia x e un punto di equili¬ 
brio della (2), ossia un vettore costante tale che: 

(4) F(x e , t) = 0, vf£t 0 . 

Si pone la seguente definizione di stabilità. Il punto di equilibrio 
x ( , è stabile qualora dato ad arbitrio un numero positivo e piccolo a pia¬ 
cere è possibile trovare un 8 tale che 

(5) se II Xp — x e II < 8 allora II x(t)-x e \\ < e, V t > t 0 . 

La norma II . Il che interviene nella definizione è solitamente la nor¬ 
ma euclidea: 


(6) 11x11= fjx r 2 

ma può essere scelta in modo diverso a seconda dei casi. 

Il profondo significato fisico di questa definizione di stabilità sta nel 
fatto che in realtà nei problemi dinamici i dati iniziali sono sempre af¬ 
fetti da qualche incertezza, e possono fluttuare entro certi limiti, peral¬ 
tro molto ristretti di regola. Per un effettivo controllo del fenomeno è 
essenziale essere certi a priori che pur in questa situazione di incertezza 
dei dati iniziali le traiettorie dinamiche nello spazio delle fasi si manten¬ 
gono sempre prossime al punto di equilibrio x e , al crescere del tempo 
fino all’infinito. Più precisamente si vuole essere certi che tali traietto¬ 
rie dinamiche si mantengono sempre interne ad una ipersfera di centro 
x e e raggio e arbitrariamente piccolo. Quale prezzo si deve pagare per 
avere questa certezza? 

Nei punti di equilibrio stabile secondo Liapunov il prezzo è invero 
molto modesto, perché basta che i dati iniziali x 0 stiano anch essi den¬ 
tro un’ipersfera di centro x e e raggio 8, di regola dipendente d e. 

È importante osservare che in questa definizione: 

1) nessun’altra condizione è richiesta ai dati iniziali x a , i quali pertan¬ 
to sono delle costanti arbitrarie, purché interni all’ipersfera di rag- 
gio 8. 

2) Nessun’altra condizione è richiesta alle traiettorie dinamiche fuor¬ 
ché di restare contenute nell’ipersfera di raggio e al passare del tem¬ 
po. Dentro questa specie di cappello di Liapunov è consentito alle 
traiettorie dinamiche di fare tutto quello che vogliono, di avere anche 
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i comportamenti più bizzarri, senza che con questo venga intaccato, per 
definizione, il carattere di stabilità del punto di equilibrio. 

Esistono però dei casi, anche molto semplici, in cui pur pagando il 
prezzo suddetto non è possibile raggiungere lo scopo di costringere le 
traiettorie dinamiche a stare entro una ipersfera di raggio arbitrariamente 
piccolo. Entriamo allora nel campo minato della instabilità. Un esem¬ 
pio tipico è dato ancora dal pendolo, nella configurazione di equilibrio 
d = ir, punto più alto A. Basta infatti uno spostamento, o una velocità 
iniziale, piccola fin che si vuole a partire da A perché il punto P si allon¬ 
tani di quantità (e con velocità) non arbitrariamente piccole. 

Se il sistema è conservativo, con potenziale U(qi , ..., q„), è ben no¬ 
to il teorema di Dirichlet che fornisce una conduzione sufficiente per 
la stabilità nel senso sopra precisato. Precisamente posto: 

(7) x r = q r -q r ; x r+n = q r (r- 1, 2, ..., n) 

dove le q r sono i valori delle q r in condizioni di equilibrio statico, allo¬ 
ra il punto di equilibrio x e = 0è stabile secondo Liapunov se in corri¬ 
spondenza dei valori q r il potenziale U presenta un massimo effettivo. 
Se in luogo del potenziale U si introduce Venergia potenziale: 

(8) W=-U 

si avrà stabilità nei punti di minimo. 

Nel caso del pendolo l’energia potenziale W(d) vale: 

(9) W(à) = - mgl cos d 

ed ha l’andamento indicato in Fig. 3a con minimo in B (# = 0, equili¬ 
brio stabile) e massimo in A (d = ir, equilibrio instabile). 

È interessante e utile per lo studio dei moti caotici un’altra forma di 
energia potenziale W(x), con un massimo e due minimi adiacenti, che 
rappresenta una doppia buca di potenziale, v. fig. 3b. 

Accenniamo infine ad un ulteriore, fondamentale aspetto della stabi¬ 
lità o instabilità. Supponiamo che un sistema dinamico del tipo (2) am¬ 
metta una certa soluzioine x=x(t). Supponiamo poi che tale situazione 
venga perturbata modificando i dati iniziali, e sia: 

(10) x=x(t) + òx(t) 


la soluzione perturbata. Quale comportamento presenta la perturbazione 
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8x(t) al crescere del tempo t? Anche qui si tratta di introdurre un’op¬ 



portuna norma per 8 x(t), e di investigare in quali condizioni essa si man¬ 
tiene arbitrariamente piccola, o comunque limitata, o addirittura tende 
a zero per t -» + oo. Se la x(t) è periodica ci si riconduce allo studio di 
equazioni a coefficienti periodici, come ad esempio l’equazione di 
Mathieu: 


(11) x+(X + /i cosai 0^ = 0 

le cui soluzioni possono avere comportamenti del tutto diversi per 
t-> + oo, a seconda dei valori assunti da Xe g. 


2. L’instabilità 

Le due tipiche configurazioni di equilibrio del pendolo, o della dop¬ 
pia buca di potenziale, evidenziate nel numero precedente in A e in B, 
possono anche essere riguardate nel modo seguente. 

Ammettiamo, come ipotesi di lavoro, che le traiettorie dinamiche sul 
piano delle fasi {x, x) possano essere contenute in cerchi, con centri ri¬ 
spettivamente ( x A , 0) e (x B , 0), di raggio arbitrariamente piccoli, e quin¬ 
di che sia lecita la linearizzazione dell’equazione del moto nella forma: 

(12A) y-a 2 y = 0 in prossimità di A, con y = x-x A 

(12B) y + b 2 y = 0 in prossimità di B, con y=x-x B 
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Queste equazioni ammettono le soluzioni semplici 
(13) y = e ±at in zi; y = e ±,bt in B 

e pertanto nel caso stabile B le traiettorie dinamiche, in forma di ellisse, 
sono certamente racchiudibili nel cerchio suddetto, purché i dati iniziali 
siano contenuti in un conveniente intorno del suo centro. Dunque Vipo¬ 
tesi di lavoro era corretta e la linearizzazione lecita, nel senso che forni¬ 
sce un modello corretto del fenomeno. Viceversa nel caso instabile A, 
essendo l’esponenziale reale e al divergente per t -*■ + oo, non è possibile 
trovare un intorno di (x A , 0) tale che qualora i dati iniziali siano conte¬ 
nuti in esso le corrispondenti traiettorie dinamiche stiano tutte rinchiu¬ 
se in un cerchio di raggio finito. Dunque in tal caso l’ipotesi di lavoro 
non era corretta, e la linearizzazione non è lecita, perché fornirebbe un 
modello del fenomeno in contraddizione con l’ipotesi fatta per ricavarla. 

Le cose poi si complicano nei casi con più gradi di libertà in cui il 
potenziale può, in un suo punto di stazionarietà, non essere né massimo 
né minimo (v. ad esempio i punti di sella). 

Ci troviamo di fronte ad una prima domanda, solo apparentemente 
banale. Un punto di equilibrio non stabile è instabile ? Una risposta ade¬ 
guata richiederebbe ulteriori approfondimenti, che qui non possono tro¬ 
vare spazio. Bisognerebbe dare una definizione in positivo del concetto 
di instabilità, diversa dalla negazione della stabilità. Pur tuttavia osser¬ 
viamo che l’elemento chiave della discussione precedente è che nel caso 
della stabilità ci troviamo di fronte a soluzioni esponenziali con espo¬ 
nenti caratteristici immaginari puri ± ib; nel caso della instabilità con 
un esponente positivo + a. Basta la presenza anche solo di un esponente 
caratteristico positivo, o meglio con parte reale positiva, perché tutto 
il castello della stabilità crolli, nel senso che perturbazioni piccole fin 
che si vuole dalla configurazione di equilibrio x e portano nello spazio 
delle fasi a punti infinitamente lontani da x e per t -> + ° °. L’esponente 
caratteristico con parte reale positiva ha dunque il ruolo di amplificato- 
re infinito di perturbazione. Di questa situazione si accorge anche il «com¬ 
puter» cui venga affidato il compito di integrare numericamente il sistema 
dinamico in esame, segnalando « overflow». 

Per porsi al riparo da questi guai si assume solitamente come criterio 
sicuro per la stabilità quello che tutti gli esponenti caratteristici abbiano 
parte reale negativa. In questo contesto giuoca un ruolo basilare il crite¬ 
rio di Routh e Hurwutz. 

In tutto quanto detto finora sulla instabilità giova ancora osservare 
che è necessario avere a che fare con soluzioni di tipo esponenziale, pre- 
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senti soltanto nei sistemi lineari. Nel caso non lineare, se è possibile ricon¬ 
dursi in forma approssimata alla stessa situazione, si parla in instabilità 
locale. Possono però presentarsi situazioni del tutto nuove, come quella 
del ciclo limite , di cui un esempio significativo è dato dall’equazione di 
Van der Poi: 

(14) x-fi{\-x 2 )x + x=0 0* > 0) 

L’origine nel piano delle fasi (x, x) è localmente instabile; però per 
t -> + oo le caratteristiche non vanno all’infinito, ma tendono ad una spe¬ 
ciale linea chiusa, anch’essa caratteristica, detta appunto ciclo limite, 
v. Fig. 4. 



Fig. 4 
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3. Il caos questo sconosciuto 

Io non so se abbia molto senso accingersi a spiegare cosa sono i moti 
caotici, cioè cercare di esporre con la razionalità e la chiarezza necessa¬ 
rie per trattare un argomento scientifico, dei fatti e dei concetti che si 
presentano come negazione di razionalità e chiarezza. Sfogliando la let¬ 
teratura sul «caos» copiosissima in questi ultimi tempi, si avverte l’esi¬ 
genza di segnalare una tematioca sconcertante, se pur tanto viva e 
interessante, mettendo però in guardia il lettore dall’illudersi di avere 
capito cosa sia veramente il caos dopo la lettura dei lavori ad esso dedi¬ 
cati. Così ad esempio nell’articolo di apertura del volume «Directions 
in Chaos» [3, Voi. 1°], di J. Ford, dal titolo «Directions in classical 
chaos», leggiamo: «in questo articolo noi cercheremo di indir izz arci al¬ 
le conseguenze fondamentali che si ottengono rispondendo alle questio¬ 
ni globali: ‘Cos’è il caos... Cosa significa... Chi ne ha bisogno... A quale 
scopo serve?... ’ Specificamente noi pensiamo che rispondendo definiti¬ 
vamente a queste “ unsolved questions of thè past” noi possiamo facili¬ 
tare il nostro viaggio nel futuro. 

Tuttavia, risposte non ambigue a tali quesiti possono essere accura- 
mente formulate solo se affrontiamo direttamente Vantico paradosso re¬ 
lativo al problema di definire il caos». 

Una seconda citazione, tratta dalla prefazione del volume «Chaos and 
Complexity» [4] è la seguente: «Sebbene questi nuovi campi di ricerca 
non possano essere collocati nell’area ben definita della meccanica sta¬ 
tistica, è pur evidente che non esiste ancora una chiara definizione di 
essi in termini di comuni metodologie e concetti largamente applicabili». 

Tuttavia nelle stesse fonti ora citate si trovano anche degli spunti atti 
a gettare qualche luce sulla comprensione dei fenomeni caotici. Infatti 
dice ancora J. Ford: «Per secoli la dinamica classica è stata riguardata 
come prevedibile deterministicamente; eppure l’esistenza del caos sem¬ 
bra richiedere che la dinamica classica produca anche dei comportamenti 
imprevedibilmente vagabondi ( unpredictably erratic). La grossolana e 
facile risoluzione di questo dilemma afferma che la imprevedibilità os¬ 
servata del caos deriva dalla dipendenza esponenzialmente sensitiva dello 
stato finale dallo stato iniziale (esponenti di Liapunov positivi, metri¬ 
che di entropia positive, o simili). Questo concetto di dipendenza esa¬ 
sperata dai dati iniziali viene peraltro discusso criticamente nel prosieguo 
dello stesso articolo. Viene però ripreso con maggiore fermezza da un 
articolo di S. Isola, contenuto nel già citato volume [4], dal titolo signi¬ 
ficativo «Understanding complex behaviour. Some remarks on method 
and interpretation». Si legge in tale articolo: «Le nonlinearità presenti nel 
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problema possono infatti produrre una estrema sensibilità rispetto allo 
stato iniziale del sistema, cosicché anche un piccolo errore nella sua co¬ 
noscenza sarà amplificato esponenzialmente al crescere del tempo. Un 
sistema con questa proprietà (cioè sensibile dipendenza dalla condizio¬ 
ne iniziale) è caotico : nonostante la sua natura deterministica esso subi¬ 
sce una sorta di perdita di memoria sulla sua storia passata, cosicché 
diventano impossibili previsioni precise sul futuro. Se un’evoluzione tem¬ 
porale caotica si avvicina asintoticamente ad un attrattore, allora que¬ 
sto attrattore diventa strano ». 

Dunque il caos viene presentato come esasperata sensibilità ai dati 
iniziali, fonte di non-prevedibilità nei sistemi dinamici non lineari, sep¬ 
pur deterministici. Questo concetto appare molto importante per chi stu¬ 
dia il caos, tanto che la sua paternità viene attribuita a Poincaré, il padre 
della moderna impostazione di studio della meccanica. Si legge infatti, 
sempre nella prefazione del volume [4]: «Agli inizi di questo secolo H. 
Poincaré ha già osservato che una dinamica completamente determini¬ 
stica non implica necessariamente esplicite predizioni sull’evoluzione di 
un sistema dinamico. Ciò può essere considerato una pietra miliare nel¬ 
l’approccio allo studio del caso dinamico ». 

Queste citazioni mi suggeriscono un commento personale, e precisa- 
mente la seguente domanda: La suddetta esasperata sensibilità ai dati 
iniziali implica che la condizione di stabilità di Liapunov non sia soddi¬ 
sfatta? La risposta affermativa, e cioè assenza di stabilità secondo Lia¬ 
punov, si presenta come spontanea, d’accordo con l’impostazione logica 
del titolo della presente relazione, che propone il caos come figlio della 
instabilità. Una riflessione più attenta però può fare propendere anche 
per una risposta diversa. Infatti si può pensare a moti caotici, nel senso 
di esasperatamente dipendenti dai dati inziali, e pur tuttavia tali che nello 
spazio delle fasi siano costretti entro al «cappello di Liapunov». Il «cap¬ 
pello» non viene forato, purché i dati iniziali stiano in una conveniente 
sfera, ma dentro ad essa regna il disordine più completo, con dipenden¬ 
za esasperata dai dati iniziali, cioè appunto il caos. In ogni modo però 
si deve tener presente che nei moti caotici le dimensioni di questo «cap¬ 
pello» sono limitate, e non infinitamente grandi come per le instabilità 
del tipo descritto nel numero precedente. 

Ancora una osservazione di carattere generale. Gli studi più recenti 
hanno mostrato che per incontrare moti caotici non è necessario consi¬ 
derare sistemi complessi con molti gradi di libertà, ma bastano semplici 
sistemi con un solo grado di libertà. Bisogna però che essi presentino 
grandi non linearità, contrariamente a quanto veniva fatto per lungo tem- 
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po nello studio dei sistemi non lineari, di regola condizionati dalla pre¬ 
senza del «piccolo parametro» di non linearità. 

Segnaliamo infine che una trattazione molto ampia ed interessante 
dei moti caotici, con ricca documentazione bibliografica, è stata effet¬ 
tuata nel luglio 1986 in un Corso CISM ad Udine, dal titolo: «Chaotic 
motions in non linear dynamical Systems». Ad esso rimandiamo per svi¬ 
luppi ed approfondimenti delle tematiche che qui possiamo soltanto ac¬ 
cennare, spesso ispirandoci al Corso stesso. 

Passiamo ora a mostrare con un esempio concreto come si manifesta 
la suddetta esasperata dipendenza dai dati iniziali, fonte di comporta¬ 
mento caotico, sia sotto forma di time histories sia sotto altre interes¬ 
santi forme geometriche. L’esempio è tratto dal volume di Thompson 
[5], che riporta ricerche basilari di Ueda [6]. Si consideri l’equazione di 
Duffing {non linearizzabile ): 

(15) x+kx+x 3 = B cos t 

Le sue soluzioni x=x{t) dipendono in modo essenziale: 

— dai dati iniziali x=x 0 ; x=x 0 per t = 0; 

— dai parametri ( k, B). 

Un primo gruppo di risultati è ottenuto tenendo costanti ( k, B) e fa¬ 
cendo variare ( x a , x a ). L’equazione differenziale (15) viene integrata nu¬ 
mericamente con l’ausilio del calcolatore che consente di ottenere una 
enorme quantità e varietà di risultati in tempi limitati e costi contenuti. 
Riportiamo qualche risultato. Attribuendo a ( k, B ) i valori: 

(16) £ = 0.05; B = 1.5 

e ai dati iniziali i valori molto prossimi tra loro: 


(17) 


fx 0 = 3 f x 0 = 3.01 
[ x 0 = 4 [i o = 4.01 


si hanno i diagrammi orari di Fig. 5, e sul piano delle fasi di Fig. 6. Ap¬ 
paiono ben evidenti: 

— la complessità delle forme d’onda; 

— la dipendenza dai dati iniziali, fonte di progressiva distorsione di tali 
diagrammi. La Fig. 7 indica come si genera la corrispondente mappa 
di Poincaré beccando stroboscopicamente dopo ogni periodo 2-n- i suc¬ 
cessivi punti A, B, C,... In Fig. 8 è riportata una mappa di Poincaré 










40 


SILVIO NOCILLA 


finale, ottenuta mediante un numero grandissimo di cicli. Riguardo a 
quest’ultima osserviamo che si tratta di una struttura molto complessa, 
stratificata, conseguenza di ripetuti stiramenti e ripiegamenti (stretching 
and folding) dell’insieme delle traiettorie caotiche stazionarie. I punti 
che costituiscono la mappa si avvicendano in modo molto irregolare, 
ma non eccessivamente dispersivo, sempre con ascisse positive, costituen¬ 
do quello che viene denominato « strange attractor» o «chaotic attrac- 
tor», non privo di un certo carattere di varietà ed eleganza di figura. 
È difficile però dare una definizione rigorosa di tali attrattori. Nel volu¬ 
me già citato di Thompson [5] è detto anzitutto che è preferibile il ter¬ 
mine di attrattore caotico anziché strano, e poi si precisa: «Un attrattore 
caotico può essere definito come un attrattore limitato che stira e ripie¬ 
ga il pacchetto ( bundle ) di traiettorie finali stazionarie, producendo sen¬ 
sibile dipendenza dalle condizioni iniziali e imprevedibilità a lungo 
termine». Definizione espressiva ma certo non matematica. 

Riguardo alle mappe di Poincaré per gli attrattori caotici c’è ancora 
da rilevare quanto segue. Già abbiamo osservato nel n. 1 che se il moto 
è periodico, ed il tempo impiegato nell’operazione stroboscopica è uguale 
al suo periodo, la mappa si riduce ad un solo punto, cioè una figura 
geometrica con dimensione zero. Se detto tempo è un sottomultiplo del 
periodo, come capita quando si è in presenza di risonanza subarmoni¬ 
ca, la mappa si riduce ad una successione finita di punti che si sovrap¬ 
pongono ciclicamente. Se invece abbiamo un moto caotico i vari punti 
si susseguono in modo disordinato formando un insieme con struttura 
stratificata indefinitamente {infinitely layered), del tipo recentemente stu¬ 
diato in modo sistematico da Mandelbrot [7], il quale coniò per tali og¬ 
getti il termine di fractals' , sia a causa della loro struttura infinitamente 
‘broken’, sia perché, in base ai metodi matematici standard per la de¬ 
terminazione spaziale della dimensione di un arbitrario insieme di pun¬ 
ti, si ottengono per i frattali valori non interi. Si ha cioè qualcosa di 
intermedio tra dimensione zero (un solo punto), dimensione uno (una 
linea), dimensione due (una superficie), etc. Non possiamo qui adden¬ 
trarci in questo campo difficile della geometria, che trae radici dall’in¬ 
sieme di Cantor, che qualche decennio addietro era inteso da qualcuno 
come una «invenzione perversa» dei matematici. 

Rivolgiamo invece molto velocemente la nostra attenzione ad un se¬ 
condo gruppo di risultati ottenuti per l’equazione (15) facendo variare 
i parametri (k, B). Si tratta di distinguere, sul piano dei parametri (k, 
B), regioni diverse corrispondenti a comportamenti diversi delle sue so¬ 
luzioni e precisamente: 

— soluzioni periodiche con lo stesso periodo della forzante, dette soli- 
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tamente armoniche, anche se formate da diverse (infinite) armoni¬ 
che vere e proprie; 

— soluzioni periodiche con periodo multiplo della forzante, dette an¬ 
che subarmoniche; 

— soluzioni caotiche. 

Questo tipo di indagine, a parte l’interesse concettuale e teorico, può 
anche avere significative applicazioni pratiche qualora l’equazione dif¬ 
ferenziale in studio sia il modello matematico di realtà fisiche o tecniche 
di cui i coefficienti dell’equazione costituiscono parametri di funziona¬ 
mento o di progetto. Si riesce allora a prevedere per quali valori di tali 
parametri la situazione di caos si può presentare o no. La strada seguita 
di regola per queste indagini è lo studio delle biforcazioni, ossia il calco¬ 
lo dei valori dei parametri in corrispondenza dei quali si ha un cambia¬ 
mento radicale della struttura delle soluzioni, di pari passo con l’analisi 
delle strutture stesse. Nel caso dell’equazione [15], in base ai risultati 
di Ueda [6], il piano ( k, B ) viene suddiviso in modo assai complesso in 
molteplici regioni, talvolta con ricoperture parziali tra loro, caratteriz¬ 
zate ciascuna da una delle proprietà suddette. Non è qui posibile ripor¬ 
tare e discutere tali interessanti risultati. 

Segnaliamo invece un gruppo di ricerche sistematiche della Stupnic- 
ka [8] su equazioni ancora del tipo di Duffing, sia con che senza simme¬ 
tria dispari in (x, x). In particolare ricordiamo l’equazione non più 
simmetrica: 

(18) x+hx+x 3 =P 0 + P\ cos v t 
con: 

(19) P\ = 0.16; P,, = 0.03; ù = 0.05 

e riportiamo alcuni significativi risultati relativi alla dipendenza delle so¬ 
luzioni dalla frequenza dalla forzante v, in particolare al fondamentale 
fenomeno del ‘period doubling’, spesso segno premonitore della transi¬ 
zione al caos. La Stupnicka cerca soluzioni subarmoniche della (18) del 
tipo: 

(20) X® (0 = A 0 + A, cos v t + Al cos(j>/2 + <p) 

2 

e determina i parametri incogniti (A 0 , A 1; Al, <p) mediante bilancia¬ 
mento armonico, v. ad. es., Hagedorn [1], In Fig. 9 sono riporati i 
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valori di tali parametri in funzione di v. Si vede che per v > v 2 non si 
hanno subarmoniche 2 T, essendo T-2-ir/v, mentre per v = v 2 si presen¬ 
ta il fenomeno del raddoppio del periodo, potendosi avere subarmoni¬ 
che 2 T accanto ad armoniche T. Diminuendo ulteriormente la frequenza 
v si hanno ulteriori raddoppi del periodo T, e cioè subarmoniche con 
periodo AT ,..., 2 n T, finché dopo tutte queste biforcazioni si hanno so¬ 
luzioni caotiche, in corrispondenza di un certo intervallo di valori, inve¬ 
ro non molto ampio, della frequenza v. 

Riguardo al period doubling ci pare interessante osservare quanto se¬ 
gue. Nel caso della doppia buca di potenziale, riportata in Fig. 3b, è 
facile riconoscere che sono possibili essenzialmente due tipi di vibrazio¬ 
ni libere periodiche: 

1) con ampiezza «piccola» nelPintorno di ciascuna delle due configura¬ 
zioni di equilibrio stabile B e B' ; 

2) con ampiezza «grande» e passaggio successivo per B e B'. 

Si può dire grosso modo che il periodo di oscillazione di queste se¬ 
conde sia il doppio del periodo delle prime. Se allora si eccita il sistema 
con una forzante sinusoidale (e con raggiunta di un piccolo dissipato¬ 
re), secondo un modello matematico del tipo: 

(21) x+7Ì-yx(l -x 2 )=fcosut 

qualora si passi da oscillazioni forzate del 1 ° tipo ad oscillazioni del 2° 
tipo, sia la simulazione al calcolatore, sia l’esperienza, v. ad es. le ricer¬ 
che di Moon [9], indicano l’insorgere di moti caotici. Ecco dunque un 
altro tipo di rout to Chaos, collegato ad un raddoppio di periodo, con 
la complicità della instabilità in A. 

Diamo infine un rapido cenno ad alcune ricerche su di un sistema con 
due gradi di libertà, che generalizza l’equazione di Duffing, perché ad 
essa si riduce per b t =0: 

( x + K t x + KjX* + a { y = 0 (K u h ì )eR + -, {K ì ,a u b l )€R 
^ [y + h l y + b l x= 0 

In due memorie [10] e [11] abbiamo studiato la stabilità delle solu¬ 
zioni periodiche (approssimate). Trattandosi di un sistema autonomo si 
presenta nuovamente, come già per il pendolo, il fatto fondamentale che 
in campo non lineare non esiste alcuna frequenza privilegiata, come invece 
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sono le frequenze delle forzanti nelle equazioni (15), (18) e (21). Trat¬ 
tandosi poi di un sistema con due gradi di libertà assumono rilievo anche 
i modi di vibrazione, ossia i rapporti delle ampiezze delle soluzioni pe¬ 
riodiche x(t) e y{t), ammesso che esistano. Sono stati individuati, per 
diversi insiemi di valori dei parametri (K u K 3 , a u b u h x ), i valori delle 
possibili frequenze di vibrazioni periodiche per cui si ha stabilità o in¬ 
stabilità, riscontrando altresì che si possono presentare i due tipi fonda- 
mentali di instabilità: 

— con {x, y) tendenti ad infinito per t -* + oo 

— con (x, y) sempre limitati, ma con diagrammi orari molto irregolari. 

Per questi ultimi sono in corso ricerche al calcolatore in collabora¬ 
zione con P. Niezgodsky, per metterne in evidenza le tipiche proprietà 
di struttura. Si tenga presente che il tracciamento delle mappe di Poin¬ 
caré presenta l’incognita essenziale dell’intervallo di tempo impiegato 
nell’operazione stroboscopica per individuare i successivi punti sul pia¬ 
no (x, x) oppure (y, j), appunto perché non esiste un particolare pe¬ 
riodo privilegiato. In Fig. 10 sono riportati alcuni risultati preliminari 


ottenuti per: 


f /q = 0.5 


(23) f 

A', = 1 | 

| a, = 1.2 


^3 = 0.533 1 

1 

c b x = 0.3 

tv-l-i l 



che sembrano indicare carattere caotico del moto. 


4. Una riflessione finale 

Noi che lavoriamo in campo scientifico siamo abituati a pensare che 
2 + 2 faccia sempre 4, e ce lo sentiamo ripetere, in contrapposizione ad 
altre manifestazioni dell’ingegno umano come la letteratura, le arti, la 
medicina stessa, in cui 2 + 2 può fare qualunque cosa oltre che 4. Quando 
si parla di «caos» nel modo scientifico si ha la sensazione di trovarsi in 
questa stessa imbarazzante situazione. Dobbiamo allora pensare che la 
Natura, quella ««pàcris» che pur fornisce esempi mirabili di precisione cro¬ 
nometrica — si pensi ai moti dei corpi celesti — di fronte ad una operazio¬ 
ne elementare come 2 + 2 spesso e volentieri non sappia che pesci pigliare 
e preferisca risposte poeticamente incerte anziché inesorabilmente precise? 



OSCILLAZIONI, INSTABILITÀ E CAOS 


Vi lascio con questo interrogativo, su cui certamente l’amico Gianni 
Jarre avrebbe saputo ricamare qualcuna delle sue argute battute. 
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Nuove ricerche sulla teoria della turbolenza 


Massimo GERMANO* 


1. Introduzione 

Le teorie e le tecniche applicate negli ultimi tempi allo studio della 
turbolenza sono moltissime. Esse vanno dalla teoria delle biforcazioni 
alla teoria degli attrattori strani, dalla teoria dei frattali alla teoria della 
rinormalizzazione delle scale [1], Lo sviluppo dei calcolatori ha stimo¬ 
lato nuovi approcci che non si limitano allo sviluppo di algoritmi nume¬ 
rici ma comportano nuovi modelli rappresentativi dello stato fluido, 
diversi da quello tradizionale di Stokes-Navier, quali i cellular automa- 
ta e la dinamica molecolare [2], Per quel che riguarda i campi di appli¬ 
cazione occorre segnalare il grande interesse per la stabilità non lineare 
e la genesi della turbolenza, campi nei quali si sono raccolti importanti 
risultati, mentre per quel che riguarda la turbolenza sviluppata occorre 
dire che il quadro fenomenologico classico è sostanzialmente conferma¬ 
to [3]: essere la turbolenza un processo dissipativo nel quale tipicamen¬ 
te i processi di interazione non lineare trasferiscono energia da un moto 
di larga scala o moto medio alle scale inferiori. La dissipazione viscosa 
impedisce infine l’accumulo di energia nelle piccole scale e stabilisce l’e¬ 
quilibrio tra il trasferimento inerziale ed il drenaggio viscoso, e non è 
fuor di luogo sottolineare che questa visione d’insieme costruita con tec¬ 
niche molto elementari e con grande intuito fisico principalmente dal 
Taylor e dal Kolmogorov rappresenta sempre un sicuro punto di riferi¬ 
mento ed è sempre fonte di sicuri suggerimenti all’atto pratico. 

Sempre per quel che riguarda la turbolenza sviluppata si nota tutta¬ 
via una certa stasi relativamente in particolare ai modelli di chiusura delle 
equazioni mediate di Reynolds, sempre ampiamente utilizzati in sede pra¬ 
tica [4], ma sempre esposti a grosse critiche di fondo [5]. Considerazioni 
di scala indicano del resto che i gradi di libertà di un flusso turbolento 
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scalano in ogni direzione come R 4 dove R è il numero di Reynolds [6]; 
questo fatto può anche essere espresso nei termini di dimensione dal- 
l’attrattore strano associato ad un flusso turbolento [7] e tutto ciò prati¬ 
camente si traduce nel fatto che una simulazione diretta al calcolatore 
richiede R 4 punti di calcolo se vuole calcolare in ogni dettaglio il flusso 
turbolento in esame sino alla microscala dissipativa. Se a ciò si aggiun¬ 
ge che per avere una adeguata risoluzione temporale occorre procedere 
per almeno R 2 passi nel tempo si ha una idea della assoluta attuale im¬ 
possibilità di simulare direttamente flussi turbolenti di interesse pratico. 

Tra le equazioni di Reynolds, nelle quali tutte le scale coinvolte nel 
fenomeno turbolento sono mediate, e la simulazione diretta, nella qua¬ 
le ogni dettaglio è risolto sino alla scala dissipativa, si sono ultimamente 
inseriti i cosiddetti metodi LES, (Large Eddy Simulation), nei quali si 
cerca di calcolare esplicitamente la larga scala, o i grossi vortici, e si me¬ 
diano solamente i piccoli vortici di natura più stocastica e più universa¬ 
le. Lo sviluppo di questi metodi, che a parere dell’autore sono tra i più 
interessanti e promettenti dal punto di vista delle applicazioni, compor¬ 
ta sia un ripensamento di concetti classici quali i processi di media e di 
filtraggio delle equazioni, sia una conoscenza approfondita degli algo¬ 
ritmi di discretizzazione numerica che molto spesso si comportano di 
per sé come veri e propri impliciti operatori di filtraggio attraverso per 
esempio le loro proprietà dissipative. Appare evidente che i metodi LES 
assumono il ruolo di tecniche tese essenzialmente a modellare tutto ciò 
che la numerica non può risolvere; ad un estremo stanno le equazioni 
di Reynolds e all’altro la simulazione diretta, ed il numero di Reynolds 
e ciò che apre la forbice tra le due. Nei paragrafi che seguono si cerche¬ 
rà di illustrare più in dettaglio questi punti fondamentali. 


2. Operatori di media 

L’idea di applicare al segnale turbolento degli operatori di media al 
fine di separare in qualche modo la larga dalla piccola scala, è di vec¬ 
chia data: come ricorda Kampè de Feriet [8] le prime medie proposte 
rispettivamente dal Boussinesq e dal Reynolds furono le seguenti 


7 [ 


(1) 


f(x,t')dt' 
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(2) /(x, /;/?)= ^ J j J f(x',t)dx'dy'dz' 

dove d è «un temps assez petit» e /i 3 «a certain volume». Sempre di vec¬ 
chia data è anche la constatazione che mediare dei segnali è relativamente 
chiaro e semplice, mentre mediare delle equazioni non lineari è un’ope¬ 
razione per molti versi oscura e misteriosa. Il coraggio di saltare il pon¬ 
te e di mediare le equazioni di Stokes-Navier è dovuto al Reynolds; il 
prodotto di ciò sono le relative equazioni che nel caso di un fluido in- 
compressible sono date da 


(4) da ‘ | 1 d P | yV u d(a/ u' k ) 

dt dx k q dXj 2 ' Sp¬ 
essendo Uj = %f u- e p =p + p ' e a lungo si è cercato di capire quali so¬ 
no effettivamente gli operatori di media che le producono. Certamente 
non è un operatore spazio-temporale del tipo (1) o (2) e Reynolds stesso 
declinò le «regole del gioco» che le producono e che sono definite a li¬ 
vello di algebra operatoriale come 

V df . w _ S/ 

dt dt dx k dx k 

(6) /' = 0 ; fg=fg+fg' 

essendo/=/+/' eg = g + g'.Èdel resto evidente che queste «regole del 
gioco» sono automaticamente soddisfatte interpretando f come media 
statistica in uno spazio probabilistico oj 


(7) f(x,t)= 1 f(x,t,m) ji(w)dm 

Jo 

dove è la probabilità del verificarsi di /in un dato dominio di solu¬ 
zioni 12, ed è anche evidente che f è pure interpretabile come 
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(8) f{x)= lim ( ' f(x,t')dt' 

La strada maestra era segnata: abbandonate le vaghe proposte di medie 
spazio-temporali delle origini la trattazione della turbolenza si inqua¬ 
drava nei campi della meccanica statistica e delle funzioni aleatorie mentre 
il teorema ergodico assicurava, nel caso di turbolenza stazionaria, stret¬ 
ti rapporti con le medie temporali misurate direttamente in laboratorio, 
e non è il caso ora di entrare nel dettaglio sulle luci e sulle ombre che 
tanti anni di speculazione in questa direzione hanno prodotto nello stu¬ 
dio della turbolenza. È interessante invece seguire lo sviluppo parallelo 
di quelle prime idee, sviluppo che, all’ombra della rappresentazione sto¬ 
castica dei fenomeni turbolenti e parafrasando il titolo di un recente film 
di successo, è stato neccessariamente quello di «figlio di un dio mino¬ 
re». La documentazione in proposito è abbastanza scarsa e anche qui 
è basilare il già citato lavoro di Kampè de Fèriet [8] dal quale ricaviamo 
che nel 1924 Keller e Friedmann [9] propongono una media spazio¬ 
temporale fusione della (1) e della (2), mentre è del 1939 la proposta di 
A.A. Isakson [10] di una operazione di filtraggio nel tempo della/(f) 
decomposta in armoniche 

(9) /(?) = | A (X) cos (Kt - a)) d\ 

(10) A (X) cos (\t - a (X)) d\ 

dove E è l’insieme delle frequenze conservate. 

Ritornando a quanto detto prima, che un conto è mediare o filtrare 
dei segnali e un conto è mediare o filtrare delle equazioni, occorre dire 
che la preccupazione principale di queste ricerche, compresa quella ci¬ 
tata dal Kampè de Fèriet, è essenzialmente giustificare le equazioni di 
Reynolds su basi diverse da quelle strettamente statistiche, e non tanto 
quella di derivare e interpretare le relative nuove equazioni mediate, com¬ 
pito tanto più complesso quanto più vengono a mancare le basilari «re¬ 
gole del gioco» espresse dalle (5) e (6). In particolar e il fa tto che f 
non sia più nullo e che non sia piu vero che fg=fg+f g' complica 
notevolmente le cose. Una traccia di questo imbarazzo la si rileva nel 
libro di Monin e Yaglom sulla Meccanica della Turbolenza [11]. Nel 
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capitolo in cui è affrontato il problema delle medie si esamina innanzi¬ 
tutto la generica media spazio-temporale di tipo convolutivo 

(11) f(x,t; h, j?)= | f (x ', t ') G (x -x ', t-1' -, h, &)d 3 x' dt' 
con 

(12) [G(x-x',t-t'\h,d)d ì x'dt' = \ 


Al riguardo, e citando la traduzione inglese della MIT Press, pg. 209, 
Voi. 1, viene poi testualmente aggiunto: «The use of time, space or time- 
space averaging defined by some equation of thè form (11), is very con- 
venient from thè practical viewpoint but leads to a great many unavoi- 
dable analytical difficulties in theoretical calculations. Moreover, this 
type of averaging has thè great disadvantage that thè question of thè 
form of thè function G suitable for thè given problem must be resolved 
each time before use. For all these reasons it is desirable in thè theory 
of turbulence to avoid thè use of this type of averaging altogether, and 
to adopt instead some other method of defining thè mean value, a me- 
thod that has simpler properties and is more universal. A convenient de- 
finition of this type, which we shall use throughout this book, is found 
in thè probability-theory treatment of thè fields of fluid dynamic varia- 
bles in a turbulent flow as random fields». Queste parole sembrano ab¬ 
battersi sulle medie spazio-temporali come una pietra tombale e 
concludere così il nostro excursus storico, ma così non è. Il calcolatore 
si insinua tra le tecniche analitiche e quelle sperimentali proponendo nuovi 
punti di vista, suggerendo nuovi metodi e rivalutandone vecchi. 


3. La Large Eddy Simulation 

Il calcolatore fa vedere molte cose in luce diversa. In particolare le 
«great many unavoidable difficulties» che Monin e Yaglom prevedono 
nell’applicazione delle medie convolutive (11) «in theoretical calculations» 
non sembrano così terribili nella loro traduzione numerica grazie alle 
enormi capacità di calcolo sviluppate dai moderni calcolatori. Ci sono 
inoltre interessanti analogie tra l’applicazione di «filtri» convolutivi al¬ 
le equazioni di Stokes-Navier e le discretizzazioni numeriche: entrambe 
introducono negli algoritmi che ne derivano delle nuove scale di spazio 
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e di tempo, le lunghezze caratteristiche di griglia Ax e At e le lunghezze 
caratteristiche di filtro nello spazio e nel tempo h e d. È del resto noto 
che uno schema numerico che traduce per esempio alle differenze finite 
le equazioni di Stokes-Navier risolve in realtà una «equazione modifi¬ 
cata» [12] nella quale possono essere presenti termini dissipativi e di¬ 
spersivi scalati con Ax e At; inoltre è caratteristico del troncamento di 
griglia e degli schemi conservativi applicati a problemi non lineari il blocco 
della cascata inerziale con relativo immagazzinamento sulle scale risolte 
dell’energia non dissipata, fino all’esplosione del calcolo stesso [13]. Al 
fine di separare nettamente il processo di filtraggio dalla numerica è sem¬ 
brato opportuno [14] ricorrere alle vecchie idee di medie spazio-temporali. 
Ribattezzate come filtri le medie convolutive (11) sono state riproposte 
come mezzo per separare i valori di grande scala / da quelli di piccola 
scala /', ed è opportuno notare che ora i termini grande e piccolo assu¬ 
mono un valore ben preciso in quanto spontaneamente legati alle scale 
di griglia di calcolo Ax e At a disposizione. Indipendentemente dalla for¬ 
ma del nucleo convolutivo G notiamo che le (11) rispettano ancora le 
condizioni 


(13) 


àf_ = àj_. V _ d f 

dt dt dx k dx k 


e possiamo quindi scrivere formalmente che 


(14) 


dUj 

dXj 


= 0 


(15) 


dUj d(UjU k ) _ 1 dp 

dt dx k q dx, 


anche se in generale non è più vero che /' =0 e che fg=Jg+f g '. 
Naturalmente si ripresenta anche qui il classico problema della chisura, 
ben noto dalle equazioni di Reynolds, e cioè che i termini u { u k non 
possono essere calcolati direttamente. Questo problema è stato affron¬ 
tato nella Large Eddy Simulation con due diversi approcci, l’approccio 
«implicito» e l’approccio «esplicito». L’approccio implicito, sviluppa¬ 
to in particolare da Lilly [15], Deardorff [16], Schumann [17], si basa 
più genericamente sull’idea che la discretizzazione di per sé attui una 
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implicita operazione di media: si può per esempio notare dalla semplice 
relazione 


(16) 


f(x + Ax)-f(x-àx) 
2Ax 


— (— 
dx \2Ax 


f(x')dx' 


che un operatore discreto di derivazione sulla scala Ax è interpretabile 
come un operatore analitico di derivazione di un valore mediato sulla 
stessa scala. Sulla base di queste considerazioni nell’approccio empirico 
si scrive più semplicemente 

(17) u^u J = u i Uj+T u 

dove la soprallineatura indica il processo di media operato dalla discre- 
tizzazione numerica e Ty i relativi termini di sottogriglia. Nell’approc¬ 
cio esplicito, che è più coerente con quanto detto sinora, e che è stato 
sviluppato in particolare da Léonard [14], Ferziger [18], W.C. Reynolds 
[19] si scrive invece esplicitamente, essendo w, = fi, e w, = fi, = fi y + w/ 

(18) uJUj = Ujùj + (fi,fi,- fi,- fi, ) + «/ fi, + u'j fi, + fi/ w,' 
e si scompone 7/, nei tre termini 

(19) T g . •/.„+C„ + R, y 
dove 

(20) L,, = fi/fi/— fi, fi, 

sono i cosiddetti termini di Léonard, 

(21) C iJ = u' i Uj + Uj Ut 
i termini incrociati, e 

(22) Ry = u 'j uj 

i termini di sottogriglia, e l’operatore di media indicato con la sopralli¬ 
neatura è un filtro esplicitamente dato, tipicamente una Gaussiana di 
ampiezza h 
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(23) 


G (x - x ' ; h) = 



exp 


6(x-x') 2 \ 

h 2 / ‘ 


È chiaro che, mentre i termini di Léonard possono essere calcolati espli¬ 
citamente, sia i termini di sottogriglia fg dell’approccio esplicito che i 
termini Cg e Rg di quello implicito vanno parametrizzati. A tale scopo 
si sono sviluppati i vari modelli di sottogriglia, 


4. Modelli di sottogriglia 

I modelli di sottogriglia adempiono essenzialmente al compito di tra¬ 
sferire energia dalle scale risolte a quelle non risolte, simulando al me¬ 
glio la cascata inerziale e la dissipazione viscosa. Ciò viene realizzato 
nel caso più generale attraverso una parametrizzazione di tipo dissipati- 
vo dei termini di sottogriglia. Nel caso dell’approccio esplicito si pone 


(24) 

Tn - T kk Ò;:= — 2 V t Sjj 

J 3 J 

dove 

(25) 

V 1 (ii 

+ 1 ±) 

2 \ dxj 

dx, ) 


dove v, è una opportuna viscosità turbolenta e dove è opportuno notare 


che il termine — T kk 8 if rappresentante l’energia cinetica di sottogriglia 


viene praticamente conglobato nelle equazioni mediate con il termine 
di pressione senza avere alcun effetto dinamico diretto. Per quel che ri¬ 
guarda la viscosità turbolenta v, la proposta di Smagorinsky [20] è quel¬ 
la di scalarla attraverso la scala di lunghezza Ax di griglia 


(26) v t — (c s Ax) 2 S 

dove S = 'JS ij Sg e dove c s è una opportuna costante. Passando ora al¬ 
l’approccio esplicito occorre innanzitutto ricordare che in esso appare 
esplicitamente il termine di Léonard L it che viene calcolato direttamen¬ 
te e che di per sé già opera un certo trasferimento di energia dalle scale 
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risolte a quelle non risolte. Per i termini C y e R,j dell’approccio esplici¬ 
to due proposte sono state fatte e studiate numericamente: la prima è 
quella di utilizzare anche per essi una parametrizzazione alla Smagorin- 
sky ponendo 


(27) O,, C,j- R 
e 

1 / dUi 

(28) Q u -— Q kk ÒU=-v, (—+ 

dove ora 

(29) v, = (c s h) 2 S 


dGj \ 
dx, I 


con h larghezza caratteristica di filtro. Notiamo che una precisa separa¬ 
zione della numerica dal processo di media richiederebbe h > Ax, ma 
in pratica h = 2 Ax è il valore usato con il filtro Gaussiano. Il principale 
svantaggio del modello dissipativo di Smagorinsky è quello di non tra¬ 
sferire sufficiente energia dalle grandi scale alle piccole scale. Essendo 
in genere nelle Large Eddy Simulation il numero d’onda di taglio 
k c - (Ax) ~ 1 situato nella cascata inerziale la simulazione del processo 
turbolento resta alterata. Un modello di sottogriglia che provvede un 
notevole trasferimento inerziale è quello proposto da Bardina, Ferziger 
e W.C. Reynolds [21] dato dall’espressione 

(30) Cjj^CiXUjU, u.Uj) 

con c b costante numerica, mentre i termini R fJ restano parametrizzati al¬ 
la Smagorinsky. Questo modello, detto modello «misto», è stato recen¬ 
temente provato numericamente [22] al fine di verificare la consistenza 
tra modello di sottogriglia e filtro esplicito usato, ed ha dato buoni ri¬ 
sultati in particolare nel caso di un filtro Gaussiano. 


5. Nuovi sviluppi 

L’approccio esplicito, pur essendo più articolato di quello implicito, 
presenta tuttavia notevoli complicazioni di calcolo e alcuni punti oscu¬ 
ri. Uno di questi è relativo alla invarianza Galileiana dei vari termini 
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in cui è decomposto il tensore degli sforzi turbolenti T i j=u i u j -u i u j di 
cui è immediato verificare detta invarianza. Se infatti operiamo la sosti¬ 
tuzione 

(31) w, = uf + Vj 

con Vj velocità costante di traslazione, si ha Tjj=Tfj. Nell’approccio 
esplicito si pone invece L ij + C ij +Rjj, ed è facile verificare [23] che 
in questa decomposizione i termini 

(32) Ljj = ù~àj- ajUj 
e 

(33) Cjj = Uj Uj + Uj Gj 

non sono, presi singolarmente, invarianti Galileiani. Per i termini L, r 
calcolati esplicitamente, non esistono inconvenienti, ma per i termini Cy 
si pone il problema della coerenza tra i termini parametrizzati e il relati¬ 
vo modello. È interessante notare che nel caso del modello misto in cui 
Cy è dato dalla (30) il modello segue le stesse leggi di trasformazione 
non Galileiana se sottoposto ad una velocità di traslazione uniforme a 
patto che si ponga c b rigorosamente uguale a 1, ed è interessante nota¬ 
re che i test numerici avevano già ottimizzato empiricamente questa co¬ 
stante al valore 1.05, ma tuttavia resta qualche dubbio sulla effettiva 
bontà di questo modo di procedere, e ci si può chiedere se un diverso 
modo di decomporre gli sforzi di sottogriglia non possa portare ad una 
impostazione di fondo più soddisfacente. Notiamo che, dal punto di vi¬ 
sta operatoriale, gli «sforzi turbolenti» t(u ì3 w,) associati ad un opera¬ 
tore di media sono definiti come 

(34) t(m,-, Uj) = UjUj GjQj 

e nell’approccio diretto si ha ovviamente 

(35) T ìj =t(u ì , Uj) . 

Nell’approccio esplicito una decomposizione di r(w„ Uj) coerente con la 
definizione (34) darebbe [24] 

t(u.j, Uj) = t(U„ Uj) + t(u', Uj) + t(U„ u-) + t(u[, Uj) 


(36) 
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dove i nuovi termini 

(37) y.j-Tdi, il,) 
con 

(38) t(u„ Uj) = UjUj - UjUj 
i nuovi termini 

(39) V)j=T(u;,Uj) + T(a i , U-) 
con 

(40) t {uj, Gj ) = uj Uj - u[ ùj 

(41) Tifli, u'j) = Uj Uj-Uj Uj 
e i nuovi termini 

(42) (//,', uj) 
dove 

(43) r(uj ,u-) = u- uj - u- uf 

vanno confrontati con i vecchi Ly C/j e Ry dati dalle (20), (21) e (22), 
ed è evidente che nel caso di operatori di media alla Reynolds ( J:\j = 0 
YJy = 0 e YY jj = u'j uj . Una decomposizione di questo tipo, oltre ad essere 
coerente con la generale definizione di sforzo turbolento data dalla (34), 
produce automaticamente termini di sottogriglia invarianti Galileiana¬ 
mente. È opportuno inoltre notare che il modello misto (30) precedente- 
mente citato viene ad essere interpretato in questa luce più semplicemen¬ 
te come un modello in cui il termine y y viene calcolato esplicitamente, 
mentre il termine C‘ , u = Yy + Yy è modellato alla Smagorinsky. 11 model¬ 
lo misto verrebbe quindi ad essere null’altro che il più semplice model¬ 
lo esplicito compatibile con la nuova decomposizione espressa dalla (36). 

Un altro punto che a parere dell’autore merita qualche approfondi¬ 
mento è quello relativo allo studio delle equazioni filtrate in quanto ta¬ 
li, ed indipendentemente dalla parametrizzazione dei termini di sottogri¬ 
glia. Notiamo che, almeno formalmente, i filtri definiti dalla (11) sono 
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invertibili, e prova ne = è il fatto che / f contrariamente ai filtri di 
Reynolds per i quali/=/. È facile vedere [25] che esistono filtri sem¬ 
plici che applicati alle equazioni di Stokes-Navier producono sforzi di 
Reynolds esplicitamente e facilmente calcolabili, come per esempio nel 
caso del filtro dato dall’operatore integrale 

( 44 ) /(x, t; h)= j /(*' , t) G(x -x'; h) d\x'di ' 
con 

exp 

(45) G(x-x';h)= - 

4 

soluzione fondamentale dell’operatore differenziale inverso 

(46) f-f- h 2 V 2 / 

È chiaro che nei termini della (46) la correlazione ujuj è esprimibile 
esplicitamente come 

(47) UfUj = (Uj- h 2 V 2 «,) (Uj - h 2 V 2 Uj) 
ed essendo sempre dalla (46) 

(48) a i aj=u i aj-h 2 v 2 ùiùj 
si ottiene esplicitamente 


(49) u i Uj = Q i Uj+ 2/? 2 V Uj- V Uj+h A V 2 w, V 2 m ; 

Sottolineiamo il fatto che la (49) rappresenta formalmente una espres¬ 
sione esatta, coerente con il filtro usato e che quindi può essere usata 
indipendentemente da modelli di sottogriglia. Notiamo che la classe di 
filtri di questo tipo è molto grande, e potenzialmente ogni soluzione fon¬ 
damentale di un operatore differenziale spazio-temporale può essere uti¬ 
lizzata alla bisogna [26]. L’uso di questi filtri nelle simulazioni numeriche 
è attualmente oggetto di studio da parte dell’autore. 
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1. Il «Filo» nei programmi spaziali e nelle loro prospettive 


1.1. Scenario attuale 

Il «filo» è stato concepito come «distanziatore» di lunghezza varia¬ 
bile di due masse, che rimangono tuttavia fra loro collegate e con la pos¬ 
sibilità di recuperare interamente una delle due nell’altra. La prima 
applicazione a cui si era pensato riguardava la realizzazione di un siste¬ 
ma che, ospitato in un veicolo orbitante a quote opportune per evitare 
decadimento, potesse essere calato a quote molto inferiori, da cui ope¬ 
rare osservazioni nella atmosfera terrestre e sulla Terra stessa. 

Il «veicolo madre» cui, al momento attuale, fanno riferimento i pro¬ 
grammi avviati è lo «Space Shuttle», ovviamente. La fig. 1.1.1 presenta 
una vista schematica di un tale tipo di impiego, ivi compreso uno dei 
modi previsti per lo sviluppo e la retrazione del filo, che prevede un brac¬ 
cio rigido orientabile, con cui dare inizio e terminare le operazioni. 

Ci si rese immediatamente conto che un satellite può essere posizio¬ 
nato mediante il filo, tanto verso la Terra quanto verso lo spazio e furo¬ 
no quindi concepite varie missioni, aventi scopi anche molto diversi da 
quelli che avevano condotto originariamente all’idea. 

In fig. 1.1.2 è rappresentato schematicamente l’assetto di un sistema 
destinato ad una «missione elettrodinamica». La missione è prevista per 
il 1991 e compirà studi di natura elettrodinamica sulla ionosfera e sul 
campo magnetico della Terra. In particolare si indagherà sulla possibi¬ 
lità di generare potenza per interazione del filo (conduttore) con il cam¬ 
po magnetico. Il filo di 20 km di lunghezza sarà usatq in missioni di circa 
36 ore. 

Un’altra missione basata sullo Space-Shuttle è quella «Atmosferica», 
fig. 1.1.3. Si svolgerà due anni dopo la precedente con un filo di ben 
130 km, allo scopo di studiare l’atmosfera e gli effetti aerodinamici a 
quote in precedenza raggiungibili solo con razzi. 

Studiando la dinamica di un sistema, composto da due masse colle¬ 
gate da un filo, nel campo gravitazionale, appare evidente che il filo può 
anche essere utilizzato per ingenerare nel sistema moti (oscillatori) du¬ 
rante i quali ciascuna delle masse è dotata di quantità di moto variabile 
nel tempo. In particolare tale quantità di moto può oscillare rispetto ad 
un valor medio corrispondente alle condizioni orbitali. Liberando quin¬ 
di dal loro vincolo a filo le due masse in un istante particolare, esse pos¬ 
sono trovarsi ciascuna su un’orbita che consenta particolari risultati (ad 
es. un trasferimento). È questo il caso del programma relativo al Siste¬ 
ma di Rientro dallo Spazio Iniziato con il «Tether». (T.I.S.R.S.), pre- 
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sentato in fig. 1.1.4. La capsula che deve rientrare oscilla rispetto allo 
«Shuttle». Liberata in un punto di minima velocità, si trova automati¬ 
camente su un’orbita di trasferimento che la conduce vicino alla Terra 
(lo «Shuttle», data la differenza fra le masse, subisce variazioni di orbi¬ 
ta modesti). Il frenamento aerodinamico e eventuali paracadute faran¬ 
no poi il resto per il rientro e l’atterraggio. 



Fig. 1.1.1 - Tipologia in una missione basata sullo SHUTTLE. (Per cortesia 
di AERITALIA G.S.S.). 
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Fig. 1.1.2 - Sistema di Satellite a Filo (TSS). Missione Elettrodinamica. (Per 
cortesia di AERITALIA G.S.S.). 
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Fig. 1.1.3 - Sistema di Satellite a Filo (TSS). Missione atmosferica..(Per corte¬ 
sia di AERITALIA G.S.S.). 




66 


ETTORE ANTONA 



Fig. 1.1.4 - Sistema di Rientro Spaziale con Inizio a mezzo di Fili (TISRS). (Per 
cortesia di AERITALIA G.S.S.). 


1.2. Scenario futuro 

Per quanto attiene a programmi futuri, in fase di più o meno avanza¬ 
ta analisi, molte sono le applicazioni previste per il «filo». È sufficiente 
citarne alcune per rendere l’idea dell’interesse dei programmi in cui esse 
sono inserite nonché della varietà di problemi che è possibile risolvere. 

In fig. 1.2.1 è descritta una applicazione multipla per il potenziamen¬ 
to delle capacità di una stazione spaziale a lungo termine. 

Il filo, qui azionato sia verso l’alto che verso il basso, consente: 

— l’immagazzinamento di riserve di combustibile, 

— il rifornimento di combustibile, 

— il mantenimento di una piattaforma per scopi scientifici e 

— la realizzazione di un «elevatore» di comunicazione fra la stazione e la 
piattaforma, con il quale fare manutenzione, realizzare microgravità 
modulabile, posizionare una capsula di rientro, e altri programmi. 
La fig. 1.2.2 presenta una applicazione nell’ambito di una missione 

diretta ad altro pianeta (es. Marte). Il modulo per l’equipaggio, contro- 
bilanciato da masse unite con il «filo», costituisce un «manubrio» che, 
dotato di rotazione, simula gravità di entità variabile con la lunghezza 
del filo e la velocità angolare. 

Altre applicazioni, su cui non ci si sofferma, riguardano ad esempio 
l’effettuazione di missioni destinate a riportare a Terra campioni di «co- 
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mete», nelle quali il filo consente ad esempio l’impatto con la cometa 
di un dispositivo atto a prelevare il campione. 



Fig. 1.2.1 - Applicazione multipla e potenziamento delle capacità di stazione 
spaziale a lungo termine. (Per cortesia di AERITALIA G.S.S.). 
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Fig. 1.2.2 - Gravità artificiale per missione su Marte. (Per cortesia di AERI- 
TALIA G.S.S.). 


2. Aspetti meccanici nei progetti dei satelliti a filo 
2.1. Discipline coinvolte 

Un esame anche superficiale delle configurazioni e delle manovre ot¬ 
tenibili con il «filo» quale distanziatore (flessibile) di masse in program¬ 
mi spaziali conduce alla constatazione che molte discipline devono esservi 
coinvolte. Fra queste in particolare la Meccanica (anche celeste), la Flui¬ 
dodinamica, l’Aerotermodinamica, l’Elettrodinamica, la Geofisica e la 
Fisica del sistema solare. 
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2.2. Principali obiettivi di natura meccanica delle missioni 

Restrigendo il campo agli aspetti meccanici, fra gli obiettivi cospiranti 
al compimento di missioni di satelliti a filo, appaiono importanti i se¬ 
guenti: l’eliminazione o la riduzione del decadimento in orbita bassa, 
la creazione di ambienti mocrogravitazionali di vario livello fino alla gra¬ 
vità artificiale, e, infine, i trasferimenti orbitali. 


2.3. Schemi per Vapproccio dei problemi di meccanica 

Uno degli aspetti della presente nota è quello di considerare, accanto 
al più frequentemente adottato modello discreto: (il filo non ha massa 
ed elasticità distribuite), anche un più sofisticato modello continuo: (il 
filo ha massa ed elasticità distribuite). 


2.4. Scelta di un adeguato sistema di riferimento 

Occorre distinguere fra il caso in cui il complesso «veicolo madre»- 
«satellite a filo» opera nelle immediate vicinanze di un pianeta (es. Ter¬ 
ra) e quello in cui tale complesso si trovi più genericamente nel sistema 
solare. 

Nel primo caso, poiché è trascurabile l’effetto degli altri pianeti, co¬ 
stituisce riferimento corretto (e inerziale per quanto concerne le accele¬ 
razioni lineari, una terna della quale si tenga conto della velocità angolare, 
con un asse coincidente con la congiungente Sole, Pianeta, e due assi 
a questo perpendicolari, uno nel piano dell’eclittica, , nel quale sistema 
operi l’attrazione del pianeta. 

Il gradiente della attrazione del Sole può essere trascurato se ci si tro¬ 
va ad una distanza dal pianeta molto inferiore di quella del primo pun¬ 
to di Lagrange. 

Nel secondo caso, se l’influenza di ogni pianeta può essere trascura¬ 
ta, sarà sufficiente analizzare il moto del sistema «veicolo madre»- 
«satellite a filo» nel suo complesso come un unico corpo attratto dal 
Sole e la dinamica del satellite rispetto al veicolo rispetto ad un sistema 
con un asse sulla congiungente Sole-sistema e gli altri rispettivamente 
tangente e normale alla traiettoria del sistema. 

Tali risultati, per altro noti, possono agevolmente desumersi da quanto 
ad es. esposto in [5], para. 3, nonché, con qualche accorgimento, da quan¬ 
to esposto anche nel presente lavoro al para. 3.1. 
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3. Problemi di meccanica nel modello discreto 

3.1. Dinamica del satellite a filo 

3.1.1. Derivazioni usuali 

Il veicolo madre e il satellite, fra loro tenuti a distanza obbligata al 
filo, costituiscono un «manubrio», il cui baricentro descrive un’orbita. 
Nei comuni approcci si analizza la dinamica di ciascuna delle due masse 
del «manubrio», ponendo nel baricentro l’origine di una terna di assi 
ortogonali, dotata di velocità angolare attorno al proprio asse perpen¬ 
dicolare all’orbita,. Non si considerano gli effetti delle accelerazioni li¬ 
neari del baricentro e si trascurano le mutue attrazioni Newtoniane delle 
due masse. 

La fig. 3.1.1 introduce all’equilibrio fra forza gravitazionale e forza 
centrifuga nel caso di orbite circolari. Se Me m sono le masse del piane¬ 
ta e del sistema «veicolo-satellite» di volume V, si ha: 

-j ^^-rdm-V | u 2 rdm = 0, 

dove G è la costante di gravitazione e r è il raggio (vettore) di ogni dm. 
Se r 0 è il raggio del baricentro, si ha: 

Fg = Fe , 

ovvero 



e anche 



Se una massa dm non è situata nel baricentro ed è dotata di velocità 
rispetto alla terna d’assi, essa risulta soggetta a forze (di d’Alembert), 
che vengono esaminate nelle figg. 3.1.2, 3.1.3 e 3.1.4. 
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Fig. 3.1.1 - Sistema orbitale-traiettoria circolare. Equilibrio per forze gravita¬ 
zionali F g forze centrifughe F c . 
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Nella fig. 3.1.2 si introducono le forze dovute a gradiente di gravità 
verticale. Nel piano orbitale z—x, componenti z = r-r 0 di posizione ri¬ 
spetto al baricentro generano 

1^1 * \F C \ , 

e quindi una componente di forze 


F z = 


/ dF g dF c \ 
\ dr dr ) 


z . 



Fig. 3.1.2 - Gradiente di gravità verticale. 
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' (7 f g - t 2 r 

Il gradiente gravitazione e-— = 2 dmc * k ' 

dF c j 2 r 

Il gradiente centrifugo e —^—=dm co k. 

Risulta 

F z = 3 dm v 2 ìc. 

Nella fig. 3.1.3 si analizza il gradiente di gravità ortogonale al piano 
orbitale. Si considerano componenti Y di posizione rispetto al baricen¬ 
tro. Esse generano 



Fig. 3.1.3 - Gradiente di gravità ortogonale al piano orbitale. 
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Ne risulta una forza, secondo y 
F y =-F g + F c , 

la quale può essere espressa, come risulta dalla figura, nei due modi se¬ 
guenti 


y ^ y — 

F y = F g —j ; F y =-F c — j. 

r r Q 


Ne risulta: 


F y = -dm c o 2 y j . 

Nella fig. 3.1.4 si analizzano le forze dipendenti dalla velocità rispet¬ 
to alla terna di riferimento, dovute all’effetto di Coriolis. La combina¬ 
zione di x con a; genera: 

A2 dm wxk . 



Fig. 3.1.4 - Forze dipendenti dalla velocità. 
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La combinazione di z con ai genera: 

AF=2dm uzT. 

Tutti gli effetti sopra analizzati vengono introdotti (fig. 3.1.5) nelle 
equazioni di equilibrio di una massa dm : 

F e =dm\(x+2wz)T+ly+ù> 2 y)J+(z~2o}x-'5& 2 z)ìc] , 

dove F e è la risultante di tutte le forze esterne applicate a dm (compre¬ 
se le vincolari). 



Fig. 3.1.5 - Equazioni del moto di un sistema spaziale a filo. 
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3.1.2 Approccio con il problema di tre corpi 

Un altro approccio che si avvale sostanzialmente deH’impostazione 
del problema dei 3 corpi semplificato, che si propone qui all’attenzione 
del lettore, appare assai più pregnante e capace di dar conto di tutte le 
scelte operate nell’approccio precedentemente descritto. 

L’asse costituito dalla congiungente il Sole con il pianeta è dotato di 
velocità angolare, che risulta piccola e si trascura. Allo stesso modo si 
trascurano i termini di gradiente di gravità solare, limitandosi a trattare 
casi in cui la distanza fra il veicolo e il pianeta è molto piccola rispetto 
alla distanza del 1 ° punto di Lagrange fra pianeta e Sole. 

Resta, quindi, da considerare la gravità del pianeta. 

Il pianeta, il veicolo e il satellite costituiscono, quindi, un sistema di 
tre corpi, nel quale certamente il satellite non influisce sui restanti. An¬ 
che la massa del veicolo è molto piccola: utilizzando l’approccio dei tre 
corpi semplificato, essa potrà poi essere fatta tendere a zero (La distan¬ 
za del 1 ° punto di Lagrange fra pianeta e veicolo è di pochi millimetri 
dal baricentro del veicolo). 

Le figg. 3.1.6 e 3.1.7 presentano riferimenti e geometrie rispettiva¬ 
mente per moti del satellite nel solo piano orbitale del veicolo e con com¬ 
ponenti anche ad esso normali. 

Su tali basi si può quindi procedere ad una analisi della dinamica del 
satellite a filo. Siano: 

M = massa del pianeta, 

m = massa del veicolo madre, 

£, ri, f = coordinate su assi ruotanti intorno a f con velocità w, 

U = potenziale dell’unità di massa del satellite. 

Le equazioni della dinamica sono le seguenti: 


£-2cot7-g 

2 , dU 

J t= - 

as 

r) + 2 o) £ — c 

2 dU 
° Sr, 

dU 
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Valgono le seguenti espressioni: 


U= , 

r q 


(se G = K 2 = 1), 


M+n 

Z > 3 


/• 2 = £ + —- D +7j 2 ; f 2 = r 2 + f 2 ; 

\ M+ m 


q 2 =U-~ — D) + v 2 ; è 2 =e 2 +r ■ 

\ M+ ir 1 


2 , =M| _ r+ i) +m (_| r+ |) 


si ha per le equazioni della dinamica l’espressione: 


d( 


rj + 2 = 


r o‘ .. . 

Nella espressione di 2<t>, i termini e j- rappresentano gli ef- 

2 2 , 

fetti delle forze centrifughe, ì termini — e — rappresentano invece 

r Q 

gli effetti delle forze gravitazionali. 
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Si effettui ora il seguente combinamento di coordinate: 

d+É' , 

M+m 


V = V * 

r=r • 

Le equazioni della dinamica assumono la forma: 


w - -, 90 

£ — 2 cor? =- , 

a? 


^+2cor = -v, 






Valgono le seguenti relazioni: 

e 2 =r 2 +v 2 ; f 2 =r 2 + f 2 ; 

r 2 =(r+^) 2 + t)' 2 =2ri)+f) 2 +e 2 ; 

LiLLi(i) 2 + i(nM(l 

r DÌ D 2 \ Dì 2 D 2 \D 


2 _ 3 ,, r 1 f 2 


e, quindi, ricordando che —- = co 2 --, si ottiene anche: 

Z ) 3 £) 3 


rj ' 2 2w 


a meno di una costante. 
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Con unità di misura per cui M= 1, co= 1, posto r 
D 3 = m* + 1 




In definitiva le equazioni della dinamica diventano: 




r?'+2 — 

\ Q 1 +n 


v'=o 


f,+ i~T + T^ 


r=o. 


Ove m*^> 0 (m~* 0), si ha: 

r-2v-«r=o, 
i}'+ 2r =0 , 
r+ r=o. 

Il termine - 3 £' nella prima equazione,-già trovato nell’altro approc¬ 
cio, sembra discendere, dalle sole forze gravitazionali. Un piu attento 
esame rivela però che nella espressione di 


r i . 

+ — compaiono ì termini 
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il secondo proveniente dalle forze centrifughe, e gli altri dalle gravita¬ 
zionali. Il secondo ed il terzo termine si elidono, ma se si vogliono gli 
effetti separati della centrifugazione e della gravitazione, si può osser¬ 
vare che a meno di infinitesimi di ordine superiore | q \ = | f ' | e quindi 
l’effetto totale consta di due parti gravitazionali e una parte centrifuga. 

3.1.3. Equilibrio del satellite 

Le relazioni ottenute nel precedente paragrafo sono le accelerazioni 
cui, il satellite è assoggettato rispetto al riferimento £', 7? ', f dipendenti 
dalle posizioni e dalle velocità rispetto allo stesso sistema, per effetto 
delle attrazioni Newtoniane, degli effetti di Coriolis dovuti al moto del 
sistema di riferimento. In effetti il legame rappresentato dal filo intro¬ 
duce forze esterne al satellite e determina certi valori delle componenti 
di accelerazione. 

In fig. 3.1.8 ad esempio si può osservare l’equilibrio. Nel piano 77 ', 
£ ' il satellite è soggetto ad un campo di accelerazione £ ' mentre è tratte¬ 
nuto dal filo. Nel piano f, £' esso è soggetto ad un campo di accelera¬ 
zioni £' e f, mentre è trattenuto dal filo. 

La fig. 3.1.9 rappresenta una via intuitiva per rendersi conto che nel 
moto di librazione nel piano trasversale le forze di richiamo sono piu 
intense che nel piano orbitale. Un satellite, disturbato longitudinalmen¬ 
te in posizione, procede sulla traiettoria, sia pure anticipando o ritar¬ 
dando. Se disturbato trasversalmente si pone su una traiettoria che rispetto 
alla precedente ha moti oscillatori, perché su un piano inclinato rispetto 
a quello della precedente. 

Nella trattazione esposta in 3.1, si è supposto che l’origine del siste¬ 
ma di riferimento £', 77 ', compia attorno alla Terra orbite circolari 
(con velocità angolare costante). Il baricentro di un sistema satellitare 
composto di più parti ma privo di forze esterne applicate può compiere 
una tale orbita (come anche orbite di altra natura geometrica). In quel 
che segue faremo ancora riferimento ad orbite circolari, lasciando al let¬ 
tore di immaginare i necessari ampliamenti dell’analisi atti a compren¬ 
dere altre orbite, e penseremo all’origine del sistema di riferimento come 
al baricentro del sistema orbitante. 

3.2. Moti di librazione 

Il sistema costituito dal «veicolo madre», dal «filo» e dal satellite co¬ 
stituisce un «manubrio», nel quale il filo per effetto dei gradienti gravi¬ 
tazionali risulta teso se 
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Fig. 3.1.9 - Perturbazione nel piano trasversale. 

(v'=o e rr=o 
1^ = 0 t t = 0 

(lo si deduce anche dalle equazioni della dinamica). Tale stato è an¬ 
che stabile, nel riferimento £', r?' eze sono possibili attorno ad es¬ 


so moti detti di librazione, le cui equazioni, posto 0 = arc sin 


/ 
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</> = arc sin -y , sono, (fig. 3.1.10): 

f 6= -3w 2 sin0 cos0 
1 ij)= -4 ù) 2 sin<f> cos (f> 

rispettivamente nel piano orbitale e nel piano normale. 

È interessante notare come, per grandi oscillazioni in 0 e </>, il filo può 
risultare in bando (cessare di essere teso), generando così intuibili com¬ 
portamenti qui non analizzati. Nella fig. 3.1.10 sono riportati da [3], 
i valori della tensione nel filo in funzione di 0 per vari 0 max . Le zone con 
valori minori di zero darebbero luogo al «bando», che si verifica per 
0 max = 65° e <t> max = 60° circa. 


3.3. Trasferimento di quantità di moto 

Il filo è un dispositivo capace di provocare trasferimenti di quantità 
di moto angolari, nel moto rispetto al pianeta, fra le due masse del ma¬ 
nubrio. 

1) TRASFERIMENTO per cambiamento di lunghezza 

Il distanziamento di due masse lungo il filo, attuato a 0 = 0, ovvia¬ 
mente, riduce la q.d.m. angolare della massa che si avvicina al pianeta 
e, di altrettanto, aumenta l’altra. 

2) TRASFERIMENTO per oscillazione 

Durante una oscillazione, nel piano orbitale, la q.d.m. delle due masse 
cambia continuamente, a somma costante, realizzandosi un trasferimento 
periodico di q.d.m. fra le due masse. 

3) TRASFERIMENTO per rotazione 

La rotazione del manubrio provoca anch’essa un trasferimento pe¬ 
riodico di q.d.m. fra le due masse. 


3.4. Gradiente gra vitazionale 

Si consideri il sistema satellitare nella configurazione stabile r? ' = 0, 
r)' = 0, J - = 0, f = 0. Ogni punto di £, lungo il filo, è sede di un valore 
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Fig. 3.1.10 - Librazioni di un manubrio in orbita circolare. (Da bibl. [3]). 


locale di gravità-apparente, la quale, come già evidenziato, ha origine 
da due contributi: 

1) la diversa attrazione gravitazione a g del pianeta rispetto alla quota 
del baricentro, 

2) la diversa accelerazione centripeta a c , a parità di velocità angolare, 
rispetto alla quota del baricentro. 

Alla quota del baricentro (£ = 0) l’attrazione del pianeta è esattamen¬ 
te uguale alla accelerazione centripeta. Aìtrove l’accelerazione apparente 
è la differenza a a = a g — a c (le accelerazioni si intendono positive se di¬ 
rette verso M) 

a=GM - l —— = Du 2 M-2-) 

g (D + O 2 \ D) 

« c = u 2 (Z> + £) = u 2 ^1 + ^ 


essendo GM/D 2 = w 2 D. 
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Quindi, 



Poiché £ può essere compreso fra zero e valori in modulo molto pic¬ 
coli rispetto a D, possono essere realizzate lungo £ stazioni a «microgra¬ 
vità», molto utile per vari esperimenti scientifici. Stazioni a gravità ridotta 
possono anche servire per l’immagazzinamento di combustibili o altri 
materiali. 


4. Principali operazioni da attuarsi per l’effettuazione delle missioni 
4.1. Svolgimento e retrazione 

Si consideri ancora una massa dm del satellite appartenente al siste¬ 
ma veicolo-satellite a filo, collocate a distanze / dal baricentro di detto 
sistema, ove hanno origine £' e r?', ad una anomalia 6 rispetto a £' in 
senso antiorario (fig. 4.1.1). 

Ci si limita qui a verificare l’esistenza di particolari traiettorie di dm, 
ed esattamente traiettorie a 6 costante. 

Valgono le seguenti relazioni geometriche e cinematiche 

f £' =/ cos 6 r k*~l cos# 

[ ij' -/ sin# , l v' =1 sin# • 


Valgono inoltre le seguenti espressioni delle componenti di forze di 
gravitazione e complementari applicate a dm (nel piano £', r/'): 

(F k ' = 3 dm oi 2 £' +2 dm on)' ; 

(/y= —2 dm o>£' . 
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4.2. Strategie di controllo della configurazione 

Durante le manovre di svolgimento e rotazione e durante la fase di 
impiego, possono manifestarsi oscillazioni di varia natura. Alcune di esse, 
principalmente nel piano orbitale, possono essere attenuate con oppor¬ 
tune «strategie». L’azione di controllo può svolgersi secondo varie for¬ 
me, tra cui: 

— tensione nel filo, 

— spinta, 

— generatori di spinta, 

— masse mobili, 

— «tether» rigido, braccio mobile, 

— aerodinamica. 

Un esame dettagliato di ciascuno di questi mezzi e delle sue possibili¬ 
tà esce dallo scopo della presente analisi, e per esso si rimanda alla lette¬ 
ratura. Le azioni esercitabili attraverso la tensione nel filo sono invece 
spiegabili con gli stessi fenomeni dinamici fin qui illustrati. Le oscilla¬ 
zioni nel piano orbitale (librazioni) possono essere smorzate rilasciando 
ulteriormente il filo mentre il filo stesso è nella fase regressiva e richia¬ 
mandolo quando è in fase progressiva. Le oscillazioni longitudinali so¬ 
no smorzate dalla dissipazione di energie nel filo. 

Alla dm è applicata inoltre la trazione del filo. Essa ha componente 
nulla perpendicolarmente al filo stesso (e se si trascura la differenza fra 
C.G. e punto di vincolo del filo perpendicolarmente a / ). In tale dire¬ 
zione l’equilibrio è quindi dato da 

F r sin0-F v ocs0 = O . 

Tenuto conto delle espressioni di e F,-, delle relazioni geometriche 
e di quelle cinematiche, si ottiene 


Posto 


— = — — u sin# cos# . 
/ 2 


— cosinflcos#, esiste la soluzione 
2 


l=l 0 

con a > 0 nel I e III quadrante (traiettorie di retrazione) e a < 0 nel II 
e IV quadrante (traiettorie di svolgimento). 
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La componente della direzione del filo della forza di gravitazione è 
F / = i r r cos# + /y sin# 
e sostituendo: 

F,= 3 dm co 2 / cos 2 # . 

D’altra parte, l’accelerazione di dm è 
d 2 l 9 

-— = a 2 l =— co 2 / sin 2 # cos 2 # = - a / 

dt 2 4 

mentre la sua velocità è 


di 3 

— = - a 1= -co sin# cos# / . 

dt 2 


Indicata con T la trazione nel filo, deve essere 


d 2 l 

iFi~T) =dm —y , 
dt 


T „ , d 2 l 

T=F,-dm —r- . 

dt 2 


In definitiva si ha 


T= ^3 - sin 2 #j dm co 2 cos 2 # / . 

Tale valore è certamente positivo e risulta funzione di #. Se dm si trova 
ad anomalia # con lunghezza l 0 del filo, alla velocità 


3 _ 

2 


sin# cos# l 0 
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e sottoposto alla trazione del filo predetta, la sua traiettoria sarà una 
esponenziale di / nel tempo, a d costante. 

In pratica si realizzano traiettorie (fig. 4.1.1) di cui quelle a 6 costan¬ 
te sono asintoti. 


4.3. Amplificazione e attenuazione di moti 

I moti di librazione'e di «spin» (cioè di rotazione attorno al baricen¬ 
tro del sistema «veicolo satellite a filo», possono essere amplificati (o 
attenuati), agendo attraverso la trazione nel filo. L’attenuazione è un’o¬ 
perazione già accennata al paragrafo precedente, almeno per la libra¬ 
zione. L’amplificazione si ottiene esercitando trazione (e recupero) nel 
filo quando la trazione è massima, e rilasciando quando essa è minima. 
Si ottiene così un effetto di «pompaggio» nel quale l’energia trasmessa 
quando la trazione è più elevata è minore di quella asportata quando 
la trazione è meno elevata. 


5. Alcuni aspetti della dinamica di un modello continuo 


5.1. Concetti generali e problema esempio 

Il modello discreto, con masse concentrate per il veicolo madre e per 
il satellite e con filo privo di massa (distribuita), consente di studiare 
molti aspetti della dinamica del sistema, fino a ben precisi livelli di ap¬ 
prossimazione. 

Se si prende in esame un modello che consideri la massa distribuita 
del filo e la sua elasticità, il panorama dei problemi e di comportamenti 
che possono essere esaminati si allarga enormemente. 

A titolo di introduzione a tale panorama, si riportano qui alcune con¬ 
siderazioni riguardanti i fenomeni dinamici che intervengono nel filo, 
in particolare mentre se ne opera lo svolgimento o la retrazione, mentre 
cioè il suo punto vincolato al veicolo madre cambia nel tempo. 

La dinamica dei «tether» (per i problemi considerati in questo capi¬ 
tolo, la ref. [3] costituisce un valido aiuto), così come quella di ogni al¬ 
tra struttura, può essere considerata in una analisi lineare come la 
sovrapposizine di onde di propagazione. Ciò richiede lo studio della di¬ 
namica della propagazione lungo il tether con le opportune condizioni 
al contorno. Con riferimento ad un sistema di riferimento lagrangiano 
(un sistema, cioè che introduce una corrispondenza biunivoca tra una 
coordinata longitudinale ed ogni sezione del tether), durante il rilascio 
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ed il recupero il tether può essere considerato come composto da una 
parte a riposo e da una soggetta al fenomeno di propagazione. Queste 
due parti sono separate da una sezione che cambia col tempo, cioè la 
sezione del tether che limita la parte vincolata al riposo cambia col tem¬ 
po (la restante parte è libera di muoversi e di vibrare ed ha la sezione 
opposta soggetta alle condizioni al contorno imposte dal satellite). 

I fenomeni di propagazione di maggior interesse da un punto di vista 
pratico sono i seguenti: 

1) Onde trasversali, seguenti un comportamento di corda vibrante, do¬ 
ve sia le forze inerziali che la tensione nel tether — in combinazine 
con la curvatura locale — sono gli elementi più importanti dell’equi- 
librio dinamico. 

2) Onde longitudinali, rappresentate matematicamente dalla equazione 
della corda vibrante, dove le forze inerziali e le forze interne longitu¬ 
dinali dovute alla deformazione elastica, sono gli elementi piu impor¬ 
tanti dell’equilibrio dinamico. 

La dipendenza dal tempo della sezione che vincola la parte a riposo 
richiede l’analisi della riflessione dell’onda che viaggia verso la parte a 
riposo. Tale riflessione genera un’onda che viaggia verso l’esterno e le 
cui caratteristiche devono essere determinate. 

Questo capitolo tratta della necessità di avere una parte esterna a ri¬ 
poso (non solamente le sezioni di contorno). Perciò le sezioni di contor¬ 
no si comportano come delle superfici dove i fenomeni dinamici interni 
subiscono una riflessione. 

Si fa ovviamente riferimento ad una struttura monoassiale di sezione 
costante, entro la quale categoria il tether può essere qualificato. Cari¬ 
chi longitudinali addizionali dovuti alla microgravità permettono di te¬ 
nere conto di piccoli carichi di tensione o di compressione senza incorrere 
in fenomeni critici. 

5.2. Fondamenti fisico-matematici 

Dapprima richiamiamo alcuni concetti con riferimento soltanto al caso 
delle onde longitudinali. Il lettore potrà facilmente estendere l’analisi 
alle onde trasversali. Il problema può essere analizzato per mezzo di una 
equazione di continuità e di una condizione di equilibrio. 

Se <t = o(x, t) è lo stress nel punto x al tempo t e u = u{x, t) è la veloci¬ 
tà, l’equazione linearizzata dell’equilibrio è 


do du 

~dx^ Q di ‘ 


(1) 
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D’altra parte l’equazione linearizzata della continuità è 


( 2 ) 


du _ 1 da 
dx t dt 


Se poniamo c 2 = E/q, otteniamo, 


d') 


d 2 a 1 d 2 
dx 2 c 2 dt 2 


2 U _ i d 2 

dx 2 c 2 dt 2 


Le soluzioni generali delle equazioni (!') e (2') sono 


< 3 > a - e , ('-f) +62 ( ,+ f) 

«4, „.„(,-£) + „( (+ f). 


Le equazioni (3) e (4) indicano che il moto della sbarra è composto da 
due onde: una che viaggia nella direzione delle x positive ed un’altra che 
viaggia nella direzione delle x negative. 

Sulla base di queste considerazioni è possibile ottenere una relazione 
tra p ed u per ogni singola onda. 

Se fi =» t - — e t 2 + 1+ —, da 3) e 4) si ottiene: 


da _ 1 / dp 2 dpi \ 

dx c \ dt 2 dt x / 

du _ du x du 2 

dt ~ dt x dt 2 


(6) 
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/ dp, 

dp 2 \ 


+ ilC| 

l Sii ' 

" at 2 ) 

U<i 

Uh 1 


Considerando l’equazione (1), da (5) e da (6) abbiamo la relazione 

(7) 

Se p 2 = 0, u 2 = 0, l’equazione (7) fornisce la relazione 

dp, dw, _ E du } 

dx t QC dti ~ c a/, 

D’altra parte se p, = 0, u x = 0 abbiamo la relazione 


dp 2 du 2 

- pc —m 
dt 2 dt 2 


E du 2 
c dt 2 


Con le condizioni iniziali u x (t 0 ) = 0, Pi(t 0 ) = 0 o, separatamente, 
p 2 (/ 0 ) = 0, u 2 (t 0 ) = 0 valgono le seguenti relazioni per ogni onda viag¬ 
giante: 

( 8 ) 


(9) 

Tali relazioni ci mettono in grado di determinare uno dei due valori w, 
oppure p, quando ne conosciamo l’altro (equazioni 8) ed inoltre di de¬ 
terminare uno dei valori u 2 oppure p 2 quando ne conosciamo l’altro 
(equazione 9). 

Questi sono rispettivamente in relazione con il comportamento sepa¬ 
rato delle due onde in propagazione. 



5.3. Riflessione sulla sezione dipendente dal tempo 

Le condizioni al contorno di solito considerate nelle applicazioni pra¬ 
tiche per sbarre monoassiali elastiche sono l’estremo libero (o = 0) e l’e¬ 
stremo incastrato (u = 0). Alle estremità dove avviene il rilascio, si verifica 
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che u(x c ) = 0 dove x c = x c (t) è la sezione dipendente dal tempo che si può 
considerare come fissa e dx c /dt = v è la velocità di spostamento del vin¬ 
colo. In questa notazione, x c è una coordinata della barra non de¬ 
formata. 

Per una prima analisi del comportamento di una sezione fissa dipen¬ 
dente dal tempo possiamo considerare il problema di una barra elastica 
che ha un estremo libero ad x = 0 ed un estremo vincolato con u(x c ) = 0 
ad x c = l 0 + cv, dove l 0 sono costanti. La sbarra è soggetta in x = 0 ad 
una forza specifica /, estensionale esterna indipendente dal tempo t. 

Una tale forza produce un’onda che viaggia nella direzione delle x 
positive; quando tale onda raggiunge la sezione x c di vincolo mobile vie¬ 
ne generata un’onda riflessa di forza specifica p 2 che viaggia nella di¬ 
rezione delle x negative. 

Determineremo ora le caratteristiche dell’onda riflessa, prima che es¬ 
sa raggiunga la sezione x = 0. 

Allo scopo di determinare l’onda riflessa per mezzo di un bilancio 
di energia è necessario disporre di una valutazione dell’energia scam¬ 
biata all’estremo vincolato. Il vincolo ha un bilancio di energia diverso 
da zero e l’onda riflessa può essere determinata come segue. È già stato 
introdotto lo spostamento sezione dove è imposto il vincolo durante il 
tempo. Tale spostamento corrisponde alla deformazione interna della 
sbarra. 

L’intervallo dx=vdt durante dt, è sottoposto ad una variazione di 
lunghezza di- - («, + u 2 )dt con un allungamento E = dl/vdt. 

D’altra parte la tensione finale in vdt deve essere o=p x +p 2 = 
= (£/, - U 2 )E/c il che significa un allungamento e — (u ì -u 2 )/c. Ugua¬ 
gliando le due espressioni dell’allungamento si ottiene: 

c-v 

(10) u 2 =-u, -. 

C+ V 

Questo risultato coincide con quello dell’applicazione della soluzione del 
problema di Goursat, Daboux e Cauchy, con la condizione u = 0 in 
x=l„ + vt. Oltre alla coincidenza con i ben noti risultati della analisi del¬ 
la corda vibrante, il modello proposto presenta il vantaggio di essere ap¬ 
plicabile a problemi più complessi quali i sistemi dispersivi. 
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5.4. Estensione ad una velocità di propagazione dipendente 
dalla posizione longitudinale 

Il presente paragrafo riguarda l’estensione dell’analisi al caso in cui 
le equazioni (1 ') e (2') diventano come segue 


0 ") 


d 2 a 1 d 2 a 
~dx 2 ~ C 2 (X) ~W 


( 2 ") 


3 2 u _ 1 d 2 u 

dx 2 c 2 (x ) dt 2 


Questo caso comprende la dinamica di onde trasversali in tethers do¬ 
ve il carico longitudinale dipende dalla posizione per effetto della mi¬ 
crogravità. Siccome c(x) non è costante ma una funzione di x, c = c(x) 
le espressioni (3) e (4) non sono più valide, ed è necessario trovare un 
modo appropriato di soluzione. In un articolo precedente, [6], l’autore 
proponeva che la soluzione generale delle equazioni (1 ") fosse composta 
da due onde che viaggiassero in opposte direzioni con velocità dipen¬ 
dente da x. Tali onde si riducono alle equazioni (3) e (4) quando c(x) 
diventa una costante C. se <t> indica a oppure u quando necessario, si 
può adottare la seguente espressione (quando c(x) = costante e K= 1 l’e¬ 
spressione proposta si riduce a (3) e a (4)): 


OD 


(,_ 1 

[ dx \ 

1 +02 | 

LA 

! dx \ 

l ! 

| Kc(x)) 

r! 

I Kc(x)J 


dove K è una costante arbitraria e ±Kc(x) potrebbe essere considerata 
come la velocità di propagazione in x. Ovviamente, non tutte le funzio¬ 
ni 0, e <t> 2 sono utili per soddisfare la (1 "). Il problema si riduce ora al¬ 
la determinazione delle funzioni 0, e 0 2 , qualora dovessero esisterne, 
che soddisfino alla (1 "). Le funzioni 0, e 0 2 possono essere esaminate 
separatamente. 

f dx 

Se poniamo s = t- \ -, si ottiene 

] Kc(x) 
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<30 ! 

d<j) { 

ds 

_ d<t>i | 

( 1 

\ 

dx 

ds 

dx 

“ ds 1 

V Kc{x) 

1 

d 2 fa 

d 2 4> 

i ! 

1 \ 

2 + dfa! 

1 c‘(x) \ / 

dx 2 

ds : 

' [ 

Kc(x)) 

ds \ 

K c\x) ) ’ V 


d<j> i 

ds 

dfa 

9 2 0i 

d 2 (t> i 

dt 

ds 

~dt 

ds 

dt 2 

ds 2 


Richiamando ora la (1"), otteniamo 


(12) 


d2( t» 1 , dfa 1 c , d %. _ Q 

tfs 2 tf 2 dsK ds 2 


La funzione 0,, se esiste, deve soddisfare l’equazione (12). Siamo por¬ 
tati alla stessa conclusione per quanto riguarda fa. 

Come esempio, possiamo restringere l’analisi al caso c = c Q (x). Ab¬ 
biamo che c‘(x) = c 0 e che 


( 12 ') 



Tìl^!=o. 

K ds 


Questa è una equazione differenziale lineare del secondo ordine. 
La sua soluzione è del tipo fa = e ds e, più precisamente 


(13) 4>ì = exp j^+ -—— (-KC 0 + ln x)j 

Considerando ora fa, otteniamo che 
(13') fa = exp U ( KC 0 + In *)j . 


Se lasciamo che k assuma tutti i valori da zero all’infinito otteniamo per 
</>, e per fa un insieme completo di funzioni, che ci permettono di esten- 
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dere il risultato a qualunque funzione tramite integrazione. Tutti i feno¬ 
meni dinamici in una struttura dove le equazioni (1") sono valide e 
c = CgX possono essere analizzati da (si vedano le equazioni (13) e (13')), 

4>=j C] (K )• exp |+ i l - (-KC 0 + lnx)^dk + 

+ | C 2 (K) ■ exp 1 1 (KC 0 + Inx )j dx. 


Ovviamente, l’equazione (14) non rappresenta l’unico modo tramite cui 
espandere per mezzo di un integrazione un tale fenomeno dinamico, ma 
questo modo di procedere ci permette di considerare delle funzioni com¬ 
ponenti quando la velocità di propagazione di ogni x è nota. Come sap¬ 
piamo, tale velocità è un dato fondamentale per la valutazione dell’onda 
riflessa, in particolare nel caso di condizioni di vincoli dipendenti dal 
tempo. 
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Modelli e pregiudizi 
nella meccanica dei fluidi odierna 


Gino MORETTI* 


Ogni tentativo di correlazione generale fra entità fisiche è un model¬ 
lo. L’analisi matematica ha permesso di cristallizzare i modelli in espres¬ 
sioni sintetiche, alcune delle quali brillano per semplicità (e quindi per 
forza espressiva) come la classica equazione F = ma e la popolarissima 
E = mc 2 . 

Nessun modello, tuttavia, è corretto. Il vero, il giusto, sono asintoti 
ai quali tentiamo di avvicinarsi, senza sapere di quanto ne distiamo. Ta¬ 
lora ci accorgiamo che un modello proprio non va, e lo releghiamo nel¬ 
la categoria dei pregiudizi. Così passò dopo tanti secoli la teoria 
aristotelica di proporzionalità tra velocità e forza; così dopo alcuni anni 
la credenza dell’isentropia degli urti. 

Si ha quindi, lavorando su un modello, una lenta evoluzione di pensiero 
e un crescente disagio intellettuale, finché di colpo un nuovo modello vie¬ 
ne a sostituirsi al vecchio; la rivoluzione è tanto più benefica e duratura 
quanto più semplice è l’espressione del nuovo modello. Si pensi alla rivolu¬ 
zione copernicana, al trionfo delle coniche su cerchi e sulle epicicloidi. 

In quella rivoluzione c’è anche molto di più. Infatti, il sistema plane¬ 
tario come lo vediamo noi (o almeno, come l’abbiamo visto fino ad al¬ 
cuni anni fa) è «spiegato» dal modello in termini di cause ed effetti, 
mentre il sistema tolemaico è puramente descrittivo. 

È bene notare che: 1) l’insoddisfazione crescente per un modello, alme¬ 
no dal punto di vista dei fisici, può essere solo generata da nuove osserva¬ 
zioni empiriche, 2) fintantoché si rimane nell’ambito del modello prescel¬ 
to, ogni nuovo risultato dev’essere una conseguenza logica delle premesse. 
Questo secondo punto è essenzialmente il tema del discorso che segue. 

Nell’era classica della meccanica dei fluidi, i modelli erano pochi e 
di grande bellezza, e i loro limiti erano chiaramente stabiliti. Per esem¬ 
pio, la teoria dei moti stazionari dei fluidi incompressibili non viscosi 
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insegnava che le velocità agli spigoli erano nulle o infinite, e che un moto, 
uniforme all’infinito, intorno a qualunque figura chiusa, era definito 
meno di un vortice di intensità arbitraria. Queste erano ovvie violazioni 
del buon senso fisico, ma erano accettate come necessarie conseguenze 
del modello prescelto. Tuttavia l’insoddisfazione dianzi menzionata ge¬ 
nerò due nuovi modelli alla fine del secolo scorso, a integrare il modello 
precedente. Uno diede luogo ai campi di moto a pelo libero, brillante- 
mente descritti da Kirchhoff e purtroppo alquanto dimenticati. L’altro 
fu la geniale invenzione di Kutta e Joukowski che è alla base di ogni 
studio di profili alari. 

Le limitazioni del modello incompressibile non viscoso, dal punto di 
vista pratico, sono molte e non sempre gradevoli. La viscosità naturale 
si cura di eliminare le singolarità agli spigoli, e di conseguenza anche 
di generare campi di moto simili a flussi di Kirchhoff e di giustificare 
il modello di bordo di fuga di Kutta e Joukowski. L’importanza di que¬ 
sti due modelli risiede nella loro capacità di presentare certe soluzioni 
irrotazionali come possibili casi limite di moti viscosi a numero di Rey¬ 
nolds infinito. 

In questi esempi, le equazioni indefinite del moto sono rese definite 
da un’opportuna scelta di condizioni al contorno e, per ciascuna scelta, 
la soluzione può essere ottenuta per via analitica inequivocabilmente, 
sicché la natura stessa delle condizioni al contorno e la dipendenza da 
esse della soluzione non possono essere oggetto di controversie. 

Un modello totalmente diverso è alla base della meccanica dei fluidi 
compressibili. La propagazione istantanea, tipica dei fluidi incompres¬ 
sibili, è sostituita da una propagazione a velocitò finita, e questo getta 
una nuova luce sui rapporti tra moti stazionari e non. La densità non 
è costante; compressioni ed espansioni si sviluppano con tendenze op¬ 
poste e le compressioni, presto o tardi, diventano incompatibili con una 
distribuzione di valori differenziale. Modelli di stati stazionari si posso¬ 
no concepire per fluidi compressibili, ma sembrano molto più restrittivi 
che per fluidi incompressibili. Infatti, per un fluido incompressibile non 
viscoso, la non stazionarietà è una successione di stati stazionari. Al con¬ 
trario, nel modello compressibile, la non-stazionarietà sembra essere lo 
stato naturale, con stati stazionari solo raggiungibili in forma asintoti¬ 
ca. Le equazioni di Eulero per stati stazionari sono un modello; il cre¬ 
dere alla loro validità fisica è un pregiudizio. 

Nell’ambito dei fluidi compressibili, troviamo un gran numero di mo¬ 
delli più ristretti: moti monodimensionali, quasi-monodimensionali, bi¬ 
dimensionali, subsonici, supersonici, transonici, con o senza discontinuità. 
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In quei rari casi in cui si può ricavare una soluzione analitica, le condi¬ 
zioni al contorno variano radicalmente da un modello all’altro e, se il 
problema non è stazionario, condizioni iniziali intervengono in modo 
essenziale. 

Ancora una volta, se il trattamento del problema è analitico, esso ri¬ 
mane fedele al modello; le condizioni iniziali e al contorno sono scelte 
di necessità in accordo con il modello e la soluzione del problema ne 
consegue coerentemente. Purtroppo, soluzioni analitiche in campi com¬ 
pressibili, anche non viscosi, si trovano solo per rari problemi di grande 
semplicità e servono essenzialmente di guida all’interpretazione fisica dei 
fenomeni più complessi. L’analisi numerica è ormai accettata come me¬ 
todo di studio di problemi della vita reale. Anche qui, tuttavia, si deve 
scegliere un modello di base. Anche qui bisogna scegliere un insieme di 
condizioni iniziali e di contorno. 

La difficoltà nuova proviene dalla discretizzazione necessaria per usare 
il calcolatore. Nel fare ciò, è facile fuorviare il modello o le condizioni 
iniziali e di contorno. È compito dell’analista di far sì che ciò non av¬ 
venga. L’analista deve quindi al tempo stesso conoscere la fisica del pro¬ 
blema e i metodi numerici: anzi, io credo, dev’essere più un fisico che 
un matematico! 

È doloroso constatare che, ai nostri tempi, il calcolo tende a rimpiaz¬ 
zare il pensiero. La massificazione delle scienze, stimolata dal consumi¬ 
smo tecnologico e accelerata dallo sviluppo delle macchine calcolatrici 
e, d’altro canto, la furia competitiva hanno messo al bando lo studio 
della filosofia della natura. Comunicare è ormai pressoché impossibile, 
paradossalmente, perché i mezzi di comunicazione sono troppi. In pra¬ 
tica, tutti scrivono e nessuno legge. Le rare gemme sono destinate a ri¬ 
manere nascoste. I lavori citati (ma non necessariamente letti, e tanto 
meno digeriti) provengono dalle scorie superficiali della letteratura tec¬ 
nica. Inoltre è difficile comunicare per varie ragioni, di volta in volta 
o tutt’insieme presenti: L’autore, per troppa familiarità col soggetto, sor¬ 
vola su informazioni importanti; l’autore non sa scrivere o scrive in una 
lingua straniera in modo approssimativo ed arruffato; l’autore crede di 
aver capito il suo problema ma non l’ha sufficientemente approfondito 
(quanti di noi hanno il coraggio di riconoscere i propri limiti?); l’autore 
è volutamente oscuro per timore di plagi! E infine, un discorso di qua¬ 
ranta pagine dev’essere compresso in otto per motivi editoriali. D’altro 
canto, ai giovani immaturi ed inesperti spesso manca, e non per colpa 
loro, la guida solerte del docente. Responsabile ne è la rilassatezza di 
costumi nelle scuole di ogni rango dove, per timore di tarpare le ah a 
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geni in potenza, si trascurano i puntini sulle i, si dà troppa importanza 
alle definizioni formali che sono difficili da interpretare e, paradossal¬ 
mente, si stimola una generale tendenza ad accettare luoghi comuni, in¬ 
vece di forzare gli studenti aH’autocritica. Di conseguenza, certi modelli 
diventano facilmente pregiudizi. 

Che cosa è infatti un pregiudizio? Un pregiudizio è un modello cre¬ 
duto valido nel seno di un altro, quando non lo è, oppure un modello 
a validità limitata, esteso arbitrariamente al di fuori del suo campo di 
validità. Ecco alcuni esempi. 

1) Un’idea generale, abbastanza corretta, dei moti di fluidi compres¬ 
sibili, stazionari e non, in condotti ad area variabile si ottiene attraverso 
calcoli quasi-monodimensionali. Il modello è, anzitutto, limitato a va¬ 
riazioni di area dolci. Inoltre è assai poco rappresentativo nel caso di 
flussi supersonici, dove il gioco delle caratteristiche e di onde d’urto obli¬ 
que è essenzialmente multi dimensionale. Il pregiudizio di «potersi fare 
qualche idea» su come il condotto si possa reagire a stimoli esterni con¬ 
duce soltanto a «qualche idea sbagliata». Pregiudizi di sorta non sono 
solo errori di principianti. Il pubblico denaro di molti paesi ricchi è so¬ 
vente sprecato in varie imprese del genere. 

2) Un’altra fonte comune di pregiudizi è la confusione tra fenomeni 
stazionari e non. Lasciamo da parte errori volgari, come la credenza nella 
costanza di valori d’arresto in moti non stazionari. In una sfera supe¬ 
riore c’è ben di più. I modelli stazionari, che in un primo tempo aveva¬ 
no permesso di semplificare le equazioni del moto e ottenere alcune 
soluzioni in forma chiusa, rimasero come paradigma anche quando l’av¬ 
vento dei computers suggerì l’uso di tecniche numeriche. La struttura 
stessa delle macchine impone formulazioni iterative sicché anche un fe¬ 
nomeno stazionario è numericamente ottenuto attraverso un’evoluzio¬ 
ne temporale, sia essa interpretabile fisicamente o no. In entrambi i casi, 
il cammino è infestato da trappole che il più delle volte sono ostiche a 
chi è condizionato da pregiudizi di tipo stazionario, anche se sono facili 
da capire in termini di meccanica dei fluidi classica. 

Considerando ogni calcolo numerico come non stazionario, la prima 
domanda da porsi è se la soluzione numerica possa dipendere dalle con¬ 
dizioni iniziali, e se non vi sia contraddizione col modello teorico pre¬ 
scelto. Questo implica analizzare il processo numerico d’integrazione nel 
tempo in rapporto con il modello teorico rappresentato dalle equazioni 
differenziali del moto e, in misura ancor più importante, analizzare il 
trattamento dei contorni in rapporto con il modello prescelto per il mondo 
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esterno e con l’interazione fisica di quest’ultimo con i fenomeni occor¬ 
renti nella regione di calcolo. 

Un breve commento, sul primo punto. La meccanica dei fluidi com¬ 
pressibili, procedendo dalle idee di Riemann, sa che un’evoluzione fisi¬ 
ca è il risultato dell’interazione di segnali provenienti da regioni diverse 
e convergenti in un punto. Questo concetto dev’essere rispettato da ogni 
codice di calcolo; in altre parole, il codice deve, per quanto possibile, 
aggiornare i valori fisici in un punto usando soltanto informazioni pro¬ 
venienti dal dominio di dipendenza del punto stesso. Uno dei pregiudizi 
più diffusi nel primo ventennio di storia della meccanica dei fluidi nu¬ 
merica permetteva di approssimare le derivate spaziali con differenze finite 
prese all’avanti, all’indietro, centrate, o comunque definite, purché l’e¬ 
spressione tendesse alla derivata spaziale al decrescere dell’intervallo. Solo 
nell’ultimo decennio l’importanza del dominio di dipendenza e della coe¬ 
renza fisica, e non solo formale, tra il modello e la sua interpretazione 
numerica, è stata riconosciuta. Con gran piacere rammento che alcuni 
fra i primi segnali furono lanciati da questo Politecnico nel 1976, insie¬ 
me con alcuni dei più efficaci metodi numerici scevri dai pregiudizi di 
cui sopra. 

Quanto al secondo punto, non si dimentichi che le equazioni diffe¬ 
renziali del moto (a buon diritto chiamate «indefinite» al tempo di Levi- 
Civita) descrivono qualunque problema nell’ambito di un certo model¬ 
lo, ma che un problema specifico è definito dalle sue condizioni inziali 
e di contorno. 

Ci sono contorni rigidi e, sovente, contorni tracciati arbitrariamente 
in un campo di moto, per necessità pratiche di definizione di una regio¬ 
ne finita di calcolo (contorni permeabili). Il modello dei contorni rigidi, 
per moti non viscosi, prescrive l’annullarsi della componente normale 
della velocità. Incredibilmente, per molti anni i codici di calcolo aggiun¬ 
sero a tale condizione altre condizioni affatto arbitrarie, come un certo 
«principio di riflessione» applicato a tutte le variabili, senza capire che 
l’unica grandezza fisica che si «riflette» è il segnale procedente lungo 
la caratteristica ben definita e che la «riflessione» è la reazione della pa¬ 
rete al segnale incidente. Il panorama è migliorato notevolmente negli 
ultimi anni. 


Non altrettanto si può dire per i contorni permeabili, dato che siano 
subsonici. Qui è necessario avere un modello del mondo esterno, e un 
metodo per rappresentarlo. Tanto per fare un esempio, si consideri un 
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corpo investito da una corrente subsonica, uniforme all’infinito. Del mon¬ 
do esterno si conoscono con precisione solo due cose: la temperatura 
totale, che in condizioni stazionarie si conserva inalterata in tutto il cam¬ 
po, e l’entropia che, in assenza di urti, pure è costante in tutto il campo. 
Checché avvenga all’interno della regione calcolata, tali due quantità fi¬ 
niranno per penetrarvi, attraverso l’arco del contorno «di entrata», quan¬ 
do le perturbazioni generatesi a distanza finita abbiano raggiunto una 
distanza sufficientemente grande. È quindi logico interpretare il model¬ 
lo del mondo esterno interagendo i segnali provenienti dall’interno con 
opportuni segnali provenienti dall’esterno, che trasportino la tempera¬ 
tura totale e l’entropia al contorno permeabile. Purtroppo, in problemi 
bidimensionali, un altro dato va prescritto sul contorno di entrata. Qui 
un’aggiunta al modello del mondo esterno si impone. La quantità più 
innocua da prescrivere sembra sia la pendenza delle linee di corrente. 
Guardando il problema con la mente di un ingegnere, si può dire che 
ciò equivale a disporre sul contorno una serie di guide così fatte da de¬ 
viare la corrente di quanto si vuole senza infirmare l’energia totale. Ov¬ 
viamente, si cercherà di scegliere una distribuzione di pendenze 
ragionevole, per esempio quella ottenuta da una soluzione incompressi¬ 
bile a potenziale. 

È interessante osservare che uno studio dettagliato del problema mo¬ 
stra che ogni altro modello è un pregiudizio, e che la soluzione conse¬ 
guente non è «la» soluzione! 

Quanto sopra riguarda il modello esterno per ciò che si riferisce alla 
scelta delle variabili significative. Non migliore è la situazione nei con¬ 
fronti dei metodi proposti (e, purtroppo, sovente accettati) per calcola¬ 
re punti di contorno (indifferentemente, d’entrata o d’uscita) 
minimizzando le informazioni dal mondo esterno. Qui regna la confu¬ 
sione fra concetti stazionari e non. È ovvio che, se il contorno è stato 
arbitrariamente tracciato attraverso il campo di moto, il mondo ester¬ 
no, vicino al contorno, deve ancora comportarsi come il mondo inter¬ 
no. Se una perturbazione è prodotta nell’interno del campo calcolato, 
essa tende a propagarsi verso l’esterno. Dovrebbe essere altrettanto ov¬ 
vio che, arrivando dal contorno, la perturbazione interagisca con ciò che 
già esiste (e generalmente varia esso stesso) nel mondo esterno, e quindi 
una nuova perturbazione si muova verso l’interno, e così via. Il pregiu¬ 
dizio alla base delle cosiddette «condizioni al contorno non riflessive» 
consiste nell’ignorare le onde procedenti dal contorno verso l’interno, 
creando quindi la condizione paradossale di un flusso che è autodeter¬ 
minante, qualunque sia il mondo che lo circonda. 
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Essendo il male assai diffuso, mi sono domandato quali ne fossero 
le origini. Credo di aver trovato una spiegazione, che si può dividere 
in due parti. La prima mi è stata suggerita dal fatto che sovente queste 
condizioni sono attribuite a Sommerfeld il quale, purtroppo, non è più 
in grado di rifiutarne la paternità. Lavorando in forma analitica sull’e¬ 
quazione di D’Alembert in spazi infiniti (e questo è essenziale), si pos¬ 
sono individuare onde che procedono da punti al finito verso l’infinito 
e viceversa. L’intensità di queste ultime dev’essere nulla all’infinito, in 
assenza di singolarità. Pertanto, o non esistono affatto o sono affette 
da coefficienti che tendono a zero all’infinito; considerazioni di questo 
tipo, insieme con analoghe condizioni a distanza finita, permettono di 
definire la soluzione del problema dato. Salvo il caso molto particolare 
in cui non esistano onde provenienti dall’infinito, tali onde sono pre¬ 
senti nella soluzione a distanza finita. 

Inoltre, alcuni matematici amano dimostrare certe proprietà sotto con¬ 
dizioni restrittive, generalmente non verificate in casi reali, salvo uno, 
che essi scelgono ad esempio di controllo del loro teorema. Si pensi ad 
esempio ad un condotto trattato in forma quasi-monodimensionale; sia 
il condotto ad area costante nelle zone contenente i contorni permeabili 
e al di fuori di esso fino a distanze infinitamente grandi. Ogni perturba¬ 
zione generata nell’interno del condotto si propaga verso l’esterno e fuo¬ 
riesce dalla zona di calcolo come un’onda semplice: dall’esterno non arriva 
nessuna perturbazione. Ovviamente, le condizioni al contorno sono «non 
riflessive» e la soluzione è in ogni istante completamente definita dalle 
condizioni iniziali e dalle eventuali perturbazioni successive create nel¬ 
l’interno della regione di calcolo. Ogni altro problema, per esempio l’a¬ 
nalogo problema in due dimensioni, deve avere differenti condizioni al 
contorno. Onde circolari provenienti dall’interno di un contorno circo¬ 
lare non sono onde semplici ed esprimono propagazioni verso l’esterno 
come pure verso l’interno. Di conseguenza, condizioni al contorno «non 
riflessive» che sopprimono le onde centripete travisano il problema. 

3) Molto si è discusso sui meriti ed inconvenienti di due metodi di 
calcolo per le equazioni di Eulero, l’uno basato sulle equazioni di con¬ 
servazione di massa, quantità di moto ed energia e l’altro basato sulla 
teoria delle caratteristiche ed il calcolo esplicito degli urti («shock-fitting»). 
È interessante notare il pregiudizio dei «conservatori», che un calcolo 
numerico del secondo tipo non conservi quanto dev’essere conservato. 
Questo non è vero, se il calcolo è fatto bene; d’altra parte, i «conserva- 
tori» dimenticano spesso che i loro codici di calcolo non conservano l’en¬ 
tropia, ovunque essa va conservata! 
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4) Un esempio di pregiudizi sorti dal voler oltrepassare i limiti di un 
modello venne alla luce qualche anno fa attraverso i tanti calcoli intesi 
a trovare metodi efficienti per il calcolo di moto transonici intorno a 
profili alari. Usando il modello di flusso compressibile a potenziale, ri¬ 
sultò che per un dato profilo la portanza, negativa per incidenze negati¬ 
ve, restava negativa aH’aumentare dell’incidenza fino a circa 2/10 di grado 
e poi scattava di colpo a valori positivi. 11 fenomeno opposto si verifica¬ 
va partendo da valori positivi e diminuendo l’incidenza. Per un breve 
periodo di tempo, si parlò di «isteresi» nel rapporto incidenza-portanza. 
Questo è ovviamente un pregiudizio, dato che un fenomeno assai com¬ 
plesso come l’isteresi non può essere collegato a un modello così sempli¬ 
ce come il flusso a potenziale. Ben presto, infatti, Salas scoprì che la 
curva portanza-incidenza poteva avere tre soluzioni nell’intorno dell’in¬ 
cidenza nulla, e che tali soluzioni potevano essere determinate senza am¬ 
biguità pur di prescrivere la portanza e calcolare l’incidenza, anziché 
procedendo alla rovescia come fatto dianzi. Che le anomalie riscontrate 
nel modello a potenziale siano prive di senso fisico fu poi mostrato stu¬ 
diando lo stesso problema nel modello euleriano, il quale provvede una 
soluzione unica e continua della portanza in funzione dell’incidenza. 

5) In calcoli di problemi stazionari, forzare uno o più valori noti pe 
la soluzione stazionaria a rimanere costanti nel corso del calcolo, men¬ 
tre le equazioni non stazionare li trattano come variabili, è un errore 
abbastanza comune. Il pregiudizio sta nel credere che ciò acceleri la con¬ 
vergenza verso la soluzione asintotica, quando invece si può creare un 
danno irreparabile. Infatti, la forzatura contraddice il modello non sta¬ 
zionario. 

6) Interessanti fenomeni si scoprono attraverso il calcolo numerico, 
ma la loro interpretazione in termini del modello prescelto può essere 
affetta da pregiudizi. Si consideri, per esempio, il moto, uniforme al¬ 
l’infinito, intorno a un cerchio. Siano x ed y due assi ortogonali, pas¬ 
santi per il centro del cerchio, con l’asse x parallelo alla velocità 
all’infinito. A bassi numeri di Mach il moto non può che essere irrota¬ 
zionale e quindi simmetrico sia rispetto all’asse x che all’asse y. Ogni 
scostamento da queste simmetrie o, peggio, la separazione del flusso a 
valle e la formazione di scie vorticose non è che la rivelazione dell’ina¬ 
deguatezza della tecnica numerica. All’aumentare del numero di Mach 
all’infinito al di sopra di 0,405 circa, il moto diventa transonico; due 
urti si formano simmetricamente rispetto all’asse x. La rotazionalità della 
corrente dietro gli urti può creare due zone di separazione e ricircolazione. 
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La quantità di vorticità in esse intrappolata dipende dalla storia dell’e¬ 
voluzione. Inoltre vi sono due soluzioni possibili: una stazionaria e sim¬ 
metrica rispetto all’asse x, ma instabile, ed una stazionaria in cui la 
separazione (non viscosa) avviene in modo alternato e periodico nei se¬ 
mipiani superiore ed inferiore, con formazione di vortici che produco¬ 
no una scia di Karman con un numero di Strouhal appropriato. Il moto 
degli urti è stato elegantemente spiegato da Pandolfi in termini fisici e 
l’insieme dei fenomeni è un esempio di «ciclo limite» con chiare e sem¬ 
plici interpretazioni fisiche. 

7) In altri casi, tuttavia, l’uso dell’espressione «ciclo limite» è un pre¬ 
giudizio, non avendo alcuna interpretazione fisica nè nel modello pre¬ 
scelto nè in qualunque altro modello. Questo avviene, per esempio, 
quando si calcola un flusso contenente urti, siano essi espressi esplicita¬ 
mente come discontinuità ovvero, come suol dirsi, «catturati». Nell’av- 
vicinarsi a stati stazionari asintoticamente, gli urti possono vibrare in 
maniera impercettibile sullo schermo di un computer, ma tale da far oscil¬ 
lare il residuo fra valori di diverso ordine di grandezza, in modo perfet¬ 
tamente periodico e con una frequenza di pochi cicli iterativi. Il fenomeno 
è puramente numerico; nel caso di «shock-fitting» le sue cause e il rime¬ 
dio sono facilmente individuati: piccole variazioni nella posizione del¬ 
l’urto lo portano a dipendere da altri punti della griglia di base, e quindi 
a spostarsi in direzione opposta, ripetendo il ciclo. Nel caso di urti cat¬ 
turati, credo che la causa sia analoga. Il limitarsi a battezzare il fenome¬ 
no con un nome prestato non è buona scienza. 

Si può, da questa breve analisi, limitata a un campo assai ristretto 
di una sola delle tante branche della scienza, ricavare qualche conclu¬ 
sione? A costo di essere tacciato di pregiudizio, mi permetto di pensare 
che gli esempi che ho portato non siano fatti isolati, da curarsi nell’am¬ 
bito di una certa tecnologia, ma siano invece tipici esempi di un male 
diffuso. Se ciò è vero, la conclusione deve uscire dal campo ristretto della 
meccanica dei fluidi, press’a poco in questa forma: 

I pregiudizi sono figli della superficialità, che è sempre associata con 
l’arroganza e la presunzione. La nostra epoca ne è particolarmente ric¬ 
ca, e tocca alla scuola rimettere in auge modestia, autocritica e disciplina. 
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Introduzione 

La realtà fisica dei flussi in regime ipersonico è molto complessa. Qui 
di seguito ci si propone di esaminarne un aspetto: l’interazione fra flui¬ 
dodinamica e fenomenologia di non-equilibrio. 

Per poter mettere meglio in risalto tale interazione, si assume una mo- 
dellizzazione alquanto cruda per la fluidodinamica. Il fuido è conside¬ 
rato compressibile e non stazionario. Inoltre, trascurando i fenomeni 
di trasporto, esso risulta pure non viscoso e non conducibile termicamente. 

Dei tre fenomeni tipici di non-equilibrio, quello termodinamico, re¬ 
lativo alla energia di vibrazione, quello chimico e quello della ionizza¬ 
zione, si considera solo il secondo. Per l’energia vibrazionale, si può 
ipotizzare il gas eccitato parzialmente o totalmente, oppure in equilibrio. 
Qui si è assunto il livello di eccitazione corrispondente ad un grado di 
libertà. In ogni caso, tale assunzione non è importante nel contesto che 
segue. La ionizzazione è del tutto trascurata. Tali due ipotesi concorro¬ 
no alla semplificazione del modello e sono, peraltro, pienamente ragio¬ 
nevoli per un ampio campo di condizioni tipiche del rientro. Un modello 
di cinetica chimica descrive il relativo fenomeno di non-equilibrio. 

In conclusione, la fuidodinamica risulta descritta dalle equazioni di 
Eulero ed il processo chimico di non-equilibrio da una cinetica, succes¬ 
sivamente illustrata. 

Qui di seguito vengono analizzati alcuni punti caratteristici della in¬ 
terazione fra fluidodinamica e chimica in non-equilibrio. In particola¬ 
re, vengono considerate le equazioni che descrivono il fenomeno e la 
forma a cui esse tendono nelle condizioni limite, rappresentate rispetti¬ 
vamente dal flusso congelato e da quello in equilibrio. Viene quindi ana¬ 
lizzata la propagazione ondosa nel fluido reagente, in non-equilibrio e 
nella configurazione di pieno equilibrio. Infine si ricorda come l’irre- 
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versibilità propria della evoluzione in non-equilibrio comporti la forma¬ 
zione di vorticità nel seno della corrente. 

Le considerazione che seguono, sviluppate sulla base dei suggerimenti 
forniti in bibliografia classica (per esempio [1], [2]), appaiono interes¬ 
santi sia per l’indagine del fenomeno fisico che, in misura non trascura¬ 
bile, per le indicazioni fornite nel condurre una analisi critica dei risultati 
forniti dalla indagine numerica. 


Le equazioni 

Si considera un fluido compressibile, non viscoso, non conducibile 
termicamente, in regime non stazionario e miscela di gas perfetti. Con 
le notazioni usuali, le risultati equazioni di Eulero (continuità, quantità 
di moto ed energia), scritte in forma quasi lineare, risultano: 


(1) 

( 2 ) 

(3) 


d e _ 

Dt 


= 0 


Dv 

DT 

Dh 

Dt 


— V p = 0 
Q 


1 Dp 
q Dt 


La forma quasi lineare è conveniente poiché consente di mettere in 
risalto gli aspetti di propagazione ondosa che non sono direttamente evi¬ 
denziati nelle leggi di conservazione, cui, peraltro, si fa ricorso nella ap¬ 
plicazione numerica qualora si intenda catturare numericamente 
discontinuità fluidodinamiche. 

Nel modello di cinetica chimica qui assunto, si considera la miscela 
di cinque specie ( O, N, NO, 0 2 , N 2 , indicate rispettivamente con i= 1, 
2, 3, 4, 5). 

Si indica con Y, la concentrazione in massa della specie /-esima, che 
concorre alla formazione della miscela, e con g, la sua massa mole¬ 
colare. 

L’equazione di stato per la miscela è data da: 


(4) 




Y, 
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L’entalpia che compare nella Eq. 3 è legata alla temperatura ed alla 
energia complessiva di formazione ( h for ): 

(5) h = h M + £ c pì y ìT 


hfor= E w 


5 

~2 



Hi 


2 ); 



— (i = 3, 4, 5) 
M; 


La cinetica chimica è descritta dalle equazioni di produzione delle spe¬ 
cie. Tali equazioni risultano: 


In esse, la velocità di produzione delle specie /-esima bilancia la varia¬ 
zione locale nel tempo e la convezione della concentrazione, essendosi 
trascurata la diffusione. 

Tali velocità di produzione delle specie (co, espresse in g-mole-g~ { - 
- 5- 1 ) dipendono da proprietà termodinamiche locali (p. q) e dalle con¬ 
centrazioni (T,): 

(7) W, = u>,(p, e, YJ) , Vi 

In virtù della conservazione delle due specie atomiche ( O, N), due delle 
Eqq. 6 sono superflue e di queste se ne possono considerare solamente 
le prime tre. Per una valutazione quantitativa delle velocità co, si rinvia 
il lettore a [3], ove sono riportate formule e costanti necessarie. 

L’interazione fra fluidodinamica e chimica è essenzialmente conte¬ 
nuta nelle Eqq. 5 e 7. Nella prima le concentrazioni Y„ derivate dalla 
cinetica chimica (Eqq. 6), definiscono il legame tra entalpia, pressione 
e densità, soprattutto attraverso l’effetto dei calori di formazione ( h °). 
Nella seconda, la pressione e la densità (in particolare il loro rapporto, 
proporzionale alla temperatura) influenzano fortemente le velocità di 
reazione. 
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È utile rammentare due modelli alternativi di miscele reagenti, che 
si possono considerare come casi limite del flusso in non-equilibrio. Il 
primo è il modello di flusso congelato, per il quale alle equazioni di pro¬ 
duzione delle specie (Eqq. 6) si sostituisce la condizione che le concen¬ 
trazioni permangano invariate sulle linee di flusso e, per ingressi uniformi, 
in tutto il campo. In tal caso la Eq. 5 risulta semplificata per la costanza 
di h for . Il secondo è rappresentato dal modello del flusso in equilibrio, 
le cui condizioni sono qui di seguito indicate con asterisco (*). Le con¬ 
centrazioni di equilibrio sono fornite da equazioni algebriche: 

(8) Yj= Yf (p, q) 

Anche in tal caso, nella Eq. 5 l’entalpia dipende solo dap e q, esplicita¬ 
mente dalla temperatura ed indirettamente attraverso le concentrazioni 
di equilibrio che concorrono alla determinazione di h for . 

È opportuno ricordare qui la definizione della velocità del suono, co¬ 
me la velocità di propagazione di un disturbo a gradino che mantiene 
inalterata la sua forma [1], Per il flusso congelato, essa è data da: 

, h. 

(9) aj = -- e — 



e 


mentre per il flusso in equilibrio si ha: 


(10) 


h e + L,h Y; Yl 
— -h p -y.,h Y : y] p 

Q 


Nella prima definizione (Eq. 9), le concentrazioni non appaiono, poi¬ 
ché permangono inalterate attraverso il disturbo ( dp, dq). Nella secon¬ 
da definizione (Eq. 10) appaiono le concentrazioni di equilibrio, attraverso 
le loro derivate Y* e o Y* ip . Tali derivate esprimono le variazioni delle 
concentrazioni (dYf) che si accompagnano alle variazioni dp e dq che 
caratterizzano il disturbo. 

È conveniente eliminare la derivata sostanziale della densità dalla equa¬ 
zione di continuità (Eq. 1). Attraverso la definizione della entalpia (Eq. 
5) e l’equazione dell’energia (Eq. 3), si ottiene: 
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Dq _ 1 Dp 
Dt aj Dt 


dove: 

OD £ h y. CO, /X/ 

L’equazione di continuità (Eq. 1) diventa dunque: 


( 12 ) ^. + ea j^-d=aji 

Si mette in evidenza che nella Eq. 12, equazione di continuità per flusso 
in non-equilibrio, compaiono due elementi importanti. Da un lato, la 
velocità del suono, definita come quella del flusso congelato (Eq. 9). Dal¬ 
l’altro, il fattore \p, che contiene gli elementi più significativi della chi¬ 
mica, fra cui le costanti degli equilibri chimici, proprio quelle che 
determinano le relazioni delle Eqq. 8. 

Infine è importante ricordare che l’equazione di continuità per il flusso 
congelato risulta essere: 

Dp 7 

(13) — + ea;v-v=o 


molto simile alla Eq. 12, dalla quale differisce solo per l’assenza del ter¬ 
mine di sorgente legato a \p. Al contrario, per il flusso in completo equi¬ 
librio, si ottiene: 


(14) — + G a^V-v=0 


Qui la differenza con il non-equilibrio (Eq. 12) è duplice: l’assenza del 
termine di sorgente (come già per la Eq. 13) e l’introduzione della velo¬ 
cità del suono in condizioni di equilibrio (Eq. 10) in luogo di quella de¬ 
finita per il flusso congelato. 
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Flusso in non-equilibrio, congelato ed in equilibrio 

Un parametro che può definire il grado di non-equilibrio in un flusso 
di miscela reagente è il numero di Damkòhler (Da). Esso è proporziona¬ 
le al rapporto di una lunghezza l x , tipica della fluidodinamica, su di 
una lunghezza l ch , lungo la quale si sviluppa il rilassamento chimico. 
Parimenti, Da è proporzionale al rapporto dei relativi tempi t x e t ch . È 
interessante osservare la natura del flusso in non-equilibrio a differenti 
valori di Da. In tali analisi è conveniente assumere valori di riferimento 
prefissati per le proprietà termodinamiche (p■ ej e per la velocità 
(uj. Per problemi di rientro, ciò significa assumere il numero di Mach 
diavolo e l’altitudine. La variazione di Da risulta determinata da una 
lunghezza tipica del corpo ( IJ , ossia dalla sua scala fisica. 

Si ipotizzino valori molto piccoli di l x (per esempio veicoli di rien¬ 
tro di dimensione ridottissima), che portano a Da~>0. I tempi di resi¬ 
denza del gas aH’intorno del corpo sono così ridotti, rispetto ai tempi 
di rilassamento chimico, che non si formano apprezzabili concentrazio¬ 
ni di nuove specie chimiche. Il flusso congelato appare dunque il limite 
del flusso in non-equilibrio, per Da -* 0. Si noti che le velocità delle rea¬ 
zioni chimiche non sono nulle, anzi possono essere molto elevate. 

Al contrario, si ipotizzino valori molto alti di l x , indicativi di veico¬ 
li molto grandi. Ora Da~*°°. I suddetti tempi di residenza del gas di¬ 
ventano così alti, che c’è modo di sviluppare nella miscela tutte le reazioni 
chimiche sino alle condizioni di equilibrio in accordo con il campo flui¬ 
dodinamico generato dalla presenza del corpo. Il flusso in equilibrio rap¬ 
presenta ora il limite del flusso in non-equilibrio per Da -* oo. Si osservi 
che, in questo caso, le velocità di reazione tendono ad annullarsi. 

È interessante osservare come le equazioni del flusso in non-equilibrio 
si comportino in questi due casi limite. Si prenda in considerazione l’e¬ 
quazione più semplice, quella della produzione delle specie (Eqq. 6). La 
si analizzi nella scrittura originaria, scritta nel piano fisico, senza ope¬ 
rare alcuna adimensionalizzazione, e poi nella forma più consona al cal¬ 
colo numerico (che non può prevedere grandezze tendenti a zero o 
infinito), nella quale le grandezze, lunghezze e tempi in particolare, ap¬ 
paiono normalizzati e scalati con riferimento a l m . La Eq. 6, nella for¬ 
ma normalizzata, presenta nel termine di sorgente una velocità di reazione 
0,-, definita come segue: 
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Per l x -► 0, il termine di sorgente co, tende ad un valore finito, e così la 
derivata sostanziale della concentrazione nel piano fisico. Ciò non è con¬ 
traddittorio, poiché tale valore finito è il limite di variazioni delle con¬ 
centrazioni tendenti a zero su spazi e/o tempi che tendono anch’essi a 
zero. Ossia, non si hanno apprezzabili risultati delle reazioni chimiche, 
non perché le velocità di reazione siano nulle, ma perché spazi e tempi 
a disposizione sono evanescenti. Dunque, la Eq. 6, scritta nel piano fisi¬ 
co, ben riflette l’interpretazione del flusso in non-equilibrio che tende 
alla configurazione congelata per Da^O. Ovviamente alle stesse con¬ 
clusioni conduce la stessa equazione scritta in forma normalizzata, quella 
che fornisce indicazioni sulla soluzione numerica. Tuttavia, si nota che 
qui il termine di sorgente 0, -^0 per l x -*0. Questo è in accordo con dif¬ 
ferenze di concentrazione nel campo che tendono a zero, essendo, per 
contro, finite le scale degli spazi e tempi, ora normalizzati con riferi¬ 
mento a / . Dunque il processo di limite al flusso congelato non pre¬ 
senta problemi, sia per l’interpretazione fisica che per l’analisi numerica. 

Per / ->oo, il flusso in non-equilibrio appare tendere al flusso in 
equilibrio. Considerando la Eq. 6 sul piano fisico, poiché la velocità di 
reazione tende a zero, si ottiene che pure la derivata sostanziale ten¬ 
de a zero. Tuttavia ciò non avviene per la costanza delle concentrazioni; 
anzi, la configurazione di equilibrio cui si tende è caratterizzata dalle 
massime variazioni di queste. Essa tende a zero poiché tali variazioni 
finite si diluiscono su spazi e/o tempi che tendono a oo. Passando ora 
alla forma normalizzata della Eq. 6, si osserva che il termine di sorgente 
fi,. può tendere ad un valore finito, essendo co, moltiplicata per l x . Ora 
la derivata sostanziale sul piano di calcolo, che bilancia fi,, presenta va¬ 
lori finiti a numeratore (variazioni finite di concentrazione) e a denomi¬ 
natore (spazi e tempi normalizzati). Questo fatto è indicativo di problemi 
che possono insorgere nella previsione numerica. Difatti, la soluzione 
che deriva dalla integrazione della Eq. 6, scritta in forma normalizzata, 
è costruita sulla base dei valori locali di fi,. Questi sono ottenuti nume¬ 
ricamente come prodotti di un valore piccolissimo (co,) per uno grandis¬ 
simo (/ ) e risultano non affidabili ad alti Da, essendo limitato il 
numero di cifre significative disponibile nella macchina calcolatrice. Pra¬ 
ticamente, per Da troppo alti, il calcolo si arresta per problemi che se¬ 
guono alla determinazione errata di fi,. Si noti che ciò avviene là dove 
la rigidezza del sistema fluidodinamica-chimica diventa critica e che il 
problema della valutazione incerta di fi, ne è una diretta conseguenza. 
È dunque possibile concepire il flusso congelato ed il flusso in equili¬ 
brio come configurazione limite del non-equilibrio, al variare di Da. 
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Tuttavia, mentre il limite al flusso congelato si presenta chiaro, non al¬ 
trettanto si può dire del limite al flusso in equilibrio. Si osservino le Eqq. 
12 e 13 scritte in forma normalizzata (in particolare rispetto alla lun¬ 
ghezza l x ). Il termine di sorgente della Eq. 12, risulta ora definito come: 



« c 


Le due equazioni tendono chiaramente a coincidere per Da^> 0, poi¬ 
ché ^ -> 0 per I x -> 0. 

Al contrario, per l x -* oo, con tendente ad un valore finito, la Eq. 
12 non appare esplicitamente tendere alla Eq. 14 per la presenza di due 
fattori: la velocità del suono, che moltiplica la divergenza della veloci¬ 
tà, è quella congelata nella Eq. 12 e quella definita di equilibrio nella 
Eq. 14. Inoltre il termine di sorgente 4', finito nella Eq. 12 per Da -*• oo, 
non compare nella Eq. 14. Concludendo, si può dire che resistenza di 
almeno due fattori che differiscono nelle equazioni non esclude di poter 
sperare che, con indagine più approfondita, si possa riconoscere che i 
due si compensino e che le equazioni siano equivalenti. 


Propagazione ondosa in flussi reagenti 

Il problema emerso nel confronto fra le equazioni di continuità del 
flusso in non-equilibrio, a Da^> oo, e di quello in equilibrio risulta an¬ 
cor più evidente dalla scrittura, sempre in forma normalizzata, delle equa¬ 
zioni di compatibilità che descrivono la propagazione di segnali lungo 
linee caratteristiche. 

Dalla forma quasi lineare delle equazioni per il flusso in non-equilibrio 
si ottiene: 

dR X3 , dR x 3 2 T 

(15) -—-+ (U T a f ) —— = a} ¥ 

dt dx 


ove: 

(16) dR U3 = dpT Qa f du . 


Per contro, per il flusso in equilibrio, risulta: 
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(17) 

ove: 


dRl 3 

dt 


+ (u -F 


a eq ) 


ÒRI 3 

dx 


= 0 


(18) dR*ì = dpTQa ei} du. 

La differenza fra le Eqq. 15 e 17, allorché Da^> o°, e quindi ’L tende 
ad un valore finito, sono evidenti. Le linee caratteristiche su cui propa¬ 
gano i segnali differiscono per la diversa velocità del suono. Per lo stes¬ 
so motivo, i segnali stessi (Eqq. 16 e 18) risultano definiti in modo diverso. 
Infine, per il flusso in equilibrio si può parlare di invarianti di Riemann 
(almeno per flusso omoentropico), mentre nel caso del non-equilibrio 
il termine di sorgente non lo consente. 

L’equivalenza della Eq. 15, a Da~> oo, e della Eq. 17 può essere 
intuita osservando che le velocità di reazione, che compaiono nel termi¬ 
ne di sorgente >K, contengono nella loro definizione le costanti di equili¬ 
brio che determinano le relative concentrazioni (Eq. 8). Dunque, i termini 
a f e nella Eq. 15 potrebbero rendere questa espressione equivalente 
alla Eq. 17, ove compare la sola a eq . L’intuizione viene confermata dal¬ 
la trattazione condotta nell’ambito della approssimazione acustica e, per 
semplicità, di una sola reazione. Si indica qui di seguito con l’apice (') 
la perturbazione rispetto a valori costanti di base. In questa approssi¬ 
mazione la velocità di produzione della specie risulta: 


d Y' 

~di 


Y* ' - T 


Tale velocità è proporzionale direttamente allo scostamento della con¬ 
centrazione rispetto al valore di equilibrio ed inversamente ad un tempo 
di rilassamento tipico della reazione (r). È conveniente introdurre l’o¬ 
peratore che indica la derivazione lungo le linee caratteristiche nel flus¬ 
so congelato: 

D f d d 

— — — t a, — 

Dt dt 1 dx 


Definito il segnale in non-equilibrio: 
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(19) dR' {i -dp' =F Qa f du' 

l’equazione di compatibilità, corrispondente alla Eq. 15 risulta essere: 


( 20 ) 


D f R[ 3 2 Y*'-Y’ 

-— a a f - 

Dt 1 7 


Derivando ambo i membri della Eq. 20 rispetto al tempo e ricordando 
le equazioni del moto, si riconosce come dal termine di sorgente viene 
ad emergere la descrizione della propagazione di segnali nel flusso in equi- 


librio: 


(21) 

dt 1 

ove: 

D eq 


~dT 


DfRi 3 


_d_ . 
dt 


D eg R*J 


d 

dx 


ed il segnale in non-equilibrio risulta definito da: 

(22) dR* \ = dp' =f QCt eq du’ 

Si nota come la Eq. 21 esprime la propagazione ondosa in un mezzo 
reagente in non-equilibrio e come il fattore di scala l x modula l’influen¬ 
za del contributo del flusso congelato e di quello in equilibrio. Per valo¬ 
ri piccoli di /a,, ossia per Da 0, la propagazione ondosa è descritta in 
termini del flusso congelato: 


(23) 


D f Rj ì 

Dt 


-0 


Al contrario, per valori elevati di l x è il contributo del flusso in equili¬ 
brio ad emergere 


D eg Rl 3 


= 0 


(24) 


Dt 
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Dissipazione e non-equilibrio 

La produzione di entropia nell’ambito delle equazioni di Eulero può 
aver luogo solo attraverso onde d’urto. Se queste ultime presentano in¬ 
tensità variabile nel campo di moto, i livelli di entropia, a valle delle stesse, 
risultano differenti su linee di flusso limitrofe. È noto che la generazio¬ 
ne di tale gradiente di entropia comporta la nascita di vorticità nella 
corrente. 

Anche i processi di cinetica chimica in non-equilibrio portano alla pro¬ 
duzione di entropia. In particolare, la variazione di entropia nella mi¬ 
scela reagente è data da: 

Tds = dh- — dp+ V frdYt 

Q T 

ove f,(A G- Y ì) è legato ai potenziali chimici delle specie e alle concen¬ 
trazioni. 

Considerando le variazioni incontrate da un elemento di massa di mi¬ 
scela nella sua traettoria e rammentando l’equazione dell’energia (Eq. 3), 
si ottiene, nella scrittura normalizzata: 

Ds yi 5. D Y i 
Dt = T — ' Dt 

In condizioni di non-equilibrio, la produzione di entropia lungo il tra¬ 
gitto di un elemento di miscela è tanto maggiore quanto più forte è il 
disequilibrio locale e più intensa la generazione di specie. Tuttavia, tale 
produzione di entropia scende a zero nei due casi limite di Da -+ 0 e 

Nel primo caso è evanescente la generazione di specie (DY/Dt-+ 0), 
pur essendo finito Nel secondo, la derivata sostanziale della concen¬ 
trazione DY/Dt è finita, come già visto precedentemente. Ora però il 
termine f,- tende a zero, poiché C* (p, q, Y*) = 0. Dunque, per la dissi¬ 
pazione, il limite del non-equilibrio appare regolare sia verso la confi¬ 
gurazione congelata che verso quella in equilibrio, essendo per entrambe 
nulla la produzione di entropia. 

La generazione di vorticità collegata ai differenti livelli di entropia 
su linee di corrente limitrofe è una ulteriore interessante manifestazione 
della interazione fra fluidodinamica e non-equilibrio. 
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Un esempio semplice ma significativo 

Fra i numerosi problemi fisici in cui viene messa in risalto l’interazio¬ 
ne fra fluidodinamica e non-equilibrio chimico, se ne ritrova uno sem¬ 
plice per la sua definizione e molto esauriente per la fisica in esso descritta. 

Si consideri una corrente ad alta velocità, decisamente supersonica, 
ed ad alto livello entalpico, in condizioni di equilibrio, quindi con rile¬ 
vante dissociazione. Tale corrente lambisca una parete che presenta una 
brusca deviazione attraverso un angolo di parete convesso. Se ci si limi¬ 
ta all’ambito delle sole equazioni di Eulero, gas perfetto nelle condizio¬ 
ni di monte, il fenomeno è semplice e chiaro. Esso è noto come 
l’espansione di Prandtl-Meyer. Il campo di moto è formato da due re¬ 
gioni uniformi a monte e a valle, separate da un ventaglio di espansione 
delimitato da due caratteristiche. La prima, confinante con il campo di 
monte, è determinata dalle sole condizioni del flusso a monte. La se¬ 
conda delimita il ventaglio sì che la corrente in uscita dallo stesso sia 
parallela alla parete deflessa. Il campo di moto non è caratterizzato da 
alcuna lunghezza di riferimento e appare seìf-similar. 

Tuttavia il processo chimico introduce una lunghezza di riferimento, 
collegata ad un tempo di rilassamento chimico ed annulla il carattere 
geometrico di self-similarity. 

In particolare, le particelle che corrono molto vicine alla parete subi¬ 
scono una espansione del tipo alla Prandtl-Meyer, sulla base delle pecu¬ 
liarità fisiche del gas di monte. Infatti, l’espansione avviene in spazi così 
ridotti da ritenere congelato il flusso in prossimità dello spigolo. Si noti 
come, a valle della espansione sullo spigolo, la temperatura sia discesa 
e le concentrazioni di monte (stato fortemente dissociato) si trovino in 
forte disequilibrio chimico. Inizia dunque un processo di riassociazio¬ 
ne, dapprima molto lento e quindi sviluppato su notevoli distanze a val¬ 
le dello spigolo, lungo la parete deflessa. La lentezza della riassociazione 
è conseguenza della bassa temperatura. Allorché la riassociazione inizia 
ad aver luogo in entità non trascurabile, la temperatura cresce ed au¬ 
menta la velocità del processo. Più innanzi essa svanisce verso una con¬ 
dizione di equilibrio caratterizzato da pressione p, e temperatura T x . 

Ben diversa è l’evoluzione per particelle che corrono distanti dalla pa¬ 
rete. Qui il ventaglio di espansione è molto esteso sul piano fisico e le 
particelle subiscono una espansione diluita nello spazio e nel tempo, sì 
che le reazioni chimiche hanno modo di svilupparsi sino ad una confi¬ 
gurazione molto prossima all’equilibrio locale. Ragionevolmente si può 
dire che l’evoluzione segue una sequenza di stati di equilibrio sino al ter¬ 
mine del ventaglio, alla pressione p 2 e temperatura T 2 . Si osservi che 
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mentre a parete l’evoluzione della temperatura è caratterizzata da una 
istantanea discesa allo spigolo e da un recupero, dapprima lentissimo, poi 
più veloce, poi ancora asintotico ai valori di equilibrio {p u T x ), ora in¬ 
vece, a grande distanza sopra lo spigolo, la temperatura di una particella, 
lungo la sua traettoria, discende monotonicamente sino alla condizione di 
equilibrio (p 2 , T 2 ). Fra le due linee di flusso ora viste, scorrono parti- 
celle che seguono stati intermedi di non-equilibrio. Tale esempio mette 
bene in luce come, nello stesso problema, le equazioni del flusso rea¬ 
gente in non-equilibrio abbiano da descrivere situazioni estreme. 

Precedentemente si è visto come la velocità del suono abbia una dif¬ 
ferente definizione nel flusso congelato ed in quello in equilibrio. Di ciò 
si trova riscontro nel presente esempio. Con riferimento al problema iper¬ 
bolico del flusso stazionario bidimensionale supersonico, la prima ca¬ 
ratteristica del ventaglio di espansione presenta un angolo con la corrente 
a monte determinato dal numero di Mach, rapporto della velocità del 
gas a monte alla velocità del suono. In prossimità dello spigolo, il cam¬ 
po è dominato dal flusso congelato, con velocità del suono pari a a,. 
A forte distanza dallo spigolo, il flusso è praticamente in equilibrio e 
la velocità del suono è data da a eq . Poiché a f > a eq , tale prima caratte¬ 
ristica vicino allo spigolo è meno inclinata di quella a forte distanza dal¬ 
lo stesso. 

Infine è possibile osservare, in tale esempio, gli effetti dei diversi li¬ 
velli di dissipazione (e quindi della relativa vorticità) dovuti al non¬ 
equilibrio. Le particelle che scorrono a ridosso della parete, subito do¬ 
po l’espansione attorno allo spigolo, si ritrovano in condizioni di forte 
disequilibrio. Conseguentemente sono soggette a dissipazione. Lontano 
dalla parete, il flusso è sempre più prossimo all’equilibrio e l’evoluzione 
tende alla reversibiità. A notevole distanza a valle dello spigolo, sulle 
due linee di flusso ora considerate, si stabilisce la stessa pressione 
(P\=Pi)- Infatti il campo è isobaro poiché tutte le linee di corrente scor¬ 
rono parallele alla parete piana. Al contrario, la temperatura di equili¬ 
brio a parete (T,) è maggiore di quella a forte distanza dalla stessa 
parete ( T 2 ), causa il processo dissipativo del non-equilibrio. In corri¬ 
spondenza di tale profilo non uniforme di temperatura lungo la norma¬ 
le alla parete, si ritrova un profilo non uniforme della velocità, con valore 
minimo a parete. Da esso si riconosce subito la vorticità generata dal 
non-equilibrio. 

Questo esempio si presta molto bene come base per esperimenti nu¬ 
merici, al fine di controllare quanto le suddette particolarità della inte¬ 
razione fluidodinamica e non-equilibrio siano riprodotte dai risultati 
numerici. Ciò costituisce l’oggetto della attuale ricerca dell’autore. 
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Stability of Marangoni and combined free 
convection in shallow cavities 


Luigi G. NAPOLITANO and Antonio VIVIANI* 


Abstract. The paper deals with thè linear and energy stability analysis of 
Marangoni, naturai and combined free convection in rectangular cavities 
of small aspect ratio, open at thè top and with differentially heated end 
walls. The class of flow considered is that of parallel flows, away from 
end walls, for arbitrary values of thè Marangoni, Grashof and Prandtl num- 
bers. Smith (1982) and Hart (1983) have studied thè linear stability of Ma¬ 
rangoni convection and naturai convection, respectively. In this work we 
consider combined convection and give thè criticai Reynolds number (R c ) 
and thè Reynolds number of asymptotic stability in thè mean (R E ) for thè 
dynamic Bond number ranging from zero (Marangoni convection) to infi¬ 
nite (naturai convection), in a wide range of Prandtl numbers. The criticai 
curves of thè linear theory are given in thè piane Reynolds-Prandtl, and 
thè instability and stability domains are given in thè piane (R u , R ( J with 
R m and R 0 thè conditional Reynolds numbers based on thè Marangoni 
and buoyant speed, respectively. The limiting values, for Pr -* oo, of R c 
and R e are also calculated, and asymptotic laws for high Prandtl numbers 
are proposed, in thè whole range of definition of thè dynamic Bond number. 


1. INTRODUCTION 

Convection in thè absence of imposed velocites and/or pressure gra- 
dients in fluid systems with interfaces can be classified (Napolitano 1982) 
as naturai, Marangoni or combined free convection, according to whe- 
ther thè motion is caused by buoyant forces, Marangoni stresses (surfa- 
ce tension grandients) or both. 


* Istituto «Umberto Nobile», Università di Napoli, P.le V. Tecchio 80, 80125 Napoli, Italia 
e MARS - Microgravity Advanced Research and Support Center, Via Diocleziano 328, 80125 Na¬ 
poli, Italia. 
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The subject matter of this paper is free convection in a shallow rec- 
tangular cavity with differentially heated end walls and a free liquid 
surface. 

The bottoni boundary of thè cavity is rigid and insulated. The top 
liquid-air surface is assumed to be fiat and undeformable (this places 
constraints on thè values of thè contact angle and of thè capillary num- 
ber) and there is not dynamic and thermal coupling between thè liquid 
and thè surrounding air. The surface tension along thè free surface is 
assumed to be a linear decreasing or increasing function of thè tempera¬ 
ture. In most papers dealing with surface tension effects, thè surface ten¬ 
sion is assumed to be a linearly decreasing function of thè temperature. 
This behaviour is typical of a large class of fluids, but one may find Sy¬ 
stems like some alloys, molten salts or liquid crystals with a surface ten¬ 
sion growing linearly with temperature (Guyon and Pantaloni 1980). This 
last behaviour has been accounted for in various papers, particularly in 
Papers on Marangoni instability of liquid layers heated from below (see 
Lebon 1988). However, from thè theoretical point of view, it must be 
pointed out that, if thè surface tension is assumed to depend only on 
thè temperature, a model of surface tension linearly increasing with thè 
temperature implies that thè surface phase is thermodynamically unstable. 

The imposed horizontal temperature gradient induces a circulating con¬ 
vection flow into thè cavity. The motion is generated by density diffe- 
rences in thè gravitational field and by surface tension gradients along 
thè free surface. For small temperature differences between thè side walls, 
thè flow is laminar and steady, and for a sufficient small aspect ratio, 
and far away from thè end walls it is parallel to thè horizontal bounda- 
ries. In thè case of Marangoni convection we have a parabolic velocity 
profile, whereas in thè case of naturai convection a cubie velocity profi¬ 
le; for combined free convection we have different profiles (linear com- 
bination of thè previous ones) depending on thè value of thè dynamic 
Bond number. When thè surface tension increases with thè temperatu¬ 
re, thè resulting Marangoni convection is from cold to hot regions near 
thè free surface, in opposition to thè motion due to thè gravity forces 
(we do not consider thè case of a negative thermal expansion coefficient): 
in this case some particular profiles (e.g. zero velocity at thè free surfa¬ 
ce) can form. 

The evaluation, by means of thè classical linear stability analysis, of 
thè criticai conditions for thè onset of thè instability in terms of thè flow 
parameters, of thè nature (steady or oscillatory in time) of thè new state 
and of thè type of thè secondary flows (two-dimensional: transverse or 
longitudinal; tridimensional), as well as a stability analysis via thè ener- 
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gy method are thè object of this study. This latter method studies stabi- 
lity against disturbances of finite amplitude restricted only in certain 
smoothness requirements and in their behaviour in directions where thè 
geometry of thè flow goes to infinite (i.e. so that thè energy integrals 
exist), and it gives a sufficient condition of stability, whereas thè linear 
theory can give only a sufficient condition of instability. If R L indica- 
tes thè criticai Reynolds number found by thè linear theory and R E thè 
Reynolds number of asymptotic stability in thè mean found by thè energy 
theory, in thè region of thè parameter space where R E < Re < R L in¬ 
stability (called subcriticai instability) is not ruled out, but its existence 
must be an effect of finite amplitude. Sometimes thè subcriticai range 
is small, but in thè most cases thè energy stability boundary is conserva- 
tively low and thè subcriticai range very large. The limit cases are thè 
Rayleigh-Bénard problem for which thè linear and energy stability boun- 
daries coincide (Joseph 1965), and thè piane Couette and Hagen-Poiseuille 
flows for which thè linear theory seems to predict stability ( R L = oo) 
whereas thè energy method give, respectively, R E = 41.3 and i? £ = 81.49 
(Joseph 1966; Joseph and Carmi 1969). When thè subcriticai range is 
small, then one can have confidence that both thè linear and energy theo- 
ries have reasonably captured thè essential physics of thè onset. One of 
thè objectives of thè present study is to evaluate thè above range and 
to get insight into thè application of thè energy method for thè evalua- 
tion of sufficient condition of stability for thè class of basic flows con- 
sidered. 

Instability due to gravity forces and surface tension gradients are wi- 
dely studied in thè literature only when thè temperature gradient is ver- 
tical (thè base flow is zero), see thè monograph of Chandrasekar (1961), 
thè excellent review paper paper of Normand, Pomeau and Velarde (1977) 
and for late developments thè review paper of Lebon (1988), which ac- 
count for recent papers where a negative Marangoni number (surface 
tension increasing with thè temperature) has been considered. In this case 
of imposed vertical temperature gradient also stability analyses via thè 
energy method are available, and we mention those of Davis (1969), Le¬ 
bon and Perez-Garcia (1980), Lebon and Cloot (1982). 

Instead, when thè imposed temperature gradient is horizontal, i.e. thè 
basic flow is not zero, there is a lak of theoretical works on stability. 

The relevance of thè problem herein considered arises from tasks such 
as thè growing of semiconductor crystals (Hurle, Jakeman and John¬ 
son 1974; Gill 1974), tracking thè dispersion of pollutants in estuaries 
(Cormack, Stone and Leal 1975) or from meteorology (Hart 1983). Our 
particular interest comes from thè growth of crystals from melt by thè 
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horizontal zone-melting and horizontal normal-freeze techniques (Hurle 
and Jakeman 1981) both on earth and in microgravity environment. One 
of thè problems associated with these techniques is thè development of 
striations or segregations bands parallel to thè melt-crystal interface. The 
cause of this is recognized to be a time-dependent crystal growth speed due 
to temperature oscillations in thè melt, which in turn are generated by an 
instability of thè basic flow field in thè melt. Temperature oscillations 
have been observed experimentally by several authors in liquid of various 
Prandtl numbers, e.g., by Utech and Flaming (1967) and Stewart and 
Weinberg (1971) in molten alluminium (Pr = 0.01); by Hurle, Jakeman 
and Johnson (1974) in molten Gallium (Pr = 0.02), by Schwabe and 
Scharmann (1981) in a liquid of Sodium Nitrate (Pr = 9.3); by Kamotani, 
Ostrach and Lowry (1982) in a layer of fluorocarbon liquid (Pr = 64); 
by Legros, Limbourg-Fontaine and Petré (1986), during an experiment 
under low gravity, in an aqueous long chain alcohol solution (Pr = 7), 
which exhibit surface tension increasing with thè temperature. 

Theoretical works on thè subject are those of Smith (1982), see also 
Smith and Davis (1983), Hart (1983), Napolitano and Viviani (1989), 
who performed linear stability studies. 

Smith (1982) studies Marangoni convection (zero Bond number) and 
finds an oscillatory instability with respect to tridimensional disturban- 
ces for all thè values of thè Prandtl number. 

Hart (1983) analyzes thè case of naturai convection (infinite Bond num¬ 
ber); he finds an oscillatory instability but his analysis is limited to two- 
dimensional longitudinal disturbances and low Prandtl numbers. 

Napolitano and Viviani (1989) study thè instability of combined con¬ 
vection with respect to tridimensional non-stationary disturbances, and 
as a particular result they correct thè criticai curve reported by Hart (1983) 
and exend thè study of Hart to high Prandtl numbers. 

Here we give an energy stability analysis of thè class of flows under 
consideration for thè Bond number ranging from zero to infinity, which 
is also, at thè best of our knowledge, thè first energy stability study of 
convection in cavities with horizontal heating. To help intepreting and 
understanding thè present results, we summarize thè results of thè linear 
theory already presented by thè authors. 

Thus, after establishing thè set of relevant equations and thè base flow 
Solutions (sections 2 and 3), we summarize thè instability analysis via 
thè normal modes technique and thè relative results in sections 4 and 5. 
Then, in section 6 we study thè sufficient conditions of stability by thè 
energy method, and thè results are presented in section 7. Final com- 
ments and conclusions are given in section 8 
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2. THE MODEL 

Consider thè cavity shown in Fig. 1, of length L, height H and small 
aspect ratio H/L, with thè left wall at temperature T c and thè right wall 
at temperature T H =T C +AT. The liquid filling thè cavity is bounded 
from above by a liquid-gas interface, along which acts thè surface ten- 
sion a. The gravity vector g acts parallel to thè side walls and is direct- 
ed downward. 



Fig. 1 - Geometry of thè problem. 


The state equations for thè liquid density q and thè surface tension 
a are assumed to be, respectively: 


(1) e = e,[i-M7'-n)] 

(2) o = o r -o T (T-T r ) 

where T r is thè reference temperature, Q r and a r are density and surfa¬ 
ce tension for T= T r \ & T is thè coefficient of thermal expansion; a T is 
thè negative rate of change of thè surface tension with temperature 
(a r = - da/dT ), that we assume to be a positive or negative Constant. 
We shall not consider thè case of a negative coefficient of thermal ex¬ 
pansion, so that for a T > 0 thè Marangoni convection is superimposed 
to thè naturai convection. 

We assume that: thè liquid is Newtonian, non-elatic, non micro-polar, 
homogeneous, with Constant diffusion coefficients; thè Boussinesq 
approximation applies; viscous dissipation and radiation effects are 
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negligible. With these assumptions thè flow in thè cavity is governed by 
thè Boussinesq form of thè continuity, Navier-Stokes and energy equa- 
tions, that is non-dimensional form are: 


(3 a) 
(3b) 

(3c) 


„ ( dv 

Re — + v • V v 
\dt - 


(dT 

w 


+ Vp = V 2 y+GTi z 


= v 2 r 


where y, p, T denote non-dimensional velocity, pressure and tempera¬ 
ture, i z is thè unit vector of thè z-axis. 

The non-dimensional form results from scaling thè coordinates x, y, 
z by thè height H of thè cavity and thè velocity y(«, v, w) by therefe- 
rence velocity V r to be defined below. The non-dimensional tempera¬ 
ture is measured above thè mean temperature of thè side walls and has 
its scale ( AT) r = AT H/L. The scales for time and pressure, thè latter 
measured above thè hydrostatic pressure in thè reference state at thè tem¬ 
perature T r , are H/V r and g V r /H, with fi thè dynamic viscosity of thè 
liquid. 

The three parameters that appear in Eqs. (3) are thè «conditional» 
Reynolds (Re), Peclet (Pe) and Grashof (G) numbers (defined in terms 
of a stili undetermined reference velocity, Napolitano 1982: 

Re= V r H/v, with v thè kinematic viscosity; 

Pe = Re Pr with Pr = \i/a T thè Prandtl number in which a T is thè ther- 
mal diffusivity; 

G = {g^TH 2 )/{vV r )H/L. 

The Reynolds and Peclet numbers are «response» numbers and can 
be expressed as thè ratio of thè reference velocity to thè diffusivity velo- 
cities: Re = V r / V v and Pe = V r / V a with V v = v/H and V a = a T /H thè vi- 
scous and thermal diffusivity velocities. 

The Grashof number is a «driving» number, defined as thè ratio 
of thè gravitational velocity to thè reference velocity: G = V g / V r with 
V g = (g(3 T ATH 2 /v) H/L. 

For a systematic study of response, driving and conditional characte- 
ristic numbers, see Napolitano (1978, 1982). 

The relevant boundary conditions at thè free surface are deduced from 
thè generai balance equations for surface phases (Napolitano 1979), 
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under thè following assumptions: 

i) thè surface phase is a stream-surface devoid of any thermodynamic 
property apart from thè surface tension; 

ii) thè surface phase is in thermal and mechanical equilibrium with thè 
adjacent volume phases; 

iii) momentum and energy coupling between thè adjacent fluid phases, 
radiation effects and surface irreversibility are negligible; 

iv) thè capillary number V r /V c with V c = o r /n thè capillary velocity, is 
much less than one so that thè interface preserves its hydrostatic shape, 
that we assume to be planar (contact angle equal to 7 t/2). 

The boundary conditions at thè interface then read: 


(4a) 

w = 0 


z = 

(4b) 

du 

— +M 
dz 

11 = 0 
dx 

z = 

(4c) 

dv 

— +M 
dz 

H=o 

dy 

z = 

(4d) 

>L-o 

dz 


z = 


where thè parameter that appears in Eqs. (4b, c) is thè conditional Ma¬ 
rangoni number M = V m /V r , with V m = (o r AT/ii)H/L thè Marangoni 
velocity and it is a driving number. The boundary conditions at thè lo- 
wer rigid wall are: 

(5a,b,c) w = v = w = 0 
(5d) 

dz 

We do not consider boundary conditions at thè side walls because thè 
present analysis is devoted to thè parallel core-flow solution. 

The reference velocity V r appropriate to thè problem under study, 
according to Napolitano (1982), must be of thè order of magnitude of 
thè maximum characteristic driving velocity, i.e. 


Z = 0 
z = 0 


(6) 


V r = 0(max[\ V m \, V g ]) 
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where thè absolute value is necessary for a T can be either positive or ne¬ 
gative. We assume: 

(7 ) V r =\V m \ + V g 

In this way we satisfy thè criterion (6) in thè complete range 0 ^ Bo < oo, 
with Bo= V/l V m \ [=gAQH 2 /Ao; Aq = @ r @ r AT, A<j=\o t \AT] thè dyna- 
mic Bond number, that measures thè relative importance of thè gravita- 
tional and surface tension effects. 

Due to thè choice (7), thè Marangoni and Grashof numbers are at 
most unitary and not independent, verifying thè relation: 

(8) |M| + G=1 

By introducing thè parameter 


(9) 


V m 

~v7 


l/Bo 


Mand G can be eliminated in terms of X by means of (8) and (9), which 
give: 


(10a,b) 


1 + |\| 


G = 


1 

1 + |X| 


The choice of X as parameter of thè problem, instead of thè more useful 
Bond number Bo = 1/X, has been made in order to distinguish thè two 
cases G = 0 with M= - 1 and G = 0 with M= 1. 

The particular cases of naturai and Marangoni convection are cha- 
racterized by: 


(l la) X = 0 (M=0; G= 1) 

(llb) X= +oo (M= +1, G = 0) <x r > 0 

(Ile) X= -oo (M= -1, G = 0) a T < 0 


Marangoni and naturai convection are superimposed in thè range 
0 < X < oo and in opposition in thè range - oo < X < 0. 
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3. THE BASE FLOW 

If thè aspect ratio of thè cavity is sufficiently small, then a simple 
exact parallel-flow solution will exist away from thin turning region near 
thè end walls. The parallel flow has been considered in detail by Cor- 
mack, Stone and Leal (1975) in thè case of naturai convection, whereas 
Sen and Davis (1982) and Strani, Piva and Graziani (1983) studied Ma- 
rarangoni convection. 

The solution for combined convection is: 


(12a) 

(12b) 


u = a I - — (3 z 2 — 2z) + — (8 z 3 - 15 z 2 + 6 z) 
1 4 48 


+ — 2 —-5 — + z ì ) \ +b 
48 V 5 4 


(12c) 


— = a I - “M+ G I — 


The parallel flow arises in response to thè horizontal temperature gra- 
dient, thè lateral walls induce a pressure gradient that sustains thè re¬ 
turn flow; convection of internai energy by this end-to-end counterflow 
is then balanced by vertical diffusion. 

The solution (12) has been derived from thè application of thè singu- 
lar perturbations method, when thè field variables are expanded in an 
asymptotic series in terms of thè small parameter H/L. The constants 
(a) and (b) must be determined by matching thè solution with thè flow 
near thè lateral walls. Because only thè gradient of T appears in thè sta- 
bility problem, thè Constant (b) is inessential; thè Constant (a) is at a first 
approximation unitary, and we shall assume just a= 1: this is equiva¬ 
ler to take thè leading term of thè asymptotic expansion. With these 
approximations, thè base flow for thè stability analysis, expressed in terms 
of thè parameter X instead of M and G, is: 


(13a) 


1+[X| 
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(13b) 


(13c) 


1 + |X| l 2 ») 


We point out that (13) is an exact solution of thè problem (3), (4), (5) 
if one assumes a linear temperature distribution on thè bottom wall of 
thè cavity (Napolitano and Viviani 1990). In this case, thè higher appro- 
ximations of thè asymptotic expansion are identically zero, and thè la- 
teral walls effect is felt in thè core region only through thè generation 
of a pressure gradient in thè x-direction, that maintains zero max flux 
through any vertical piane. 

The base flow velocity profiles w(z) for X = 0 (G = 1; M=0); X = + oo 
(M= 1; G = 0); X = - oo (M= - 1; G = 0) are reported in Fig. 2a. In Fig. 
2b thè corresponding temperature profiles (7’-x)/Pe are shown. The 
turning points of thè temperature profiles coincide with thè points of 
zero velocity. 

In Fig. 3a thè velocity profiles are shown for different values of X 
and thè corresponding temperature profiles are reported in Fig. 3b; for 
X < 0 we can see two counter rotating cells. 



Fig. 2a - Velocity profiles for Marangoni convection, (a): M= 1 and X = + oo; 
(c): M= - 1 and X = - oo; and for naturai convection, (b): G= 1 and X = 0. 



((T-x)/Pe)*l 


Fig. 2b - Temperature difference profiles between free surface and bottom wall 
corresponding to thè velocity profiles of Fig. 2a. 


b c d e f g h i j 



u*10 2 


Fig. 3a - Velocity profiles for different values of thè parameter X: 

(a: X=0.25); (b: X=l/12); (c: X=0); (d: X= -1/24) (e: X= —1/12); (f: X= -0.125); 
(g:X= -0.15); (h: X= -0.1875); (i: X= -0.25); (j:\= -0.375); (k: X= -0.5) 
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h g f e d c b 



(T»)/Pe*10 3 

Fig. 3b - Temperature difference profiles between free surface and bottoni wall 
corresponding to thè velocity profiles of Fig. 3a. 

In Fig. 4a and 4b we illustrated three particular profiles, correspon¬ 
ding to thè values X = - 1/12, X = -0.15, X= -0.25 [profiles (a), (b) and 
(c), respectively]. For X= - 1/12 thè surface velocity is zero and thè cor¬ 
responding profile is symmetrical with respect to z = 1/2. For X= -0.15 
thè surface temperature is thè same òf that of thè lower rigid surface 
and both are less than thè temperature in thè liquid. For values of X of 
about -0.15 there are regions with positive and other ones with negati¬ 
ve stratification: this can strongly influence thè stability. For X= -0.25, 
thè lower celi due to thè gravity disappears; this case is characterized 
by du/dz = 0 for z = 0. 



Fig. 4a - Velocity profiles for X= - 1/12 (a); X® -0.15 (b); X= -0.25 (c). 
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(T-x)/Pe*10 3 


Fig. 4b - Temperature difference profiles between free surface and bottom wall 
for X= -1/12 (a); X= -0.15 (b); X= -0.25 (c). 


4. LINEAR STABILITY ANALYSIS 

To investigate for what values of thè parameters Re, Pr, M, G (or 
Re, Pr, X) thè Solutions (13) are stable, small disturbances are introdu- 
ced into thè flow. The equations for thè evolution of thè disturbances 
are obtained by writing each flow variable as a sum of a base flow quantity 
and a corresponding small disturbance: 

(14) (u, v, w,p, T) = (w, 0, 0,p, T) + (u',v f , w',p', T’); 

where thè overbar denotes thè base flow and thè apex thè disturbances. 
The equations governing thè initial growth or decay of thè disturbance 
are obtained by introducing (14) in thè Boussinesq System (3) and boun- 
dary conditions (4), (5), and dropping terms of second order in small 
quantity. thè diturbances may be expanded in thè form of normal modes: 

(15) (W, v', w’,p', T') = (u(z), v(z), w(z),p(z), è(z))e i(ax+0y)+y ‘ 

and thè amplitude functions (fife), v(z), w(z),p(z), ò(z)) are then governed 
by equations which take thè form, by eliminating thè symbol A in order 
to simplify thè notation: 
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(16a) (L-iaReG-yRe)u(z) = iap + RewDa 

(16b) (L-iaReu-yRe) v(z) = ifip 

(16c) (L-iaRea-yRe)w(z) = Dp-Gx> 

(16d) (L - /a Pe « - 7 Pe) d(z) = Re(uT x + wD T) 

(16e) Dw + iau + i(3v = 0 

(16f) u = v = w-Dù = 0 z = 0 

(16g) w = Du + Miad = Dv + Mi(3d = Dd = 0 z= 1 


where L=D 2 -k 2 \ D=—; k 2 = a 2 + /3 2 . 
dz 

T x is thè x-derivative of thè base-flow temperature field; in thè pre- 
sent study this derivative is equal to one but, for clearness sake, we re- 
tain it in thè equations. 

The parameters a and /3 are reai numbers in thè x and y directions 
and y = y r + i 7 , is thè complex frequency of thè disturbance, thè reai 
part of which is an amplification rate and thè imaginary part an oscilla- 
tion frequency. The problem (16) is a dispersion relation between a, 13 
and 7 , with Re, Pr, M, G as parameters; for fixed (a, 0) thè problem 
is an eigenvalue differential problem for 7 . According to thè linear theory, 
thè flow is stable if, for given Re, Pr, M, G, thè non-dimensional am¬ 
plification rate y r is negative for all values of a and (3. Conversely, thè 
flow is unstable if y r is positive for some values of a and /3. 

The States for which y r = 0 are called States of neutral stability and 
are characterized by 

(17) 7 r (a, 13, Re, Pr, X) = 0 

If at thè neutral state is 7 , = 0 thè instability arises with respect to sta- 
tionary disturbances and thè neutral state is called stationary; if results 
7 , ^ 0 we say that thè neutral state is oscillatory. 

Eq. (17) defines thè neutral stability Reynolds number 

(18) R N = R N (a, (3, Pr, X) 

Lastly, that neutral state for which thè neutral Reynolds number has a 
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minimum as a function of a and /3 is called criticai state, and thè corre- 
sponding Reynolds value is thè criticai Reynolds number: 

(19) R c = min R N (a,l3, Pr, X) 

a, (3 


4.1. Transverse and Longitudinal Modes 

The System (16) simplies for strictly two-dimensional transverse (/3 = 0) 
and longitudinal (a = 0) modes, thè former being independent of thè y 
coordinate, thè latter independent of x. 

For transverse modes thè v-equation decouples from thè set (16); its 
solution shows thè v-component of thè velocity to decay independently. 
The longitudinal modes depend only on y and z, but because of thè base 
flow in thè x-direction, thè w-component of thè disturbance velocity is 
not zero and thè motion is tridimensional. 

In stability analyses attention is often limited to transverse and longi¬ 
tudinal modes, and one does not consider stability with respect to tridi¬ 
mensional disturbances, as it is necessary when thè Squire theorem is 
not valid, as in thè present case. The above limitation can lead to erro- 
neous results, as our analysis has clearly shown. We shall underline this 
point in thè discussion of thè results. 


4.2 The Solution Method 

We solved thè eigenvalue problem posed in equations (16) by thè Ga- 
lerkin method after eliminating thè pressure and thè v-component of thè 
velocity by using thè continuity equation. The System results: 

(20a) [L 2 - ia Re (uL - D 2 u) - y ReL] w = k 2 G d 

(20b) [L-/aRew- 7 Re] ^u-^—Dwj =~^ R eDuw 

(20c) [L - ia Pew-yPe] $ = Pe(w T x +wDT) 

(20d) u=w = Dw = Dd = 0 z = 0 

(20e) w = D 2 w + Mk 2 d = Du + M/at9 = Z) $ = 0 z=l 
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The variables w, u, d, are expanded as 


(2 la) 

w = a,w,(z) 

i= 1 , . 

N 

(21 b) 

u = b i u i (z) 

/= 1 , 

..., N 

(2 le) 

d = C l l),(z) 

fc= 1 , 

..., N 


where thè Constant coefficients ai, b h c,, in generai complex, are unk- 
nown quantities, while thè w,, u,, are a-priori given function; in (21) 
thè summation Einstein convention is utilized. 

We assumed for w t , u h $,■ thè expressions: 

(22a) w i = z 2 (z— 1) C*_] (z) 

(22b) u i = zC‘_ l (z) 

(22c) tf, = C*_,(z) 

where C* is thè /-th modified Chebychev polynomial of first kind defi- 
ned in thè range 0 ^ z ^ 1. 

The constants a,, b,, Ci are determined by inserting thè (21) in thè 
equations (20) and then taking thè inner product with each trial func¬ 
tion over thè range 0 ^ z ^ 1; use of integration by part allows to ac- 
count for thè naturai conditions that are not satisfied by thè single trial 
functions. 

The result is a generalized complex algebraic eigenvalue problemi 
(23) Ax-yBx=Q 

where x T = (a u ..., a N , b u ..., b N , c u ..., c N ) = (a T , b T , c T ) is thè tran- 
sposed of thè column vector x of thè Galerkin coefficient a, b, c. The 
matrices A and B (Viviani 1988), are of dimension 3Nx3N, and they 
are reai in thè case of longitudinal disturbances (a = 0); for 13 = 0 thè pro- 
blem can be reduced to dimension 2Nx2N (see Appendix I). 

The eigenvalue problem (23) has been solved by using thè EIGZC su¬ 
broutine from thè IMSL library (Cowell 1984). 

For fixed Pr and X, thè neutral state (R N ) is calculated for each a and 
/3 by using thè secant iteration method to meet thè condition y r = 0, whe¬ 
re y r is thè reai part of thè principal eigenvalue, i.e. thè eigenvalue with 
thè greatest reai part; then a minimization on a and (3 determines thè 
criticai state (R c ). 
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5. RESULTS OF THE LINEAR STABILITY ANALYSIS 

The linear stability boundary has been defined calculating thè criticai 
Reynolds number R C = R C ( Pr, X), for values of thè parameter X varying 
in thè whole of its range of definition (-<»,+<») and then covering all 
thè cases of free convection (Marangoni, naturai, eombined). 

For thè Prandtl number, we analyzed thè range (10 4 , IO 3 ), thus co¬ 
vering both thè cases of liquid metals [Pr « 1] and of Silicon oils 
[Pr ^ 0(1)], often used in experiments performed in microgravity con- 
ditions, as well as thè limit cases Pr ^ 0 and Pr -> oo. 

In all thè cases analyzed, thè instability occurs always with respect to 
oscillatory disturbances. 


5.1 Marangoni Convection (X= ±°°) 

Fig. 5 shows thè criticai Reynolds number versus Prandtl, for thè ca¬ 
se of Marangoni Convection (X= ±oo). The instability occurs with re¬ 
spect to tridimensional disturbances, for any value of thè Prandtl number. 
The criticai curves for thè two cases X= ±oo coincide. The two cases 



Fig. 5 - Criticai Reynolds number (R c ) and product Pr R c versus thè Prandtl 
number (Pr) for Marangoni convection (X= ±oo). 
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X = ± oo correspond to an inversion of velocity and vertical temperature 
fields (see Fig. 2), but for X = - oo (surface tension increasing with thè 
temperature) thè flow is in direction of thè horizontal thermal gradient 
and for X = + oo it is in opposite direction. Nevertheless one must take 
into account that thè Marangoni stress of thè perturbation field is also 
inverted, so that thè physical effect that gives rise to thè onset of insta- 
bility is unchanged. From a mathematical point of view this result fol- 
low from thè property of invariance of thè System (20) under thè following 
transformation: 

(24) (M, u, w, i9, a, u,DT ) -» ( -M, - u, w, - ù, - a, - u,-DT) 

The case X = + oo (M= + 1) has been analysed by Smith (1982), who sol- 
ved thè eigenvalue problem by using a method of direct integration of 
thè equations (this method does not guarantee thè convergence to thè 
Principal eigenvalue and a numerical investigation of thè spectrum must 
be performed, see Lee and Reynolds 1967). Next table shows that thè 
criticai Reynolds number, for different values of Prandtl, obtained with 
thè two different methods, are in excellent accordance. 


Pr 

Re 

0.1 

1 

5 

Present study 

223.653 

116.028 

452.801 

Smith (1982) 

223.667 

116.022 

452.794 


We have also evaluated thè criticai value R c for Pr = 0, by assuming that 
R c remains finite for Pr -> 0. Then, thè value R c - R t (Pr = 0, X = ±oo) 
can be calculated by setting Pe = 0 in thè System (20), and thè correspon- 
ding matrices A and B of thè eigenvalue problem (23) are given in 
Appendix I. The computed value is R c = 28208. We consider now thè 
limit Pr -» oo. As Fig. 5 shows, for Pr -» oo PrR c ^ B (Constant) and 
then R c -> 0 as B/ Pr. 

The asymptotic form of thè System (20) for Pr -» oo, is obtained by 
setting PrRe = B and Re = 0. The equations (20a-c) become 

(25 a) L 2 w = k 2 Gd 

(25b) l[u-^Dw^= 0 

(L — iaBu — 7 B) $ = B(uT x + wDT) 


(25c) 
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and thè corresponding algebraic problem can be formally expressed as: 

(26) Uc = yBVc 

where c T =(c x ,... c N ) are thè Galerkin coefficients for thè temperature d. 
The expression for U is given by 

(27) U=L-BgC l D-iaB(Q + HF l NC l D) 

The different matrices that appear in (27) and thè matrix Fin (26) are 
defined as submatrices of thè matrice A and B of thè full problem (23) 
(see Appendix I). 

The asymptotic value for B = PtR c , calculated using thè (26) is equal 
to 397.8, whereas for Pr= IO 3 thè product Pr7? c is 395.3. The asympto¬ 
tic behaviour of R c for Pr ->• oo is then 397.8/Pr. 


5.2 Naturai Convection (\ = 0) 

The case of naturai convection has been studied by Hart (1983) who 
considered only longitudinal disturbances. A comparison between thè re- 
sults of thè present study and those obtained by Hart is shown in Fig. 
6, where thè criticai Reynolds number versus Prandtl is reported. The solid 



Fig. 6 - Criticai Reynolds number ( R c ) versus thè Prandtl number (Pr) for na¬ 
turai convection (X = 0). 
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line corresponds to Hart results, found here as particular case of thè ana- 
lysis performed with respect to tridimensional disturbances. For very low 
values of thè Prandtl number (Pr < 10 ~ 3 ) thè less stable «modes» are 
transverse (see Fig. 6). Increasing thè Prandtl number thè transition to- 
wards longitudinal modes occurs by means of tridimensional disturban¬ 
ces. For a Prandtl number dose to 0.5 thè longitudinal disturbances 
become very stable, so that thè analysis performed by Flart does not pro¬ 
vide for cases over this value. In fact, dose to thè value of Pr = 0.5 tridi¬ 
mensional modes again rise, joining thè curve of Hart with that of thè 
transverse modes, typical of thè instability for high values of thè Prandtl 
number. An analysis of thè various contributes to thè disturbance ener- 
gy (Viviani 1988) and of their dependence on thè Prandtl number for 
transverse, longitudinal and tridimensional disturbances, should be useful 
to explain thè type of disturbance with respect to which instability sets in. 

Then, not only did we correct thè stability curve of Hart, but also 
we pointed out that an analysis with respect to tridimensional distur¬ 
bances is essential if one is interested in thè correct definition of thè sta¬ 
bility boundary. An instability with respect to two-dimensional 
disturbances must rise from thè generai analysis, not characterizing a- 
priori thè kind of thè disturbance, but investigating thè wave number 
value at thè criticai point. For instance, for thè case, previously discus- 
sed, of Marangoni convection, since we found that thè instability always 
occurs with respect to tridimensional disturbances, thè assumption of 
two-dimensional disturbance should give rise to an over-estimation of 
thè criticai Reynolds number, for any value of thè Prandtl number. As 
Fig. 6 shows, R c approaches its limit value, for Pr -a 0, much more 
quickly than in thè case of Marangoni convection. The calculated value 
of R c , for Pr = 0 is equal to 7690. 

For Pr -> oo, R t 0. However thè behaviour is different from thè case, 
previously discussed, of Marangoni convection. The numerical results 
show that in this case R c decreases as Pr~ 1/2 . Therefore, in this case of 
naturai convection, thè proper limit for Pr -+ oo is Pr l/2 /? c ->const. 

The expression (25) for thè eigenvalue problem, for Pr->oo, is not 
more valid, because Pri? c .-» oo. The law of variation of /?,.-♦ 0 as Pr l/2 , 
for Pr-+ oo has been also found by Gill and Davey (1969) who studied 
instability of naturai convection along a vertical piate. 
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5.3 Combined Convection (X finite, different from zero) 

Fig. 7 shows thè criticai curves for three different values of thè para- 
meter X > 0. As X decreases thè criticai curves take thè shape peculiar 
to thè case of naturai convection (X = 0) (Fig. 6). However, we do not 
see yet thè appearance of transverse disturbances for very low Prandtl 
numbers. They become less stable than thè longitudinal disturbances for 
lower values of thè parameter X. This appears clearly in Fig. 8, where 
thè criticai curve for X = - 0.04 and thè stability boundary with respect 
to longitudinal modes are shown. From this figure we can note as thè 
instability to longitudinal disturbances, typical of X = 0 for Pr <1 (Fig. 
6), is replaced by an instability to tridimensional disturbances. Figg. 9 
and 10 show thè evolution of thè stability boundary towards thè typical 
case of Marangoni convection. 

The values R c (Pr = 0, X) are reported in Table 1. 

Consider now thè behaviour of thè criticai Reynolds number for 
Pr-* oo. As in thè cases of Marangoni and naturai convection, for any 
value of X, for Pr - oo R c -> 0. Except for thè case X = 0, where we 
found that R c -*• 0 as Pr ~ 1/2 , in all thè other cases investigated [ | X | € 
6 (10 2 ,oo)], for Pr-* oo, R C ^B(K). This behaviour is evident from thè 



Fig. 7 - Criticai Reynolds number ( R c ) versus thè Prandtl number (Pr) for com¬ 
bined convection with X > 0. 
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Fig. 10 - Criticai Reynolds number ( R c ) versus thè Prandtl number (Pr) for so¬ 
me values of X < 0. 


Table 1. Criticai Reynolds number R c ( Pr, X) and product R c Pr for various va¬ 
lues of X for Pr = 0, Pr = 0.0001, Pr = 0.001. 



Pr = 0.001 

Pr = 0.0001 

Pr 

= 0 

X 

R c 

R c Pr 

R c 

R c Pr 

R c 

R c Pr 

+ oo 

2119 

2.119 

13041 

1.3041 

28208 

0 

1 

3769 

3.769 

22714 

2.2714 

71769 

0 

0.5 

5089 

5.089 

30285 

3.0285 

186112 

0 

0.25 

7026 

7.026 

40960 

4.0960 

327768 

0 

0.08 

10911 

10.911 

63628 

6.3628 

97271 

0 

0.04 

13134 

13.134 

33552 

3.3552 

33564 

0 

0.01 

10112 

10.112 

10015 

1.0015 

10009 

0 

0 

7743 

7.743 

7694 

0.7694 

7690 

0 

-0.01 

6603 

6.603 

6581 

0.6581 

6579 

0 

-0.04 

5364 

5.364 

5351 

0.5351 

5348 

0 

-1/12 

6212 

6.212 

6202 

0.6202 

6201 

0 

-0.125 

22631 

22.631 

22595 

2.2595 

22593 

0 

-0.15 

29176 

29.176 

29105 

2.9105 

29098 

0 

-0.25 

21708 

21.708 

121927 

12.1927 

124910 

0 

-0.5 

8588 

8.588 

57339 

5.7339 

88751 

0 

-1 

4880 

4.880 

30850 

3.0850 

52691 

0 

-OO 

2119 

2.119 

13041 

1.3041 

28208 

0 
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analysis of Table 2, where thè values of R c and Pr R c , for Pr = IO 2 and 
Pr = IO 3 , are reported for different values of thè parameter X, together 
with thè limit for Pr -» oo, i.e. thè B(\) values calculated by (26). Then 
thè asymptotic formula for R c is: 


(28) R c ~* [|X|6(10- 2 , oo)] 
Pr 


Since for X= ±10~ 2 R c stili decrease towards zero as Pr ', thè case 
X = 0 seems to have very particular properties. We believe that this case 
should be investigated thoroughly, studying, with suitable analytical me- 
thods, thè asymptotic characteristics for Pr ^ oo, for small values of thè 
parameter X. 


Table 2. Criticai Reynolds number 7? c (Pr, X) and product R c Pr for various va¬ 
lues of X for Pr= 100, Pr=1000, Pr = oo. 



Pr = 

100 

Pr = 

1000 

Pr 

= 00 

X 

R c 

R c Pr 

R c 

R c Pr 

R c 

R c Pr 

+ OO 

3.749 

3749 

0.395 

395 

0 

398 

1 

21.501 

21501 

2.165 

2165 

0 

2267 

0.5 

32.282 

32282 

3.249 

3249 

0 

3256 

0.25 

53.621 

53621 

5.390 

5390 

0 

5391 

0.08 

130.132 

130132 

13.060 

13060 

0 

13053 

0.04 

223.653 

223653 

22.450 

22450 

0 

22458 

0.01 

644.412 

644412 

68.840 

66840 

0 

64881 

0 

8729.301 

8729301 

2813.100 

2813100 

0 

« (*) 

-0.01 

198.591 

198591 

19.870 

19870 

0 

19872 

-0.04 

44.371 

44371 

4.436 

4436 

0 

19872 

-1/12 

32.023 

32023 

3.164 

3164 

0 

3159 

-0.125 

140.352 

140352 

14.160 

14160 

0 

14184 

-0.15 

44.463 

44463 

4.490 

4490 

0 

4490 

-0.25 

9.531 

9531 

0.958 

958 

0 

960 

-0.5 

5.323 

5323 

0.536 

536 

0 

537 

-1 

4.143 

4143 

0.419 

419 

0 

420 

-oo 

3.749 

3749 

0.395 

395 

0 

398 


(*) For X = 0 (naturai convection) is /? c Pr 1/2 -► const. 


A theoretical approach appears indispensable because, for high va¬ 
lues of thè product Pr R c , one cannot be certain of thè validity of thè 
numerical results. In fact thè spectrum for thè temperature perturba- 
tion (coefficients c,) could be so wide that, from thè numerical point 
of view, thè problem should be solved by means of direct integration of 
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thè equations (if thè spectrum is very wide there are numerical difficul- 
ties, for thè solution of thè algebraic problem (23), because one obtains 
matrices of great dimensions). 

The results of thè three cases previously analysed can be illustrated 
in thè piane (Re, X), for different values of thè Prandtl number. They 
are shown in Figg. 11 and 12. 

We can see thè presence of a peak of stability for negative values 
of X, near thè value X= -0.15 (remind what we said about tbis par- 
ticular value of X when we discussed thè base flow in section 3). For 
very low Prandtl, we can see a lack of stability, near thè point X = 0. 
This occurrence can be explained if one remember that, for very low 
values of Prandtl, thè instability is due mainly to thè Reynolds stress, 
so that thè stability strongly depends upon thè velocity profile (see 
Figg. 3 and 4). 

For Pr = 0.1 thè lack of stability disappears, because thè thermal con- 
vection prevails and then thè influence of thè velocity profile shape on 
thè onset of instability becomes less important. 

Increasing thè Prandtl number a strong peak of stability for X = 0 oc- 
curs. Concerning this, two remarks are needed. First, for Prandtl num- 
bers greater than or equal to one, thè Marangoni effect considerably 
contributes to thè disturbance energy also for values of X very dose to 
zero, so that when this effect vanishes, for X = 0, a strong increase of 
thè stability boundary can occur. On thè other hand, as already remar- 
ked, thè case of naturai convection presents particular behaviours and 
it needs of further investigations; in this regard an analysis of thè de- 
pendence on X of thè disturbance energy for high Prandtl number will 
be helpful. 

In Figg. from 13 to 16 thè domain of linear stability is illustrated, 
for a more useful reading, in thè piane ( R M , R G ), where R M = V m f V r 
and R g = Vg/V v are thè conditional Reynolds number, defined in terms 
of thè Marangoni and buoyant speeds, respectively. The ratio of these 
two nondimensional numbers is just equal to thè parameter X. The cur- 
ves in thè piane ( R M , R c ) show thè peaks of stability and instability al¬ 
ready commented in thè piane (Re, X). 

The area inside these diagrams is thè domain of linear stability for 
thè basic motion under study. In fact, if we assume as reference velocity 
thè diffusive viscous speed V v , then R M and R G should appear as inde- 
pendent parameters in thè field equations and in thè base flow, instead 
of thè Marangoni and Grashof numbers, respectively. Thus thè domain 
of stability can be defined varying one parameter, for fixed values of 
thè other. The results obtained in this way are equivalent to those obtained 




Fig. 11 - Criticai Reynolds number (/? c ) versus X for Prandtl numbers Pr = 0.001, 
reported thè asymptotic values for X = ±oo (Marangoni convection). 





Fig. 12 - Criticai Reynolds number ( R c ) versus X for Prandtl numbers Pi 
ted thè asymptotic values ofr X = ±oo (Marangoni convection). 
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Fig. 14 - Linear stability domain in (R M , R G ) piane for thè Pandtl number 
Pr = 0.01. R m = V m /V v \ R c = V g /C„; V m , V g and V, are, respectively, thè Ma¬ 
rangoni, thè buoyant and thè viscous diffusive velocities. 


in thè piane (Re, X). In fact, R m =M Re, i? c = GRe, and we have found 
just ( R m ) c = MR c , and ( R g ) c = GR c . 

In thè piane ( R M , R G ) thè results can be read either for fixed X (Bond) 
(along straight lines starting from thè axes origin) or for fixed values 
of thè Marangopni and gravitational effects. 

For Pr < 1 (Figg. from 13 to 15) and R M > 0, as R G increases, thè 
corresponding criticai value of R M decreases, i.e. thè gravitational field 
acts as an unstabilizing effect on thè flow due to surface tension gra- 
dients. Viceversa, for R M < 0 (buoyancy and surface tension convection 
adverse) thè gravitational effect is stabilizing, so that thè criticai values 
of R m and R c are greater than thè criticai ones corresponding to pure 
Marangoni and naturai convection. 

For high values of thè Prandtl number (see Fig. 16), corresponding 
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Fig. 15 - Linear stability domains in (R M , R c ) piane for Prandtl numbers 
Pr = 0.05 and Pr = 0.1. R M = V m /V„; R G = V g /V„; V m , V g and V„ are, respecti- 
vely, thè Marangoni, thè buoyant and thè viscous diffusive velocities. 

to Pr = 20), thè diagram is strongly twist because of thè peak of stability 
corresponding to R M ~0 (X = 0). However, for R M < 0 it is stili possi- 
ble to see thè typical peak of stability near thè value X= -0.15, already 
pointed out and discussed in thè piane (Re, X). 


6. ENERGY STABILITY ANALYSIS 

We shall now examine thè stability of thè basic state by thè method 
of thè energy. For a wide account of this method see thè monograph 
of Joseph (1976) and for late developments that of Galdi and Rionero 
(1985). 
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Fig. 16 - Linear stability domain in ( R M , R c ) piane for thè Prandtl number 
Pr = 20. R m = V m /V v \ R c = V ? /V„; V m , V g and V„ are, respectively, thè Maran¬ 
goni, thè buoyant and thè viscous diffusive velocities. 

Let (v, p, T) be thè solution of thè basic state, and let (v' ,p', T') 
be any other solution to thè Boussinesq equations (3) and boundary con- 
ditions (4), (5). Now let 

y = y' - v 


(29a) 
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(29b) d=T'-T 

(29c) 7 T=p'-P 

(v, ir. d) represents thè difference between thè disturbed and undisturbed 
States. Since both (v, p, T) and (v', p' , F') satisfy thè System (3), (4), 
(5), we obtain thè following non-linear System governing thè evolution 
of thè disturbances (y, ir, $): 


(30a) 

(30b) 

(30c) 

(30d) 

(30e) 


Rei — + v• V v + v • V v + v • V v 

la/- 


Pe — + v.Vi5 + vVr+v-Vi? 

la/-- - - - - 


+ V 7 r = V 2 y + Gì) i z 

= V 2 i? 


u = 


a z 


for z = 0 


^a^ a^ 

W= — +M — = — + M — = — = 0 

ax 3z a^ 


Since thè ranges of x and y are infinite, we must consider that all distur- 
bance quantities remain bounded for |x| -* oo, | y \ -> so that thè energy 
integrals exist. We shall confine our attention to disturbance functions 
(u, v, w, ir, i?) in thè class 5 

S-(/|j / d* < °°| 


where / is either 

a) periodic in x and y 

b) Fourier transformable in x and y 

c) Fourier transformable in either x or y and periodic in thè other. 
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In case (a), we define thè integrai over a volume R: 

</) - j f<k= J 2 "’ ” j 2 "” ^fdxdydz 


and thè surface integrai at z = 1 : 

(■ p2ir /a |>2ir//3 

</),= /= f{x*y,\)dxdy 


where / is any disturbance function and a and (3 are thè wave-numbers 
in thè x and y directions, respectively. 

In case (b) thè integrals will extend from - oo to + oo in thè coordi¬ 
nate x and y, and in case (c) they will extend over a wavelength in coor¬ 
dinate in which thè function is periodic and from - oo to + oo in thè 
coordinate in which thè function is Fourier transformable. 

We now form thè scalar product of (30a) with y, multiply (30b) for 
and integrate over thè volume R. We employ thè Gauss theorem and 
thè fact that y and y are divergence-free and, by accounting for thè 
boundary conditions on thè surfaces z = 0 and z= 1 , obtain: 

Energy Identity 


(3la) Re-= - (Vy: Vy)-Re(y-(V y) s -v} + G(y-i,è) 

dt 

-M{ud x +vù y ) l 


Entropy Identity 

(3 lb) Pe — = - (Vf V t?)-Pe(i?y-V T) 
dt 

where 


‘ <J r> 




(31c-d) 
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(3 le) 


(V V)*=- 


"0 o a z ' 
0 0 0 
J, 0 0 _ 


where thè subscripts denote differentiation and ~ transposition. 

By introducing thè coupling parameter (Joseph 1965) A, we define 
thè energy fucntional 


(32) £ A = /f+APr9 

The coupling parameter A must be greater than zero in order to insure 
that thè condition E A -> 0 is satisfied if and only if both K and 9 ap- 
proach zero. 

We can obtain an expression for thè rate of chance of E A by adding 
(3la) and (3lb): 


(33a) Re—--A + ReP + E 

dt 

where 

(33b) A= ( Vy: Vv) + A( V ù- Vi?) 

(33c) P= -<v-(V v)*-v)-PrA(tfv- V T)+-^{v-i z ù) 

(33d) 'L=-M{ué x + vé y ) l 


The flow is said to be asymptotically stable in thè mean if E A -* 0 for 
t -> 00. 

Defined, then, thè maximum problem (following Davis and Von Kerc- 
zek 1973): 


(34) 


fi = max 
r 


-A + ReP + E 
Re£ A 


T = [v(u, v, w), d | y, ù with continuous second partial 
derivatives 


V • v = 0, v(x, y, 0) = 0, w{x, y, 1) = 0, u, v, w, ù 6 3] 
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from (33a) it follows 


1 dE A n 
(35) — —f ^ fi 

E a dt 


so that thè flow is stable if 0 < 0. 

For fixed A, thè condition 0 = 0 determines thè Reynolds number R A 
such that for Re < R A stability is guarantied, according to thè previous 
definition. 

The condition 


dR A 

(36a) —^=0 

<3A 

determines thè optimal value A opl of thè coupling parameter A and thè 
optimal «stability boundary» 

(36b) R E = R Aopl 

The solution of thè problem (34) can be obtained via thè direct method 
of thè variational calculus, by applying thè Rayleigh-Ritz technique to 
thè following functional 

(37) F= -A + ReP + E-ORe E a + 2(tt V-v) 

where ir is thè Lagrange factor for thè condition V_-v = 0, when we 
does not work in a space of divergence-free vector. 

We can also solve thè problem (34), by defining thè Eulero-Lagrange 
equations 

(38) 8F(v, ò,tt) = 0 

and then solving thè resulting differential eigenvalue problem for 0 via 
thè Galerkin technique, as we have done for thè eigenvalue problem of 
thè linear theory. 

By applying thè variational procedure to (38) we obtain thè problem: 

(39a) ORey + Rey-(Vy) s +V tt + 

+ — Red ( APr V T —— /J = V 2 v 
2 \ — Re Ì 
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(39b) fìPed + — v- (APr V T-— ùì = V 2 d 

2A“ \ — Re 7 

(39c) V-v = 0 

(40a) u = v = w= — =0 for z = 0 


(40b) 


dd 

M 

dw 

du 

dz, H 

h 2A 

dz 

dz 

dv 

M 

dd 

= 0 

iZ 

2 

dy 


The last three conditions of (40 b) and thè last one of (40a) are thè natu¬ 
rai conditions satisfied by thè variational principle (34). 

The problem (39), (40) defines an eigenvalue problem for thè reai pa- 
rameter 0. The principal eigenvalue of thè problem coincides with thè 
solution of thè maximum problem (34). 

Let us now express thè disturbance function as: 
in case (a) 


/= Reai [f(z, a, /3) e Kax+m ] 
and in case (b) 

/= ( + ( hz,*,P)e i(aX+ ^dadp 


Separate thè variables or Fourier transform thè System (39)-40). In ei- 
ther case (a), (b), (c) thè only influence of thè x and y variation is to 
introduce thè wave-numbers a and (3, and thè System for thè perturba- 
tion amplitudes u(z), v(z), w(z), rr(z), $(z) will be, by eliminating thè 
Symbol A in order to simplify thè notation: 

(41a) (D 2 -k 2 ) u = QReu + iair+ wDu+ y Re?9APr 

(41b) (D 2 — k 2 )v = ti Rev + ifì-ir 




(41e) Dw + iau + il3v = 0 

(41 f) u = v = w = Dd = 0 for z = 0 

M M M 

(41g) w = Dd-\ - Dw = Du + — ia& = Dv-\ -/j6i? = 0 

v 2A 2 2 

for z= 1 


wher qD=—\ k 2 = a 2 + P 2 
dz 

If fì max is thè maximum eigenvalue of thè problem (41): 

(42) fl max = 0 max (a, 18 , A, Pr, X, Re) 

then thè condition fi max = 0 defines thè criticai Reynolds number i? 0 , a f° r 
fixed (a, /3, A): 

(43) R 0 A = R 0 A (a, (3, A, Pr, X) 

and thè criticai Reynolds of asymptotic stability for a fixed value of thè 
coupling parameter A is given by: 

(44) /? A = mini ? 0 A (a, /3, A, Pr, X) 

a,0 


Finally, thè optimal stability boundary determined by thè energy theo- 
ry, i.e. R E , is: 

R F = maxR A = max mini ? 0 a(«> /3, A, Pr, X) 

A A a, fi 


( 45 ) 
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In thè numerical solution of (41) it is sufficient that one search only thè 
first quadrant of thè (a - 0) piane in order to find R A because of thè 
following symmetry properties of /? 0jA for fixed (A, Pr, X): 

(46a) Z? 0 ,a(A, Pr, X, -u,(3) = R OA (A,Pr,\,a,l3) 

(46b) 7? 0>a (A, Pr, A,ce, - 0) = R 0 , a (A, Pr, X, a, /3) 

The relations (46) can be obtained, by noting that U is reai, by transfor- 
ming thè ststem (41) as follows: 

(47a) ( u, v, w, ir, tf,a,f3,Q) -> (- wf, vf,- wf,- 7r*,- ,-af .jS^Q,) 

(47b) (u, v, w, 7T,t?, ce, |6,12) -> (iq,- v,, w, , i„ -i^cq,/3,, fi,) 


where * denote thè complex coniugate. 

The numerical determination of R E must be done by calculating R A 
for various values of A and then taking its maximum value by, e.g. a 
polynomial interpolation of thè points near thè maximum. 

However, it is possible to obtain, by parametric differentiation (that 
allows to deduce a generai relation satisfied by thè eigenvalue problem 


solution when thè parameters are varied), thè optimal coupling parame- 
ter A opt , for each Pr and X, in terms of thè eigenfunctions which solve 
dR A 

thè problem at -=0. 

dA 


In our case, by putting Q = 0 in (39), (40) and then performing thè 


parametric differentiation, we obtain 


G (wi?) M {ùw z ) x 

(48) A opt p £ ( u $T x +wùT z ) Pe {udT x + wdT z ) 

where u, v, w, d are to be evaluated for A = A opt . 

We note that for 7^ = 0 and T z = - 1 (liquid layers heated from be- 
low) thè expression (48) reduces to that one reported by Lebon e Cloot 
(1982) if one assumes as reference velocity thè thermal diffusivity velo- 
city V a , as Lebon and Cloot have done. 

With thè expression of A opt supplied by thè (48), one can proceed ite- 
ratively in thè search of R E , by calculating R A for a starting value of A, 
e.g. A= 1, and then modifying A according to thè (48). 
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It must be underlined that thè use of parametric differentiation im- 
plies thè assumption of analyticity with respect to A of all thè depen- 
dent variables, that can be only verified a-posteriori. For example, Shir 
and Joseph (1968) found a singularity in thè dependent variables, mani- 
fested in thè prediction of A opt = — 1, in violation of thè restriction 
A < 0. 

So, in this first analysis by thè energy method of thè basic state under 
study, we shall confine our attention to A= 1, that coincides with thè 
originai formulation of thè energy method. 

The problem (41) has been reformulated as a generalized algebraic 
eigenvalue problem, in thè same way as we have done in thè linear sta- 
bility analysis, by using thè Galerkin method. The expansion of thè di¬ 
sturbale function is also thè same and thè matrix formulation is reported 
in Appendix II, together with thè expression of A opt in terms of thè Ga¬ 
lerkin coefficients for thè velocity and temperature amplitude functions. 

Because this is thè first energy stability analysis of a problem of na¬ 
turai and Marangoni convection with a non zero basic state, we have 
verified thè solution procedure and thè numerical code on some cases 
reported in Iiterature, both with static and dynamic basic state. For ma¬ 
jor details on this point see Viviani (1988). 


7. RESULTS OF THE ENERGY STABILITY ANALYSIS 

The results of thè stability analysis via thè energy method are illu- 
strated in thè ( R M , R G ) piane. Figg. 17, 18, 19, report thè stability do- 
mains for Prandtl numbers Pr = 0.001 and Pr = 0.01, Pr = 0.1, Pr=l, 
respectively. 

The diagram show in Fig. 17 is similar to that obtained by thè linear 
analysis (Fig. 13), but thè values of thè stability are about 5-6 times less 
than those one of thè linear stability boundary. 

As thè Prandtl number increases thè stability diagrams modify until 
they become symmetrical with respect to thè axis R M = 0 (Pr ^ 0.1 in 
Figg. 18 and 19). This occurrence clearly point out that thè energy me¬ 
thod is less accurate than thè linear stability analysis, as it is possible 
to deduce by comparing thè results shown in Figg. 18 and 19 with thè 
linear one of Figg. 14 and 15. Similar results from thè application of 
thè two methods are obtained only for very low Prandtl numbers, when 
thè stability is controlled by thè velocity field. As thè Prandtl number 
increases, i.e. thè relative importance of thè thermal convection with re¬ 
spect to thè momentum convection increases, thè results show that, 
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Fig. 17 - Stability domains in (R M , R c ) piane obtained by thè Energy method 
for Prandtl numbers Pr = 0.001 an Pr = 0.01. R M = V m /V v \ R G = V g /V„; V m , V g 
and V v are, respectively, thè Marangoni, thè buoyant and thè viscous diffusi¬ 
ve velocities. 

by using thè energy method, one completely loses thè detail of thè flow 
field. This point must be object of further studies. 

For Pr= 1 (see Fig. 19) thè stability boundary in thè case of Maran¬ 
goni convection is R M = 8.984, whereas thè corresponding linear criticai 
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Fig. 18 - Stability domain in (R M , R G ) piane obtained by thè Energy method 
for Prandtl numbers Pr = 0.1. R M = V m /V v \ R G = V g /V„; V m , V ? and V v are, 
respectively, thè Marangoni, thè buoyant and thè viscous diffusive velocities. 

value of R m is 116.028; in thè case of naturai convection thè stability 
boundary is i? c = 9.142, greatly different from thè value evaluated with 
thè linear theory, which is equal to 82123. A comparison between thè 
criticai Reynolds number R c and thè Reynolds number R E of asymp- 
totic stability in thè mean is reported in Tab. 3, for Marangoni and na¬ 
turai convection. 

Results are not reported for Prandtl numbers greater than one, because 
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Fig. 19 - Stability domain in ( R M , R G ) piane obtained by thè Energy method 
for Prandtl number Pr = 1. R M = V m /V v \ R c = V g /V„\ V m , V g and V„ are, re- 
spectively, thè Marangoni, thè buoyant and thè viscous diffusive velocities. 

they can be obtained by dividing thè result of Pr = 1 by thè actual value 
of thè Prandtl number. So, with thè energy method, it si possible to con¬ 
clude that for Pr-> oo, R E approaches zero, approximately, as 10 Pr 1 , 
for all values of thè parameter X. The above asymptotic behaviour also 
holds in thè case of naturai convection, whereas by thè linear theory we 
have found an asymptotic law as Pr~ 1/2 . The different asymptotic law 
for high Prandtl numbers supplied by thè liner and energy theories in 
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Table 3 - Comparison between thè results of linear theory and energy method 
for Marangoni and naturai convection. 


Pr 

MARANGONI CONVECTION 

NATURAL CONVECTION | 

R c 

Re 

R c 

Re 

0.001 

2119.000 

175.60000 

7743 

1394.0000 

0.010 

674.680 

144.30000 

4278 

780.3000 

0.100 

223.653 

85.88000 

1781 

91.3300 

1.000 

116.028 

8.98400 

82123 

9.14200 

10.000 

27.430 

0.89310 

27183 

0.91400 

100.000 

3.749 

0.08920 

8729 

0.09140 

1000.000 

0.3953 

0.00892 

2813 

0.00914 


thè case of naturai convection, have been also encountered in thè case 
of buoyant boundary layer on a vertical fiat piate (Gill and Davey 1969: 
linear theory => Pr' 1/2 ; Dudis and Davis 1971: energy theory => Pr _1 ). 
The asymptotic law, for Pr -» oo, of R c and R E are summarized in Tab. 4. 


Table 4 - Asymptotic laws for high Prandtl numbers. 


TYPE OF 
CONVECTION 

THEORY 

LINEAR 

R c 

ENERGY 

Re 

NATURAL 
\ = 0 


C Pr 1/2 

m 10 Pr 1 

MARANGONI 

X = ± OO 


397.8 Pr 1 

m 10 Pr 1 

COMBINED 

X * 0 


B (X) Pr " 1 

= 10 Pr -1 


for X * 0 and |X| < IO 2 ; Pr^o° 


Rc 0 
Rc Pr -> B(X) 


for X = 0 


Rc ^ 0 
Rc Pr 1/2 C 


B(X) has been evaluated by thè asymptotic expression of thè eigen- 
value problem obtained by setting Re = 0 and RePr = B. The Con¬ 
stant C has not yet been evaluated. 
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The ramerkable differences, both qualitative and quantitative, of thè 
results obtained by using thè liner theory and thè energy method seem 
to demonstrate that thè latter is not suitable for thè search of sufficient 
stability conditions, for thè class of flows here considered, because they 
are conservatively low. It will be necessary to solve thè non-linear equa- 
tions, both by numerical and spectral methods, in order to validate and/or 
correct thè results of thè linear theory. However, it must be point out 
that thè present analysis has been performed by assuming thè coupling 
parameter A = 1, and further studies are needed to analyze thè influence 
of this parameter on thè stability boundary. 


8. CONCLUSIONS AND FINAL COMMENTS 

A linear and energy stability analysis of Marangoni, naturai and com- 
bined free convection in shallow cavities has been performed. 

The «base flow», thè stability of which has been analyzed, is a dose 
form solution depending on thè non-dimensional parameters Marango¬ 
ni, Grashof, Reynolds and Prandtl numbers, under thè hypotheses of 
planar and undeformablre free surface, and of surface tension linearly 
decreasing or increasing with thè temperature. 

The linear stability analysis has been performed with respect to tridi- 
mensional disturbances, and thè resulting eigenvalue problem, for thè 
evaluation of thè criticai Reynoldes number, has been solved with thè 
spectral Galerkin method. 

The criticai Reynolds number has been calculated in thè whole range 
of definition of thè dynamic Bond number, and for Pr number ranging 
from IO 4 to IO 3 as well as for thè limiting cases Pr ->0 and Pr -»■ oo. 
For all values of Prandtl and Bond numbers we always found an oscil- 
latory instability. 

The linear stability boundary, defined by thè criticai Reynolds num¬ 
ber, has been shown in thè Reynolds-Prandtl piane, where thè criticai 
curves are reported in terms of a parameter X, having thè meaning of 
thè inverse of thè Bond number. 

In thè case of Marangoni convection, we confirm thè results of a pre- 
vious study of Smith (1982); for naturai convection, our analysis with 
respect to tridimensional disturbances allowed to correct and to extend 
thè results of Hart (1983) whose study was limited to two-dimensional 
longitudinal modes. Asymptotic laws for thè criticai Reynolds number, 
for high value of thè Prandtl number, has been established. In thè case 
of naturai convection we found that R c decreases as Pr“ 1/2 ; for thè 
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other cases of Marangoni and combined convection R c decreases as 
B(\)/Pr, where thè function B(X) has been numerically evaluated by sol- 
ving thè asymptotic form of thè eigenvalue problem for Pr -> <». 

The criticai curves has been also reported in thè piane (Re, X) for va- 
rious values of thè Prandtl number, and this representation evidentiates 
a strong peak of stability, for high values of thè Prandtl number, at thè 
point X = 0 (naturai convection). In view of this result, and of thè asymp¬ 
totic behaviour of thè criticai Reynolds number, for Pr -»• oo, in thè case 
X = 0, further studies on this subject are needed. 

The linear stability domains are shown in thè piane (R M , Rg)> where 
thè results can be read either for fixed Bond number or for fixed values 
of thè Marangoni and gravitational effects. 

By thè energy method we have determined thè domains of stability 
in thè piane (R M , R c ) and also evaluated thè asymptotic behaviour of 
thè stability boundary for high Prandtl numbers, by assuming thè cou- 
pling parameter A = 1. For all thè values of thè Bond number R E de¬ 
creases for Pr->oo as Pr” 1 . 

The results obtained seem to be highly conservative, i.e. thè subcriti¬ 
cai instability range is very large and, then, thè energy method appears 
to be not appropriate to define sufficient conditions of stability for thè 
basic state herein considered. The solution of thè non-linear equations, 
both by numerical and spectral methods, will be necessary to confirm 
this conclusion. Further studies are also needed in order to evaluate thè 
influence of thè coupling parameter on thè stability boundary. 
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APPENDIX I 


Matrices of thè algebraic eigenvalue problem of thè linear stability 

The matrices of thè problem (23) are given by: 


C OD 


M 0 0 

ReE f D 

- /aRe 

— TV P 0 



Re - - 

PeG Pe^ L 


0 0 Pr Q 


R eR 0 0 


0 0 0 


aRe 


0 Re5 0 


\k 0 0 


l ^ 


0 0 PeF 


Oli 

Oli 

Oli 


where thè submatrices are reai and of dimension TVxTV; thè diagonal 
submatrices are symmetrical; 0 denotes thè nuli matrix; C, D, E, F, L, 
are decomposed as: 


C = k*C l + 2k 2 C 2 + C ì 
D= -k 2 GD l -k 2 Ml ) 2 



F= -k 2 F l -F 2 
L= -k 2 V-L 2 

N= -N' -^r- N 2 

m m k 
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The coefficient of thè above submatrices are defined by: 


cd 

= {w m 

,W $ )i 

Cd 

= ( Dw„ 

„ Dw n >; 

cd 

= (D 2 

w m , D 2 w n } 





r — ( W m 

.finti 

E>d 

, = (Dw„, 

„ Di ?„), 

Ed 

= (Duu m , W n )l 

Ed 

= ( uu m , 

Dw n )\ 

Fl.n 

= {u m . 

u„); 

Ed 

= ( Du m . 

, Du n )\ 

G mn 

= -(DT & m , w„). 

H m ,n 

= -(T x d m , u „; 

E' m ,„ 

= (d 

finì ì 

Ed 

= {E>fi m , 

, Dù n )\ 

M m , 


„{uL-D 2 u) 

™n) 



Nd, 


, Dw„); 

Nd 

, = (Du m 

„D 2 w n ; 

Pm,n 

= (0U, 

• n > Un ti 

Qm,n 

= {ùfi m 

.finì 

Rm.n 

«K 

, Lw„); 

S m ,n 

= (u m , u, 

ni 

V m , n 

= (fi m . 

finti 

K,„'„ 

= ( u m , Dw„) 


I Eì 


In thè case 13 = 0 thè matrix A can be expressed in thè following form: 


C 0 D 


M 0 0 

Rei? E 0 

- ia Re 

— N P 0 



Re = = = 

PeG 0 L 


Pr 

—- J 0 PrQ 



a ‘ 


where J m n = ( T x § m , Dw n ). 
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The problem can then be reduced to dimension 2Nx2N: 


C D 


m o 


- /«Re 

Pr 

PeG L 


— J Pr Q 

ss; 


a 2 - 


ReF 0 
0 PeF 


For Pr = 0 (and Re finite) thè solution for thè temperature field is & = 0. 
The matrices A and B become: 


c o 


m o 


- /aRe 

ì 

Re£ F 


- N P 



Re = = 


ReR 0 


o g 


. aRe 


0 ReS 


k o 
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APPENDIX II 


The algebraic eigenvalue problem of thè energy stability analysis 

The algebraic eigenvalue problem for thè evaluation of thè criticai 
Reynolds number R E of thè energy method is expressed as follows: 

(II.l) A E x=tiBx 

where B is thè same as in (1.2), while thè matrix A E is: 


C ReX D E 


M E O Z 


ia Re 

1 

R qE e E ReJ 


- N O O 


2 

Re = = = 

G e H e l 


o o o 


and thè submatrices C, F, L, N are thè same as in (1.1). 

The elements of thè matrices X, Y, Z, E e , M E , are defined by: 

k 2 APr 

X m , n = ~ u„Du) Y m n = - 13 2 — (u m , # n T x ) 

(II.3) Z nhn = AYr{w m ,Dd n T x ) E E _„= {u m , w n Du) 

M E hn = - (w m , Dw m Du) - (w m , w„D 2 u) 

while thè matrices D E , G E , H E , are decomposed as: 


U E =\ i~Gk 2 D l +Mk 2 D 2 + APr Re k 2 D 3 ) 


G e = 


PrRe ^ 
2 


2A = 2A 


P 2 


H 2 = 


Pr Re 
2 


H 


(H.4) 
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where D\ D 2 , G, H, have been already defined for thè linear problem, 
see the(1.4). For thè sake of clearity, we report again thè elements of 
these matrices together with thè elements of thè matrix D 3 : 

D' m>n = {w m , d„) D\ n = ( Dw m , #„) 

(11.5) Dl„ = (w m , § n D T) 

G m n = - {§ m , w n D T) H m n = - {ù m , u„ T x ) 


The optimal value (48) of thè coupling parameter A is given by (thè su- 
perscript * denotes thè complex coniugate): 




Problemi non lineari di statica delle strutture 


Placido CICALA* 


Summary. Three sample nonlinear problems concerning elastic structures 
show that thè perturbation analysis of postcriticai stability is not exhausti- 
ve: exploration of thè nonlinear field by a series of trajectories may reveal 
importantphenomena that investigation on thè neighborhood ofbuckling 
cannot detect. 


Introduzione 

Nell’analisi non lineare delle strutture è consueto seguire per passi l’e¬ 
voluzione della deformazione a partire dallo stato indeformato o da uno 
di piccole deformazioni: ciò è opportuno per non «perdere il filo» nel 
groviglio che possono formare le linee corrispondenti a varie sequenze 
di configurazioni, nello spazio delle variabili rappresentative dello stato 
di tensione e deformazione; nel caso di corpi non elastici ciò permette 
di usare le leggi incrementali che esprimono il comportamento di tali ma¬ 
teriali. Quindi le «traiettorie» sono individuate mediante una successio¬ 
ne di punti corrispondenti a valori, convenientemente intervallati, di una 
variabile scelta come «parametro di percorso». 

L’analisi del primo tronco della traiettoria, dall’origine alla prima, 
eventuale 0 ’ singolarità (punto critico), in molti casi può essere lineariz¬ 
zata. La linearità deve essere abbandonata nell’esame della singolarità, 
di solito costituita dall’incrocio con altra linea (percorso deviato). Per 
l’esame del comportamento del corpo in quella situazione si possono ado- 


* Dip. di Ing. Strutturale Politecnico di Torino, C.so Duca degli Abruzzi, 24,10129 Torino, Italia. 

(1) La singolarità non esiste sui percorsi relativi alla struttura con imperfezioni di forma o in 
presenza di carichi atti a provocare il tipo di deformazione relativa al cammino deviato: dunque 
la biforcazione al primo carico critico costituisce una situazione eccezionale nella realtà, ordinaria 
negli schemi della teoria. 
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perare per le variabili rappresentative sviluppi in serie di potenze del¬ 
l’incremento A del parametro di percorso rispetto al valore alPincrocio: 
nelle equazioni complete (ossia non ridotte come sul percorso fondamen¬ 
tale) i termini in A conducono a determinare gli incrementi primi delle 
variabili «incipienti» e i valori base di quelle preesistenti; così i termini 
in A 2 precisano gli incrementi secondi delle incipienti e gli incrementi 
primi delle preesistenti. Così di seguito le successive potenze individua¬ 
no l’andamento dei diversi rami, con le restrizioni note degli sviluppi 
in serie di potenze. È nota l’importanza di questa indagine, detta per- 
turbativa, e la rilevanza dei risultati. Tuttavia l’attuale disponibilità di 
potenti strumenti di indagine numerica suggerisce di varcare la soglia 
e di avanzare nel campo non lineare, come coraggiosamente per allora, 
fece Eulero per i problemi di «elastica». Qui, appunto, si accenna ad 
alcuni risultati che non sarebbe stato possibile ottenere con l’analisi per- 
turbativa, il cui uso, negli stessi campi, era risultato non meno laborio¬ 
so ma meno esauriente. 


La parete sottile cilindrica soggetta a compressione assiale 

Nel 1946 Welter <2) pubblicò risultati di prove su pannelli cilindrici 
compressi assialmente, con imperfezioni iniziali artificialmente create. 
In particolare trovò che, se un pannello ha curvatura iniziale maggiore 
di quella imposta dalle guide nelle quali è posto in compressione, il suo 
comportamento migliora notevolmente. Su contratto della Cornell Uni¬ 
versity venne sviluppata’ 3 4 ' un’indagine teorico-numerica sulla questio¬ 
ne accennata e si trovò conferma dei risultati di Welter. Quello studio 
venne riesumato in un recento lavoro (4) e completato in parte. Qui si ri¬ 
feriscono -alcuni risultati. 

Nella fig. 1, in alto, è tracciato il grafico tensioni deformazioni per 
il cilindro sottile compresso assialmente. In ascissa sono riportati gli ac¬ 
corciamenti relativi - e moltiplicati per il modulo elastico E, in ordina¬ 
ta la tensione media - a: le due co ordinate s ono riferite alla tensione 
critica a c = -2 cE/R dove è c = h/V 12(1 - v 2 ), h lo spessore, R il raggio 
del cilindro medio, v il coefficiente di Poisson. Per il cilindro perfetto si 


(2) G. Welter, J. Aero. Sci., 13, 593-596, (1946). 

(3) P. Cicala, Quart. Appi. Math., Voi. IX, 275-293, (1951). Da qui sono ricavate le figg. 1 e 3. 

(4) P. Cicala, Vecchi studi e nuove prospettive. Omaggio a Ceradini, Dpt. Ingegneria struttu¬ 
rale e geotecnica. Università di Roma, 1988. 
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arriva al punto critico O per il percorso fondamentale o = Ee. Da questo 
punto si diramano due percorsi deviati: quello «buono» che prosegue 
nella direzione di arrivo con cammino poco curvo e quello «cattivo» dove 
si ha il noto cedimento di scatto con ritorno indietro (snap back). Que¬ 
sto diagramma è ottenuto assumendo per gli spostamenti w normali al¬ 
la parete una combinazione di due autofunzioni, quella a scacchiere di 
losanghe con fattore cw, e quella assialsimmetrica con fattore cw 2 (5) . 
Nella parte inferiore della stessa figura sono tracciati i due percorsi OA, 
OB nel piano di coordinate w u w 2 , limitato a w l > 0: nel semipiano 
Wj < 0 esiste il grafico simmetrico con altri due percorsi equivalenti. 
Nella medesima figura sono illustrati i risultati relativi al cilindro im¬ 
perfetto con deformazioni iniziali definite da valori assegnati w/, vv 2 ' dei 
due fattori, indicati con circoletto sul grafico inferiore. Le stesse lettere 
sono usate in ambedue i grafici a indicare punti e curve corrispondenti. 
Si constata che i percorsi con origine situata nell’angolo compreso fra 
la curva OA e l’asse + w 2 si accostano al ramo cattivo OB. Quelli con 
origine nell’angolo supplementare si dirigono verso l’asse - vv 2 . Questo 
tipo di evoluzione, che conduce a una legge a(e) meno svantaggiosa, è 
scartato se si assume una forma più generale per gli spostamenti. Infatti 
se si introducono per gli spostamenti assialsimmetrici due componenti 
in quadratura in modo che possa verificarsi uno slittamento assiale di 
queste onde rispetto a quelle a scacchiera, la possibilità enunciata scom¬ 
pare. Ciò risulta evidente dalla fig. 2 dove, nel grafico inferiore, nel se¬ 
mipiano sotto l’asse è rappresentata l’evoluzione delle due componenti 
in quadratura. La traiettoria g' che parte da un punto assai poco di¬ 
stante dal piano tv,, w 2 vira rapidamente per accostarsi alla linea a; la 
relativa caratteristica a(e), tracciata nel grafico superiore, rivela il cedi¬ 
mento di scatto, mentre la traiettoria contenuta nel piano w u w 2 (linea 
punteggiata g, corrispondente alla punteggiata g nel grafico a, e) mo¬ 
stra un buon comportamento della struttura. Dunque, per il guscio ci¬ 
lindrico indefinito, con completa libertà nella evoluzione delle 
deformazioni, per qualsiasi germe di deviazione dal cilindro circolare 
perfetto, la deformazione si evolve verso il comportamento di minore 
efficienza strutturale. 

L’indaginesvolta per il pannello delimitato da due generatrici di lun¬ 
ghezza ira^fcR e da due archi di direttrice di sviluppo ir^^fcR, appog¬ 
giato su questo contorno, si ottennero risultati qualitativamente analoghi. 
La fig. 3 mostra, per il pannello perfetto, 4 percorsi deviati m u m 2 , m 3 , m A 


(5) Per le specificazioni sulle forme di distribuzione degli spostamenti si rinvia alle 
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che partono da due punti critici marcati con frecce nel grafico superio¬ 
re. Per pannelli con deformazione iniziale rappresentata da punti pros¬ 
simi alla traiettoria migliore m, si manifesta la tendenza a deviare da 
quella caratteristica: però il percorso n u prima di raggiungere il carico 
limite presenta un ampio tratto di eccellente comportamento, che può 
raggiungere il limite d’elasticità; quindi questa deformazione previa con¬ 
ferisce al pannello la massima robustezza. Le potenzialità degli strumenti 
numerici oggi disponibili, enormemente superiore a quella dei mezzi uti¬ 
lizzati nel 1950, può permettere un controllo delle favorevoli possibili¬ 
tà, sì da rimuovere il dubbio che queste svaniscano con l’introduzione 
di un maggior numero di gradi di libertà per gli spostamenti. 


Flessotorsione delle funi metalliche con formazione di un cappio 

Le funi metalliche, in conseguenza della loro costituzione per avvol¬ 
gimento elicoidale di fili e trefoli, hanno il difetto di attorcigliarsi quan¬ 
do si diminuisce lo sforzo di trazione che esse sono destinate a sopportare, 
in presenza di torsioni originate dal loro avvolgimento. In queste situa¬ 
zioni si può formare un cappio che, quando la fune è nuovamente mes¬ 
sa in tiro, può originare una frattura. Il fenomeno è particolarmente 
pericoloso quando le funi sono adoperate per sollevare pesi dal fondo 
del mare, perché in tal caso l’inconveniente non può essere corretto con 
intervento manuale. L’interesse tecnico di tale questione è messo in evi¬ 
denza da una serie di studi promossi dalla Marina degli Stati Uniti. 

Il carico critico della verga soggetta a trazione e torsione venne deter¬ 
minato un secolo fa da Greenhill. Posto che la verga abbia la medesima 
rigidezza flessionale EI per qualsiasi giacitura di flessione e che all’e¬ 
stremo 1 di un tronco lungo L sia applicata una forza Ti è un momen¬ 
to Mi , che - Ti , -Mi agiscano all’estremo O, essendo / un versore 
parallelo a 01, si ha la relazine 

(1) y = x 2 -l 

essendo y = TL 2 /t 2 EI, ossia lo sforzo di trazione rapportato al carico 
di Eulero, x=ML/2rEI il momento torcente riferito al momento criti¬ 
co per trazione nulla. Tale momento 2-kEI/L, non raggiungibile dalla 
verga omogenea, si trova nel campo di impiego delle funi metalliche a 
causa della loro esigua rigidezza flessionale. Per esse, dunque, ha prati¬ 
co interesse un ampio tratto della curva rappresentativa della (1), linea 
1 della fig. 4. Tale parabola costituisce il limite inferiore del campo di 
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stabilità della verga. Per rappresentare la deformazione incipiente al punto 
critico, si stabiliscono coordinate cartesiane ortogonali £,, £ 2 > $3 di cui 
la prima è misurata lungo l’asse iniziale della verga. La deformazione 
critica dell’asse della verga può essere posta in forma complessa: 

( 2 ) ^2 + i^ì = A (e tx+ i)i: -e (x ~ 1)f ) 


dove £ è la variabile immaginaria i%% x /L che assume i valori O e iir al¬ 
le due estremità. Il coefficiente A è un complesso arbitrario. 

L’arco di parabola sotto l’asse x rappresenta situazioni in cui il cari¬ 
co d’estremità y < 0 esercita una compressione. Per x = 0 le deforma¬ 
zioni finite della verga sono calcolabili secondo Eulero. In una recente, 
interessante nota (6) Rosenthal determinò l’evoluzione delle deformazioni 
a partire dalle configurazioni anzidette sotto l’azione di momenti M cre¬ 
scenti. Nella nota citata si considera T come quantità modulare, mentre 
qui si attribuisce un segno a T, sicché la forza Ti ha verso opposto a 
i quando - T=P > 0. Il versore fisso i corre dall’estremo 0 della ver¬ 
ga verso 1 quando il corpo è indeformato e non cambia quando i due 
punti si incrociano perché si crea il cappio (7) . 

Definito, come fa Rosenthal, l’angolo 7 compreso fra la forza agen¬ 
te su uno dei punti terminali e la tangente alla linea d’asse con verso 
uscente, si pone 01 = 71-7. Lungo l’arco y < 0 della parabola è o = 0: 
la forza P ha verso entrante. Sull’asse x = 0 i punti y < - 1 rappresenta¬ 
no le soluzioni di Eulero con a crescente da zero a ir mentre - y va da 
/ a 00 . Nella fig. 4 è tracciata presso il punto a la deformata per a = 80°, 
corrispondente a y= - 1.2939 e presso il punto b la deformata per 
a = 160°, relativa a y= -4.030. In questo secondo caso, come per ogni 
y < -2.18, appare il cappio della linea d’asse. Rosenthal trovò la linea 
3, segnata a tratti in fig. 4, come luogo dei punti per i quali si verifica 
quella intersezione. Se, tenendo y costante, si incrementa x da zero, l’an¬ 
golo a cresce portandosi a tt sulla linea 2 di equazione y= 1 -x 2 , tra- 
ciata con tratto e punto sulla figura. Si nota che nelle configurazioni 
rappresentate da quella Enea risultano scambiate le posizioni dei punti 
0, 1 rispetto al versore f, sicché la forza P > 0 genera trazione, come 


(6) F. Rosenthal, Greenhill’s formula and thè mechanics of cable hockling, Nav. Res. Lab. 
Rep. No. 7940. 

(7) Precisato cosi il semipiano dove il punto rappresentativo va riportato nella fig. 4, si evita 
la sovrapposizione dei due semipiani che appare nella nota citata. 
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mostra il dettaglio relativo al punto c che si raggiunge da a nell’evolu¬ 
zione indicata. Poiché le coppie M non hanno subito cambiamento di 
segno, la situazione di carico si presenta come è mostrato nel dettaglio, 
messo a raffronto con quello relativo al punto C, simmetrico di c rispet¬ 
to all’asse x. 

Le conclusioni di questa indagine risultano chiare dalla fig. 5 grazie 
alla rappresentazione adottata nella fig. 4. Se la fune viene messa in tra¬ 
zione (punto A) e poi assoggettata ad un momento inferiore a quello 
di Greenhill per trazione nulla (punto B, con x < 1), qualora si ribassi 
il carico T sino a cambiarne il segno, sussiste la configurazione diritta 
sino all’intersezione con l’arco inferiore della parabola di Greenhill. Un 
ulteriore aumento della forza -T=P dà luogo alle inflessioni trovate 
da Rosenthal, con formazione del cappio all’intersezione con la curva 
3. È probabile che lo stesso fenomeno si presenti quando si raggiunga 
lo stato critico nel ramo della curva di Greenhill con ordinate positive, 
ossia quando si segua una traiettoria discendente da un punto C (fig. 
5) di ascissa ^ 1, e si parta dalla deformata critica che nel piano norma¬ 
le al versore i ha la forma cicloidale indicata dalla (2). 

Gli stati critici per la verga soggetta a torsione e trazione (o compres¬ 
sione) sono stati determinati per varie condizioni di vincolo e si sono 
trovati grafici aventi forma simile alla parabola di Greenhill. Una di¬ 
sposizione di vincolo attuabile con relativa facilità per controlli speri¬ 
mentali è stata considerata da Cicala e da Algostino <8) : due cerniere 
cilindriche ad assi paralleli, la cui distanza può essere variata, sono col¬ 
legate alle estremità della verga con un accoppiamento rotoidale che per¬ 
mette l’introduzione controllata di sforzi torcenti. Per questo tipo di 
vincolo si è esaminata la stabilità <9) . Un’altra disposizione di vincolo at¬ 
tuabile, che consente un maggior grado di libertà agli estremi della ver¬ 
ga, è costituita da giunti cardanici. Per questa situazione, per la verga 
soggetta a pura torsione, lo stato critico venne determinato da Koiter <10) 
e la stabilità postcritica fu esaminata da Cicala 00 : in questo caso e in 
presenza di carico assiale si sono trovati valori dei momenti critici net¬ 
tamente inferiori ai valori che si ricavano dlla (l) 12 . In quella nota si è 


(8) P. Cicala, Helicoidal buckling of an elastic rod, Meccanica, 1985, p. 124-126. F. Agosti¬ 
no, Stati crìtici della verga soggetta a trazione e torsione, AIMETA Res. Rep. No. 12, 1985. 

(9) D. Bigi, P. Cicala, Elastic stability of thè rod subject to tension and torsion, AIMETA Res. 
Rep. No. 13, 1985. 

(10) W.T. Koiter, Buckling of a flexible shaft under torque loads transmitted by Cardon joints, 
Ing. Archiv, 1980, p. 369-373. 

(11) P. Cicala, Sulla stabilità elastica della verga soggetta a torsione ultracrìtica, Atti d. Acc. 
d. Scienze Torino, 1986, p. 98-116. 

(12) F. Algostino, D. Bigi, P. Cicala, Nonlinear analysis of elastic round rods. Meccanica, 1987, 
p. 203-209. 
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preparata la formulazione per l’esame di spostamenti finiti, ma non si 
è condotta la relativa indagine numerica. L’instabilità postcritica con¬ 
statata nei lavori menzionati fa ritenere che imperfezioni della fune pos¬ 
sano innescare l’attorcigliamento in uno stato subcritico: allora la 
traiettoria che scenda C, prima di arrivare alla linea degli stati critici, 
porterebbe a deformazioni che compromettono la resistenza: la prose¬ 
cuzione di queste ricerche sarebbe quindi opportuna. 


Grandi deformazioni piane di un corpo iperelastico 

Per corpi molto deformabili come la gomma, i quali, sotto sollecita¬ 
zioni statiche, possano ritenersi dotati di elasticità e isotropia, la for¬ 
mulazione di Mooney Rivlin fornisce l’interpretazione più semplice e 
plausibile del loro comportamento per grandi deformazioni. L’energia 
di deformazione viene espressa come funzione lineare dei primi due in¬ 
varianti di Cauchy: la condizione di invariabilità del volume definisce 
il terzo invariante. Nel caso di moto piano, che qui è considerato, la teoria 
si semplifica ulteriormente coincidendo con quella detta «neohookean». 
Nel riferimento ad una terna fissa di vettori unitari ortogonali i a 
(a = l, 2, 3), detto 

(3) P=Lx a 

il vettore di posizione nel corpo deformato e introdotta la matrice di Ja- 
cobi p a0 = dx a /dt; iH mediante le coordinate cartesiane £„ del corpo defor¬ 
mato, l’energia specifica, a meno di una costante additiva, è scritta 

(4) W=np nii p, i!i 

con la costante d’elasticità g. Le componenti del tensore di tensione di 
Lagrange sono ottenute per derivazione di W rispetto alle quantità p llè : 
poiché queste sono legate dal vincolo interno di incompressibilità, si ag¬ 
giunge all’energia la variazione di volume con il moltiplicatore di La- 
grange -p che introduce nel tensore degli sforzi un componente sferico, 
con il carattere di pressione idrostatica (funzione della | a ). 

Nel problema in esame, assunte A' = L = ib come coordinate mate¬ 

riali, si scrivono per il corpo deformato le espressioni 

(5) x,/A = x+U, x 2 /A = \y+V, p/A = \ + P 

(6) U= £ Uj sin jx, V=T, Vj cos jx, P=LPj cos jx 
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essendo A, X costanti: le sommatorie (6) si estendono agli interi j da ze¬ 
ro a n; i 3n+2 fattori Uj, Vj, Pj sono funzioni di y. Per questi, dalle 
3 equazioni differenziali (due di equilibrio e una di incompressibilità) 
si ricavano 3 n + 2 equazioni differenziali ordinarie per media ponderale 
con i pesi sin jx, cos jx. In un lavoro in preparazione di Algostino e Ci¬ 
cala, si presentano in forma matriciale queste equazioni per qualsiasi 
valore di n. Il sistema differenziale risulta essere di ordine 4n, separato 
da due integrali primi dove intervengono V 0 , P 0 . Analoga forma ma¬ 
triciale si dà per le condizioni al contorno, poste in numero 2n per cia¬ 
scuno dei lati di contorno y = costante. Si è programmata l’integrazione 
numerica delle equazioni differenziali e si sono elaborati metodi di «shoot- 
ing» per la riduzione al problema di valori iniziali. Qui si presentano 
alcuni risultati. 

I primi addendi nelle espressioni (5) rappresentano una deformazio¬ 
ne uniforme che, nel piano, dà solo una tensione normale nella direzio¬ 
ne /]. Per l’incompressibilità del corpo è A 2 X= 1. Per X > 1 si ha una 
compressione che accorcia nel rapporto l/VX le dimensioni secondo /, 
e allunga nel rapporto VX le ortogonali (fig. 6). Le (6),, (6) 2 rappresen¬ 
tano una deformazione nel piano la quale è simmetrica rispetto alle ret¬ 
te x=0ex=7r dove gli spostamenti U si annullano: in questi punti, quindi, 
si annullano le tensioni tangenziali nella direzione i 2 . 

L’indagine numerica è stata sviluppata per n = 2. Si è considerato il 
campo rettangolare delimitato nello stato iniziale dai lati x = 0, x=i r, 
y= ±h/2. Sul lato y= +h/2 si è supposto applicato un carico nella di¬ 
rezione 4 che dà alla componente lagrangiana dalla tensione la distri¬ 
buzione uccisa' 13 ’: sul latoj'= -h/2 si introduce la tensione - a cos a, 
con a costante. In questa situazione antimetrica, per j = 0 si annullano 
le incognite U u V 2 , P ,: la V 0 che un moto traslatorio del corpo può mo¬ 
dificare, si sceglie nulla. A rappresentare la risposta del corpo al carico 
agente si assume il valore V x sull’orlo della componente F, avente la 
stessa lunghezza d’onda del carico, e quindi implicata nel computo del 
lavoro relativo. I calcoli numerici sono stati effettuati per X = 3, h = 4. 
Dalla nota teoria' 14 ’ risulta che in questa situazione X è inferiore al va¬ 
lore critico per il quale appare una deformazione per autosoluzione. Per¬ 
ciò nella fase iniziale il corpo reagisce elasticamente (fig. 7). Presto la 


(13) Le due componenti di tensione sull’orlo, moltiplicate per dx, rappresentano le proiezioni 
della forza suH’elemento che corrisponde nella deformazione al tratto i x dx del corpo indeforma¬ 
to. Le tensioni sono riferite al modulo elastico 2/i. 

(14) K. N. Sawyers, R.S. Rivlin, Stability of a thick elasticpiate under thrust, J. Elasticity, 1982, 
p. 101-124. 
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reazione, superato un massimo, prende a diminuire con pendenze sem¬ 
pre crescenti: successivamente il diagramma indica una diminuzione nel 
valore di ascissa (snap back), scendendo a ordinate assai basse (sino a 
a= -1.7). Questa evoluzione della deformazione non è attuabile se l’a¬ 
gente che trasmette il carico non è atto a esercitare l’azione della gran¬ 
dezza e del vero prescritti dal grafico. Peraltro l’inattuabilità di questa 
evoluzione si manifesta evidente dall’esame geometrico della configura¬ 
zione. Nella fig. 8 è rappresentata la deformazione del contorno del corpo 
nella situazione di carico minimo anzidetto. Alcune delle fibre y = costante 
realizzano la condizione U=0 ancorandosi al vincolo ideale sulla faccia 
opposta alla iniziale: il materiale trabocca dietro la parete di vincolo se 
non si vieta l’attraversamento della parete a tutti i punti del corpo. Per 
formulare questa condizione basta prendere in esame la deformata dei 
segmenti di contorno y = ± h/2, i quali precederanno ogni altro punto 
del corpo nell’accostarsi alla parete. In una prima fase basta limitarsi 
a considerare i punti estremi K(x=ir, y = h/2 e x = 0, y = - h/2). In ba¬ 
se alle espressioni (6),, (6) 2 con gli sviluppi troncati a n = 2, risulta 


(7) dP/dx=f( 1-C/ 1 +2C4). 


Nel primo tratto della linea 1 di evoluzione del processo di carico, 
il fattore in parentesi nella (7), che indichiamo con 0, si mantiene posi¬ 
tivo, decrescendo dal valore 1 : la tangente al contorno in K resta paral¬ 
lela a /, e la fibra ivi situata si accorcia, accostandosi alla lunghezza 
zero. La tensione principale in quella direzione, negativa, cresce illimi¬ 
tatamente in valore: si riscontra dunque una concentrazione di sforzi 
localizzata sui due vertici. Per 0 = 0 la direzione della tengente è inde¬ 
terminata: però, se si calcolano le quantità Ax ì =x x {k)-x ì {t: - A), 
Ax 2 = x 2 (ir) — x 2 (tt — A) per A-*0, si nota che Ax x decresce come 
(A;r 2 ) 15 : la deformata assume in K la forma di parabola con tangente 
4 e esponente 1.5, ossia con curvatura infinita. Se si prosegue nel cal¬ 
colo con le condizioni al contorno invariate, si trova che il vettore tan¬ 
gente (7) si volge nella direzione - q: avviene il traboccamento 
rappresentato nella fig. 8 e la caratteristica prosegue come è indicato dalla 
linea a tratto nella fig. 7. Se invece, raggiunta in A la situazione 0 = 0, 
si impone che questa condizione permanga mentre^ dalla parete di vin¬ 
colo sia trasmessa al corpo una forza in direzione i 2 volta a contrastare 
l’avanzata, si trova che tale azione si concentra in K e la caratteristica 
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ff(K,) risale lungo la linea 2 della fig. 7 (15) . 

Ovviamente queste conclusioni sono inficiate dal troncamento degli 
sviluppi (6). Per n > 2 è da ritenere che possa ottenersi un contatto fra 
il corpo e la parete, esteso a punti distinti da K. Tuttavia questa indagi¬ 
ne approssimata può dare utili indicazioni qualitative sul fenomeno ri¬ 
scontrato. 


Conclusioni 

Nei tre esempi accennati sono indicate questioni di interesse tecnico 
risolubili oltrepassando le situazioni critiche e portandosi avanti nel campo 
non lineare. 

Assai istruttiva risulta l’analisi della struttura con imperfezioni di for¬ 
ma. Per il guscio cilindrico completo, questo esame conferma le espe¬ 
rienze che mostrano l’irraggiungibilità del carico assiale critico indicato 
dalla teoria linearizzata: per qualsiasi deformazioni previa, l’evoluzio¬ 
ne sotto carico assiale crescente porta il punto rappresentativo ad acco¬ 
starsi al cammino peggiore, quello che porta al cedimento di scatto. Invece 
per il pannello cilindrico, comunque vincolato sui 4 lati di contorno, le 
restrizioni all’evoluzione della deformazione derivanti dai vincoli pro¬ 
lungano gli effetti della deformazione previa; se questa è di tipo simile 
a quella del percorso ideale migliore, si riesce a superare il carico critico 
indicato dalla teoria linearizzata per il pannello perfetto (16) . 

Similmente informativo risulta l’esame della struttura in presenza di 
carichi «trasversali», ossia di azioni atte a produrre le deformazioni pro¬ 
prie del percorso deviato. Per il moto piano del solido iperelastico con¬ 
siderato, che raggiunge la situazione critica per X = 3.201, un lieve carico 
trasversale basta a innescare il cedimento di scatto per la deformazione 
uniforme con X = 3 (ossia quando i lati unitari del quadrato nella dire¬ 
zione di compressione e nella ortogonale assumono i valori 1/V3 e V3, 
rispettivamente). L’esame svolto per spostamenti finiti, diede informa¬ 
zioni sul comportamento della gomma assai meno pessimistiche che quelle 


(15) Nella fig. 6 sono segnati con punteggiata i lati^= ±h/2 e con tratto continuo le relative 

deformate nello stato d’equilibrio in assenza del carico trasversale (punto B, fig. 7). Si nota che, 
nell’approssimazione (6) con n = 2, invece del carico trasversale, potrebbero applicarsi le due forze 
opposte indicate, sui lati x = 0, x= ir. , , . ., 

Nella fig. 8 (a) sono segnati gli spostamenti sull’orlo che portano dalla deformazione uniforme 
allo stato relativo a a= - 1.7. , 

(16) Si è riesumato questo studio di 40 anni fa perché indirizzato sulla via che gli attuali stru¬ 
menti di computo suggeriscono. 
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fornite dall’analisi perturbativa. Quella venne sviluppata in un’accura¬ 
ta indagine 07 ’ relativa al semispazio sottoposto a compressione su una 
direzione tangente al piano terminale ed esente da tensioni su quel pia¬ 
no. In tal caso si ha il valore critico X = 3.383, sul quale si addensano 
infinite autosoluzioni. L’analisi perturbativa indica una caduta di resi¬ 
stenza al punto critico con pendenze sempre crescenti con il numero delle 
autosoluzioni tenute in conto. Dall’esame svolto a distanza del punto 
critico appare un fenomeno singolare: la superficie del corpo, là dove 
subisce un corrugamento rientrante, viene a toccare il vincolo previsto 
(oppure l’onda contigua nel corpo indefinito). Da questo contatto si ori¬ 
gina una ripresa nella rigidezza della gomma. 

L’analisi non lineare della fune metallica rivela un fenomeno che ne 
compromette la robustezza: la formazione di un cappio. Sin che questo 
non si sia verificato il superamento del carico critico non è dannoso. Le in¬ 
dagini svolte per il campo immediatamente postcritico, indicano instabilità 
per la fune il cui accorciamento sia debolmente contrastato: se il moto dei 
terminali è controllato, l’inflessione è essenzialmente reversibile. Se invece 
si forma il cappio e conseguente attorcigliamento (kinking), la deforma¬ 
zione, irreversibile per attrito interno può dar luogo alla frattura. Si è già 
osservato che l’analisi svolta da Rosenthal si riferisce a casi in cui la fune 
arriva agli stati critici per compressione: sarebbe certamente interessante 
accertare la presenza del fenomeno in esame per stati postcritici di trazione 
inoltre estendere la ricerca ad altre situazioni di vincolo più facilmente at¬ 
tuabili in un controllo sperimentale, in particolare per la fune caricata me¬ 
diante giunti cardanici (caso per cui si sono trovati risultati più pessimistici 
di quelli indicati dalla formula di Greenhill). 

I calcoli numerici occorrenti per l’analisi non lineare hanno certamente 
mole notevole. Sia nel caso in cui sia da risolvere un sistema di equazioni 
algebriche (come s’è fatto, secondo Galerkin, per lo studio menzionato sul 
guscio cilindrico) sia quando siano da risolvere equazioni differenziali or¬ 
dinarie (come per la verga) o grossi sistemi di tali equazioni (come nel caso 
di sviluppi (6) allargati) occorre procedere lungo la traiettoria con passi più 
o meno corti. L’estrapolazione dei risultati dei passi precedenti guida nella 
previsione del nuovo passo di integrazione del sistema (shooting) che poi 
viene perfezionato con processo iterativo: ma l’automatizzazione di questi 
computi, li rende oggi possibili e anzi meno deterrrenti che i passi avanzati 
del processo perturbativo, la cui elaborazione algebrica è sede di inelut¬ 
tabili errori materiali. 


(17) P. Meijers, Stability of a rubber half-space, Report of thè Laboratory of Engineering Me- 
chanics, Technische Hogeschool, Delft, Part 1 (1983), Part 2 (1985). 
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Come è ben noto, lo studio della robotica richiede apporti di scienza 
e tecnologia da settori diversi quali la meccanica, relettromeccanica, l’e¬ 
lettronica e l’informatica. 

In più, i diversi apporti si influenzano a vicenda, per cui, ad esem¬ 
pio, una buona progettazione meccanica conduce ad una più facile so¬ 
luzione dei problemi di controllo, o, viceversa, un controllo particolar¬ 
mente potente può sopperire a deficienze di carattere strutturale che pos¬ 
sono manifestarsi in certe applicazioni. 

Proprio per questo, e per il gran numero di parametri a scelta con cui ci¬ 
mentarsi per la progettazione ottimale di un robot, il relativo studio appa¬ 
re quanto mai stimolante; ciò anche se è vero che la fabbricazione di robot 
industriali è entrata da alcuni anni in una fase di relativa stabilità. 

Non così si può dire però della ricerca, italiana ed internazionale, che 
dovrà portare alla concezione dei robot della prossima generazione. 

Inoltre, ai robot cosiddetti «industriali» si stanno aggiungendo quelli 
cosiddetti «non industriali», che non hanno un loro posto fisso di lavo¬ 
ro in officina. Di questi robot esistono sinora pochi esemplari funzio¬ 
nanti, certo non per tutte le applicazioni desiderate. 

In Italia, esiste ed è al suo primo anno di applicazione il «Progetto 
Finalizzato Robotica» del CNR, mirato essenzialmente ai robot indu¬ 
striali, oltre ad un progetto dell’ENEA, riguardante i robot mobili; inoltre, 
vi è la possibilità d’accesso a programmi di ricerca europei (cito come 
esempio l’EUREKA con un programma sui robot mobili di sorveglianza). 

Verranno qui essenzialmente descritti gli sviluppi nel campo della mec¬ 
canica dei robot, anche se importanti ricerche e sviluppi sono in corso 
nei settori della applicazione dell’intelligenza artificiale, della sensori- 
stica (particolarmente della visione), della motorizzazione e dei controlli. 

I problemi meccanici, poi, possono essere suddivisi in due categorie: 
una riguarda i problemi di cinematica e dinamica dei robot, indirizzati 
a fornire ai progettisti dei controlli gli algoritmi matematici la cui el- 
borazione è necessaria per l’esecuzione delle traiettorie; l’altra riguarda 
i problemi di progettazione delle configurazioni e dei componenti (di tra¬ 
smissione, strutturali, ecc.). 


Dip. di Meccanica, Politecnico di Torino, C.so Duca degli Abruzzi, 24, 10129 Torino, Italia. 
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Per ciò che riguarda la prima categoria dei problemi, una ricerca con¬ 
siste nel tentativo di dare una formulazione in forma chiusa ai pochi ca¬ 
si d’interesse pratico in cui la cinematica inversa è risolta ancora passo 
per passo, con minor precisione e rapidità di quanto altrimenti possibi¬ 
le. Tali casi nascono, in generale, da una distribuzione non ottimale de¬ 
gli assi dei giunti, talvolta inevitabile per necessità tecnologiche, come 
accade in certi robot di saldatura con cablatura interna. Inoltre, qual¬ 
che miglioramento può ancora essere portato ai metodi di risoluzione 
delle ambiguità dovute alla non univocità delle relazioni d’inversione. 

Problemi molto più imponenti pone la dinamica dei robot. 

La questione qui riguarda la possibilità — o meno — di tener conto 
delle azioni d’inerzia nello stabilire il valore di base della corrente da 
fornire ai motori per una certa manovra, sulla base dell’indicazione dei 
sensori che, istante per istante, definiscono la geometria dei vari elementi 
dei robot nello spazio. Tale valore viene poi corretto in funzione del¬ 
l’errore istantaneo della posizione attuale rispetto a quella desiderata. 

In questi problemi, le possibilità dell’analisi dinamica si possono ri¬ 
tenere esaurite, anche se i risultati non sono soddisfacenti. 

Infatti, la risoluzione delle equazioni del moto richiede — di regola 
— un tempo di calcolo che ne impedisce l’effettuazione in tempo reale. 
Anche qui, l’applicazione del metodo delle equazioni di d’Alembert ha 
mostrato la sua superiorità rispetto al metodo delle equazioni di Lagrange, 
elegante ma troppo costoso in termini di efficienza. Partendo da una 
prima analisi che prevedeva un numero di operazioni (addizioni e mol¬ 
tiplicazioni) per ciclo di calcolo proporzionale al numero di gradi di li¬ 
bertà per le prime, ed alla quarta potenza per le seconde, si è arrivati, 
usando sofisticati metodi ricursivi, ad un risultato che, pur meno disa¬ 
stroso del primo, prevede un rapporto circa di 1 a 10 nel numero delle 
operazioni elementari. Comunque, il numero minimo è sempre troppo 
alto per una completa integrazione nel controllo delle azioni d’inerzia. 

Ci si accontenta allora di tener conto solo di una parte di tali azioni, 
e precisamente di quelle dovute ai termini della diagonale della matrice 
d’inerzia, trascurando poi anche le azioni centrifughe e di Coriolis. 

Le azioni trascurate sono considerate come un disturbo, corretto dalla 
retroazione di posizione. 

Ciò naturalmente porta ad errori istantanei tanto più grandi quanto 
più grossi sono i termini trascurati, e richiede retroazioni con forte gua¬ 
dagno, sino al livello ammesso dalla stabilità e robustezza del controllo. 

Il problema è apparentemente insolubile; potrebbe essere risolto nella 
forma presente solo con l’introduzione di processori matematici molto più 
potenti o — come è anche possibile — con memorie molto più grandi. 
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Il miglior modo di trattare un problema insolubile è però sempre quello 
di sopprimerlo. 

Si deve tener conto del fatto che il problema dinamico diventa grave 
quando il robot è spinto ad alte velocità, perché allora i termini delle 
azioni d’inerzia normalmente trascurati divengono importanti. 

Anche l’equilibramento statico delle masse dei vari bracci, che per¬ 
mette di eliminare o ridurre i termini inerziali di accoppiamento fra i 
vari gradi di libertà, risulta inefficace alle alte velocità. 

Ora gli utilizzatori dei robot richiedono proprio macchine sempre più 
veloci per minimizzare i tempi ciclo, specie nelle operazioni di manipo¬ 
lazione e montaggio. 

Non resta quindi che ridurre le masse in gioco. 

Si entra così in pieno nella seconda categoria dei problemi meccanici 
dei robot, quelli relativi alla progettazione. 

Occorre qui fare una premessa. Sin dalla loro prima concezione, ne¬ 
gli Stati Uniti, la maggior attenzione ed il maggior tempo di studio sono 
stati dedicati al sistema di controllo. La configurazione fisica ed i com¬ 
ponenti meccinici sono stati posti in secondo piano. Poiché questa linea 
di tendenza ha avuto un seguito un pò dovunque, il robot attuale, che 
ne è il risultato, è stato paragonato, con una analogia antropomorfa, 
ad una persona d’intelligenza brillante ma fisicamente menomata. 

L’analogia è sostanziata dal paragone fra le prestazioni fisiche di un 
essere umano e di un robot. L’uomo, considerato come un manipolato- 
re, può maneggiare a velocità discretamente alta carichi sino ad una de¬ 
cina di chili, a velocità un pò più bassa carichi sino a 25 chili, e può fare 
semplici movimenti anche con carichi intorno al quintale, mentre il suo 
peso proprio si aggira intorno agli ottanta chili. 

Ciò corrisponde a maneggiare con discreta facilità carichi fra un quarto 
ed un quinto del proprio peso. 

Per contrasto, i robot industriali possono maneggiare carichi di gran 
lunga inferiori, quando rapportati al peso proprio. 

Il paragone risulta particolarmente sfavorevole ai robot per bassi carichi. 

Inoltre, il costo del sistema di controllo, dei sensori, dell’hardware 
e del software è praticamente indipendente dalle dimensioni dei robot. 

Ne segue che l’inefficienza del sistema meccanico si ripercuote sia sulle 
prestazioni che sul costo dei robot. 

Un importante campo di ricerca riguarda quindi l’alleggerimento delle 
strutture e delle altre parti mobili, come modo per ridurre l’importanza 
relativa delle azioni d’inerzia ed il costo globale. 

Si deve anzitutto esaminare la natura delle masse più cospicue, e più 
prossime all’estremità dei robot. 
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Appare subito che i contributi più importanti sono dovuti alle masse dei 
motori e riduttori, che altresì contano per una parte sostanziale, ed all’in- 
circa proporzionale alla propria massa, al costo dell’apparecchiatura, e 
secondariamente alle masse delle strutture portanti dei bracci e dal polso. 

Data l’architettura seriale secondo cui tutti gli attuali robot sono co¬ 
struiti, ogni motoriduttore deve sostenere il carico del successivo, quan¬ 
do i motori sono disposti direttamente sull’asse dei giunti; ne segue che 
se i motori potessero essere portati sulla base o vicino alla base, utiliz¬ 
zando trasmissioni lunghe, anche la struttura portante potrebbe essere 
molto alleggerita senza ridurre la rigidezza strutturale. Le trasmissioni 
(mediante sistemi articolati, alberi rotanti, catene o cinghie dentate, con 
eventuali rinvìi ed ingranaggi) dovrebbero però essere in grado di assi¬ 
curare la stessa rigidezza e gli stessi limiti di tolleranza per i giochi dati 
dai motoriduttori montati sugli assi. 

Ciò richiede componenti di alta precisione e sistemi di ripresa auto¬ 
matica del gioco e del tensionamento di cinghie o catene, che non ab¬ 
bassino il valore della frequenza minima di vibrazione della trasmissione. 
Questo richiede, ad esempio, apparecchiature di tensionamento parti¬ 
colari, capaci di interdire ogni oscillazione fra le parti. 

Le strutture portanti dei bracci possono essere alleggerite a loro volta 
usando materiali leggeri come le fibre di carbonio od il titanio: tali solu¬ 
zioni hanno però limiti di costo, e, per ciò che riguarda le fibre di car¬ 
bonio, vi è qualche problema addizionale dovuto al relativamente basso 
modulo di elasticità ed al basso valore dello smorzamento strutturale. 

Un notevole contributo alla riduzione delle masse dei motori può ve¬ 
nire anche da un buon bilanciamento dei bracci principali con attuatori 
pneumatici. Uno studio in corso punta all’estensione ad ogni possibile 
posizione dei bracci degli attuali sistemi di bilanciamento, la cui effica¬ 
cia è attualmente limitata a piccole variazioni della configurazione ed 
è solo parziale. 

I risultati raggiunti sono già soddisfacenti ma sono previsti ulteriori 
sviluppi. 

II problema dell’alleggerimento del robot risulterebbe pressoché ri¬ 
solto se esso venisse costruito con architettura parallela. Allora, tutti i 
motoriduttori potrebbero essere installati sulla base fissa. 

Gli inconvenienti che si oppongono a questa soluzione sono duplici: 
da una parte vi è una inerente mancanza di destrezza dell’organo finale 
(che potrebbe essere superata dall’adozione di una soluzione mista, 
parallelo-serie); dall’altra vi è la attuale mancanza di un sistema di con¬ 
trollo adatto, e le difficoltà di trovare un algoritmo che permetta la ri¬ 
soluzione della cinematica inversa. 
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Tutti questi problemi sono stati però superati nel caso in cui siano 
necessarie unicamente micromovimentazioni attorno a date posizioni. 

Altri problemi meccanici di configurazione si trovano quando si vuole 
aumentare la destrezza di un robot, mediante un numero sovrabbon¬ 
dante di gradi di libertà. Per questo sono stati costruiti robot flessibili, 
del tipo a tromba d’elefante o del tipo a spina dorsale, la cui scarsa rigi¬ 
dezza però non li rende adatti per operazioni richiedenti una buona pre¬ 
cisione. Sono in fase di esecuzione delle strutture snodate ad elementi 
rigidi che potrebbero essere utilizzate, ad esempio, per lavorazioni al¬ 
l’interno di recipienti attraverso un passo d’uomo. 

Alla meccanica dei robot attiene anche il campo degli organi di pre¬ 
sa, per ciò che non riguardi l’eventuale sensorizzazione mediante tra¬ 
sduttori non meccanici e l’eventuale impiego di una intelligenza locale. 

Sono state costruite mani artificiali per la presa di oggetti delle più 
varie forme e dimensioni; le mani più efficaci appaiono essere quelle co¬ 
struite sul modello umano, con quattro dita di tre falangi ciascuna, ed 
un pollice contrapposto; ciascun dito può muoversi indipendentemente 
dagli altri, chiudendo le falangi, sino a bloccarsi sotto un dato sforzo. 
Per lo studio di tali sistemi complessi implicanti la movimentazione con¬ 
temporanea di molti corpi nello spazio, siano essi mani di presa o robot 
snodati con gradi di libertà sovrabbondanti, un aiuto fondamentale è 
offerto dalle tecniche di simulazione e di grafica computerizzata, ed in 
particolare dalla modellazione solida dei vari elementi chiamati ad inte¬ 
ragire. 

Esiste infine il vasto settore applicativo dei robot mobili, che appare, 
al momento, quello destinato ai più interessanti sviluppi futuri. 

Qui, oltre a qualche risultato già consolidato, sono state effettuate 
o sono in corso ricerche sia indirizzate a finalità precise che di aspetto 
esplorativo. 

Si può notare che, mentre il robot industriale è stato essenzialmente 
ispirato al modello antropomorfo, buona parte dei robot mobili si ispi¬ 
rano ad altri modelli zoologici. 

Si vedono così macchine replicanti il passo, il trotto o la corsa di un 
cavallo; il moto di un serpente; l’andatura di un insetto (generalmente 
esapode); il salto di una cavalletta, e persino il modo di muoversi di una 
piovra aggrappata ad un solido, utilizzata, in un progetto già esecutivo, 
per la manutenzione di piattaforme petrolifere marine. 

Le attuali od eventuali applicazioni coprono un campo vastissimo, 
dalle attività spaziali a quelle edili e minerarie, con interesse per i pro¬ 
blemi della difesa civile (protezione antincendio, mezzi antiterrorismo, 
sorveglianza contro le intrusioni, ecc.), della difesa ambientale (mani- 
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polazione di sostanze tossiche o radioattive, operazioni di smantellamento 
di centrali nucleari), e persino per le spesso dimenticate attività dome¬ 
stiche. 

Sono stati, per esempio, progettati e costruiti presso il nostro labora¬ 
torio un robot a ruote e zampe per uso in ambienti strutturati, capace 
di muoversi in corridoi e salire e scendere le scale, ed un robot a cingoli 
con possibilità di assumere configurazioni da insetto, per ambienti non 
strutturati. 

Tali realizzazioni sono intese a porre in pratica un’idea iniziale, risol¬ 
vendo man mano tutti i problemi che appaiono non appena si passa da 
un progetto su carta alla fabbricazione ed alla sperimentazione dell’og¬ 
getto, in modo che la sua fattibilità venga chiaramente dimostrata. 



